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1  Executive  Summary 


1.1  Overview  and  Main  Results 

This  report  describes  work  performed  from  1  December  2007  to  15  January  2011  under  AFRL  grant  number 
FA8750-08-1-0063  entitled  “Cognitive  Airborne  Networking:  Self-Aware  Communications  via  Sensing,  Adap¬ 
tation,  and  Cross-Layer  Optimization”.  We  cover  five  general  areas:  (1)  cooperative  automatic -repeat-request 
(ARQ)  relaying  in  airborne  networks;  (2)  power  assignment  for  hybrid  ARQ  protocols  in  airborne  networks; 
(3)  differential  cooperative  relaying  in  airborne  networks;  (4)  joint  power  optimization  for  multi-source  multi¬ 
destination  relay  networks;  and  (5)  cognitive  code-division  channelization  in  airborne  networks.  Some  main 
results  are  summarized  in  the  following: 

•  A  new  analytical  methodology  was  developed  to  evaluate  the  outage  probability  of  cooperative  decode- 
and-forward  (DF)  automatic -repeat -request  (ARQ)  relaying  under  packet-rate  fading  (fast  fading  or  block 
fading)  channels.  Specifically,  (i)  we  derived  a  closed-form  asymptotically  tight  (as  SNR  — >  oo)  approx¬ 
imation  of  the  outage  probability;  (ii)  we  showed  that  the  diversity  order  of  the  DF  cooperative  ARQ 
relay  scheme  is  equal  to  2 L  —  1,  where  L  is  the  maximum  number  of  ARQ  (re)transmissions;  and  (iii) 
we  developed  an  optimum  power  allocation  for  the  DF  cooperative  ARQ  relay  scheme. 

•  We  addressed  the  fundamental  problem  of  identifying  the  optimal  power  assignment  sequence  for  hybrid 
automatic-repeat-request  (H-ARQ)  communications  over  quasi-static  Rayleigh  fading  channels.  For  any 
targeted  H-ARQ  link  outage  probability,  we  determined  the  sequence  of  power  values  that  minimizes 
the  average  total  expended  transmission  power.  We  first  derived  a  set  of  equations  that  describe  the  op¬ 
timal  transmission  power  assignment  and  enable  its  exact  recursive  calculation.  To  reduce  calculation 
complexity,  we  also  developed  an  approximation  to  the  optimal  power  sequence  that  is  close  to  the  nu¬ 
merically  calculated  exact  result.  The  newly  founded  power  allocation  solution  reveals  that  conventional 
equal-power  H-ARQ  assignment  is  far  from  optimal.  Interestingly,  we  found  that  the  optimal  transmis¬ 
sion  power  assignment  sequence  is  neither  increasing  nor  decreasing;  its  form  depends  on  given  total 
power  budget  and  targeted  outage  performance  levels. 

•  Differential  Amplify- and-Forward  (DAF)  relaying  is  an  attractive  cooperative  communication  strategy 
for  airborne  networks  where  channel  estimation  is  not  feasible  or  it  is  rather  avoided.  We  found  a  new 
exact  outage  probability  expression  for  DAF  relaying  that  involves  only  a  single  integral,  which  is  much 
simpler  than  an  existing  result  available  in  literature  that  involves  a  triple  integral.  To  get  further  insight 
understanding  of  the  DAF  relaying  scheme,  we  obtained  an  asymptotically  tight  closed-form  approx¬ 
imation  at  high  signal-to-noise  ratio  (SNR)  scenario.  Based  on  the  tight  approximation  of  the  outage 
probability,  we  were  able  to  determine  an  asymptotically  optimum  power  assignment  for  the  DAF  relay¬ 
ing  scheme. 

•  We  developed  a  low-complexity  joint  power  assignment  algorithm  for  multi-source  and  multi-destination 
relay  networks  where  multiple  sources  share  a  common  relay  that  forwards  all  received  signals  simulta¬ 
neously  to  destinations.  We  considered  two  power  optimization  strategies:  (i)  Minimization  of  the  total 
transmission  power  of  the  sources  and  the  relay  under  the  constraint  that  the  signal  to  interference  plus 
noise  ratio  (SINR)  requirement  of  each  source-destination  pair  is  satisfied,  and  (ii)  Maximization  of  the 
minimum  SINR  among  all  source-destination  pairs  subject  to  any  given  total  power  budget.  Both  opti¬ 
mization  problems  involve  K  power  variables,  where  K  is  the  number  of  source-destination  pairs  in  the 
network,  and  an  exhaustive  search  is  prohibitive  for  large  K.  In  this  project,  we  developed  a  method¬ 
ology  that  allows  us  to  obtain  an  asymptotically  tight  approximation  of  the  SINR  and  reformulate  the 
original  optimization  problems  to  single-variable  optimization  problems,  which  can  be  easily  solved  by 
numerical  search  of  the  single  variable.  Then,  the  corresponding  optimal  transmission  power  at  each 
source  and  relay  can  be  calculated  directly.  The  proposed  optimization  schemes  are  scalable  and  lead 
to  power  assignment  algorithms  that  exhibit  the  same  optimization  complexity  for  any  number  ( K )  of 
source-destination  pairs  in  the  network.  Moreover,  we  also  applied  the  methodology  that  we  developed 
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to  solve  a  related  max-min  SINR  based  optimization  problem  in  which  we  determined  power  assignment 
for  the  sources  and  the  relay  to  maximize  the  minimum  SINR  among  all  source-destination  pairs  subject 
to  any  given  total  power  budget. 

•  We  solved  the  problem  of  cognitive  code-division  channelization  for  cognitive  airborne  networks.  The 
goal  is  to  improve  the  efficiency  of  spectrum  utilization  and  coexisting  of  primary  and  secondary  users  in 
heterogeneous  cognitive  airborne  networks,  by  jointly  assigning  power  and  code-channel  allocation  for 
a  secondary  transmitter/receiver  pair  coexisting  with  a  primary  code-division  multiple-access  (CDMA) 
system.  In  pursuit  of  a  computationally  manageable  and  performance-wise  appealing  suboptimal  solu¬ 
tion,  we  first  converted  the  amplitude/code-vector  optimization  problem  to  an  equivalent  matrix  optimiza¬ 
tion  problem  under  a  rank- 1  constraint.  Disregarding  (relaxing  in  formal  language)  the  rank- 1  constraint 
makes  the  problem  amenable  to  an  easy  polynomial-cost  semidefinite  programming  solution.  When  luck¬ 
ily,  a  rank- 1  matrix  happens  to  be  returned,  optimal  secondary-line  design  is  achieved.  For  the  common 
case  of  a  higher  rank,  we  developed  an  iterative  linearized  polynomial-cost  convex  optimizer  with  much 
appealing  (yet  suboptimal)  amplitude/code-vector  design  solutions  after  a  few  iterations.  Extensive  nu¬ 
merical  studies  have  validated  our  theoretical  developments  and  the  proposed  iterative  algorithm.  The 
proposed  scheme  almost  doubled  the  occurrences  of  secondary  transmission  compared  to  some  early 
work.  The  joint  power  and  sequence  optimization  executed  by  the  proposed  scheme  results  in  supe¬ 
rior  SINR  performance  for  the  secondary  receiver  compared  to  prior  work.  The  proposed  scheme  has 
great  potential  to  improve  the  efficiency  of  coexisting  of  primary  and  secondary  users  in  heterogeneous 
cognitive  airborne  networks. 

The  introduction  to  the  project  and  an  overview  of  each  area  is  given  in  Section  2.  Then  in  Section  3  we 
present  specifically  models,  assumptions,  methods,  and  main  results  for  each  area.  Section  4  contains  the 
conclusions  and  bibliography  is  included  at  the  end. 
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2  Introduction  to  the  Project 

2.1  Overview  of  Cooperative  ARQ  Relaying 

Conventional  wireless  networks  involve  point-to-point  communication  links  and  for  that  reason  do  not  guaran¬ 
tee  reliable  transmissions  over  severe  fading  channels.  On  the  other  hand,  cooperative  wireless  networks  exhibit 
increased  network  reliability  due  to  the  fact  that  information  can  be  delivered  with  the  cooperation  of  other  users 
in  networks  [1] — [13].  In  particular,  in  cooperative  systems  each  user  utilizes  other  cooperative  users  to  create 
a  virtual  antenna  array  and  exploit  spatial  diversity  that  minimizes  the  effects  of  fading  and  improves  overall 
system  performance.  Cooperative  communications,  also  known  as  relay  channels,  was  first  introduced  in  [9] 
in  which  a  three-way  channel  was  analyzed  based  on  the  capacity  region.  More  information-theoretic  analysis 
of  the  relay  channels  was  developed  in  [10]  and  recently  in  [11,  12],  Various  practical  cooperative  communi¬ 
cation  protocols  were  proposed  for  wireless  networks  [1]— [6]  and  performance  analysis  has  been  extensively 
developed  (see  [7,  13]  and  the  references  therein). 

Automatic-repeat-request  (ARQ)  protocols  for  wireless  communications  have  been  studied  extensively  in  the 
past  and  proved  themselves  as  efficient  control  mechanisms  for  reliable  data  packet  transmissions  at  the  data 
link  layer  [15]— [21].  The  basic  idea  of  ARQ  protocols  is  that  a  receiver  requests  retransmission  when  a  packet 
is  not  correctly  received.  Recently,  in  an  effort  to  increase  network  reliability  over  poor  quality  channels,  ARQ 
protocols  were  studied  in  the  context  of  cooperative  relay  networks  [22]— [25] .  In  particular,  [22]  was  among 
the  first  such  studies  to  present  a  general  framework  of  cooperative  ARQ  relay  networks.  It  was  shown  that 
cooperative  ARQ  relay  networks  have  great  advantages  in  terms  of  throughput,  delay,  and  energy  consumption 
compared  with  the  conventional  multihop  ARQ  approach  in  which  point-to-point  ARQ  links  are  concatenated 
in  network  routes.  In  [23],  information-theoretic  analysis  was  developed  and  upper  bounds  for  the  diversity 
order  of  a  decode -and- forward  (DF)  cooperative  ARQ  relay  scheme  were  characterized  for  both  slow  and  fast 
fading  channels  as  a  means  to  study  the  diversity-multiplexing-delay  tradeoff  for  the  cooperative  ARQ  relay 
scheme.  In  [25],  a  closed-form  expression  of  the  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme 
was  obtained  for  slow  fading  channels,  but  unfortunately  the  introduced  approach  cannot  be  extended  to  fast 
fading  channels. 

Outage  probability  is  arguably  a  fundamental  performance  metric  for  wireless  ARQ  relay  schemes  and  so 
is  the  diversity  order.  In  this  project,  we  developed  a  new  analytical  methodology  for  the  treatment  of  DF 
cooperative  ARQ  relay  networks  in  fast  fading  (packet-rate  fading)  channels,  which  leads,  for  the  first  time,  to 
a  closed-form  asymptotically  tight  (as  SNR  — >  oo)  approximation  of  the  outage  probability.  The  closed-form 
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expression  shows  that  the  overall  diversity  order  of  the  DF  cooperative  ARQ  relay  scheme  is  equal  to  2 L  —  1, 
where  L  is  the  maximum  number  of  ARQ  retransmissions.  The  achieved  diversity  is  partially  due  to  the  DF 
cooperative  relaying  and  partially  due  to  the  fast  fading  nature  of  the  channels  (temporal  diversity  due  to  re¬ 
transmissions  over  independent  fading  channels).  We  note  that  the  diversity  of  the  direct  ARQ  scheme  (without 
relaying)  is  only  L  and  it  is  due  to  the  fast  fading  channels.  Based  on  the  asymptotically  tight  approximation 
of  the  outage  probability,  we  arc  able  to  determine  an  optimum  resource  allocation  for  the  DF  cooperative 
ARQ  relay  scheme.  For  any  given  total  transmission  power  budget,  optimum  power  allocated  at  the  source 
and  at  the  relay  can  be  determined  which  depends  on  the  variances  of  the  channels  and  the  maximum  number 
of  retransmissions  allowed  in  the  protocol.  It  turns  out  that  the  conventional  equal-power  allocation  strategy 
is  not  optimum  in  general  and  the  optimum  power  allocation  relies  heavily  on  the  link  quality  of  the  channels 
related  to  the  relay.  We  have  conducted  extensive  numerical  and  simulation  results  to  illustrate  and  validate  the 
theoretical  developments. 

2.2  Overview  of  Optimal  Power  Assignment  for  Hybrid-ARQ  Protocols 

ARQ  communication  protocols,  in  which  a  receiver  requests  retransmission  when  a  packet  is  not  correctly 
received,  are  commonly  used  in  data  link  control  to  enable  reliable  data  packet  transmissions  as  shown  in  [15]— 
[21],  In  a  basic/simplest  ARQ  protocol,  a  receiver  decodes  an  information  packet  based  only  on  the  received 
signal  in  each  transmission  round  [15,  16].  Advanced  ARQ  schemes,  in  which  a  receiver  may  decode  an 
information  packet  by  combining  received  signals  from  all  previous  transmission  rounds,  have  been  known 
as  hybrid  ARQ  (H-ARQ)  protocols  [17] — [21].  Since  the  receiver  needs  to  save  previously  received  signals,  H- 
ARQ  communication  protocols  require  more  memory  at  the  receiver  side  compared  to  the  basic  ARQ  protocols. 
However,  the  performance  of  the  H-ARQ  protocols  is  substantially  better  than  that  of  the  basic  ARQ  protocols 
and  the  performance  improvement  is  worth  of  the  memory  increase  at  the  receiver  side  [16,  20,  21],  especially 
in  nowadays  memory  chips  arc  cheaper  and  their  sizes  become  smaller. 

In  wireless  links  formed  by  wireless  devices  with  limited  power  resources,  power  efficiency  is  a  key  research 
matter  in  the  optimization  of  ARQ  retransmission  protocols  [31]— [35],  In  [31],  a  power  control  scheme  was 
proposed  for  ARQ  retransmissions  in  down-link  cellular  systems  in  order  to  minimize  the  total  transmission 
power  of  multiple  users  where  each  user  uses  constant  transmission  power.  In  [32],  the  power  efficiency  of  vari¬ 
ous  ARQ  protocols  was  discussed  by  taking  into  account  the  energy  consumed  by  the  transmitting  and  receiving 
electronic  circuitry  in  ARQ  retransmissions.  Note  that  in  both  [31]  and  [32],  the  power  efficiency  of  ARQ  pro¬ 
tocols  was  examined  under  the  assumption  of  the  same  transmission  power  level  in  each  retransmission  round. 
In  [33],  the  transmission  power  in  each  retransmission  round  was  optimized  for  a  variety  of  ARQ  protocols  by 
assuming  that  channel  state  information  (CSI)  is  available  at  the  transmitter  side  and  CSI  takes  values  from  a 
prescribed  finite  set  of  values.  In  [34],  by  assuming  that  partial  CSI  is  available,  optimal  transmission  power  in 
each  retransmission  round  was  determined  for  an  H-ARQ  protocol  by  a  linear  programming  method  that  selects 
a  power  value  from  a  set  of  discrete  power  levels.  Recently  in  [35],  without  assuming  CSI  available  at  the  trans¬ 
mitter  side,  an  optimal  power  transmission  strategy  was  identified  for  a  basic  ARQ  protocol  where  the  receiver 
decodes  based  only  on  the  received  signal  in  each  transmission  round.  It  was  assumed  that  the  channel  changes 
independently  in  each  retransmission  round.  A  necessary  and  sufficient  condition  for  the  optimal  transmission 
power  sequence  was  found  which  indicates  that  power  must  be  increasing  in  every  retransmission.  We  note  that 
this  result  is  not  valid  to  slowly  fading  channels.  More  recently  in  [36],  without  a  priori  CSI  at  the  transmitter, 
the  authors  maximized  the  average  transmission  rate  for  an  incremental  redundancy  H-ARQ  protocol  where  the 
transmitter  sends  out  different  encoded  redundant  parity  symbols  in  each  retransmission  round.  The  average 
transmission  rate  maximization  under  optimal  power  assignment  was  also  formulated  and  numerical  results 
were  presented  for  an  incremental  redundancy  H-ARQ  protocol  with  one  maximum  retransmission. 

In  this  project,  we  considered  advanced  H-ARQ  transmission  protocols  for  airborne  networks  in  which  a 
destination  node  may  decode  an  information  packet  by  combining  all  received  signals  from  previous  trans¬ 
mission  rounds  to  increase  detection  reliability.  We  assume  that  the  source-destination  channel  experiences 
quasi-static  Rayleigh  fading,  i.e.  the  channel  does  not  change  during  retransmissions  of  the  same  information 
packet  and  it  may  change  independently  when  transmitting  a  new  information  packet.  Our  goal  is  to  find  the 
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optimal  power  assignment  strategy  that  minimizes  the  average  total  transmission  power  for  any  given  targeted 
outage  probability.  First,  we  derived  a  set  of  equations  that  describe  the  optimal  transmission  power  values 
in  H-ARQ  retransmission  rounds.  Then,  a  simple  recursive  algorithm  was  developed  to  exactly  calculate  the 
optimal  transmission  power  level  for  each  retransmission  round.  Interestingly,  it  turns  out  that  the  optimal 
transmission  power  assignment  sequence  is  neither  increasing  nor  decreasing;  its  form  depends  on  given  total 
power  budget  and  targeted  outage  performance  levels.  This  is  fundamentally  different  from  the  case  in  [35]  that 
the  optimal  transmission  power  must  be  increasing  in  retransmissions  in  the  fast  fading  scenario  (i.e.  the  chan¬ 
nel  changes  independently  in  each  retransmission  round).  To  reduce  calculation  complexity  and  obtain  more 
insight  understanding  of  the  optimal  power  assignment  strategy,  we  also  developed  an  approximation  to  the 
optimal  power  sequence  that  is  close  to  the  numerically  calculated  exact  result.  The  tight  approximation  shows 
that  the  optimal  transmission  power  in  each  retransmission  round  is  a  function  of  P\  (the  transmission  power 
in  the  first  round)  in  a  polynomial  form.  The  optimal  power  assignment  values  also  reveal  that  the  conventional 
equal-power  assignment  (using  the  same  transmission  power  in  all  retransmission  rounds)  is  far  from  optimal. 
As  an  example,  for  a  targeted  outage  probability  of  10" 3  and  maximum  number  of  transmissions  L  =  2,  the 
average  total  transmission  power  based  on  the  optimal  power  assignment  is  9  dB  less  than  that  of  using  the  com¬ 
mon  equal-power  scheme.  We  also  observed  that  the  larger  the  maximum  number  of  retransmissions  allowed  in 
the  H-ARQ  protocol  or  the  lower  the  required  outage  probabilities,  the  more  power  savings  the  optimal  power 
assignment  strategy  offers. 

2.3  Overview  of  Differential  Amplify-and-Forward  Relaying 

The  concept  of  cooperative  relaying  has  been  attracting  much  attention  recently  with  the  promise  of  increased 
wireless  network  capacity  and  thorughput/delay  performance  ([39,  53]  and  references  therein).  In  cooperative 
relaying,  a  source  node  involves  nearby  nodes  to  jointly/cooperatively  send  information  to  a  destination.  To 
reduce  implementation  complexity  and  system  requirements,  differential  cooperative  relaying  is  particularly  at¬ 
tractive  in  wireless  networks  due  to  its  advantage  that  information  can  be  differentially  modulated/demodulated 
and  channel  estimation  can  be  avoided  (at  the  cost  of  an  arguably  tolerable  error  rate  loss).  In  [41],  Differential 
Decode- and- Forward  (DDF)  relaying  was  proposed  and  optimized  in  which  each  relay  differentially  decodes 
the  information  from  a  source  node  and  forwards  it  to  a  destination.  Since  relays  may  decode  information  in¬ 
correctly,  some  reliability  control  mechanism  has  to  be  deployed  at  relays  in  the  DDF  relaying  scheme  to  avoid 
error  propagation  from  relays  to  destination.  In  [42],  Differential  Amplify-and-Forward  (DAF)  relaying  was 
proposed  in  which  a  signal  from  a  source  is  simply  amplified  by  the  relays  and  forwarded  to  the  destination.  In 
the  DAF  relaying  scheme,  implementation  at  each  relay  is  quite  simplified  since  no  decoding  and  detection  is 
needed  in  this  case.  The  outage  performance  of  the  DAF  relay  has  been  studied  [43],  however,  the  presented 
outage  probability  result  involves  a  triple  integral,  which  therefore  offers  little  insight  in  the  problem  of  DAF 
protocol  optimization. 

In  this  project,  we  applied  and  optimized  the  DAF  relaying  scheme  for  airborne  networks.  Specifically,  we 
devised  a  new  exact  outage  probability  expression  for  DAF  relaying  which  involves  only  a  single  integral. 
Furthermore,  to  provide  further  insight  in  the  process  of  DAF  relaying,  we  developed  a  simple  closed- form  ap¬ 
proximation  of  the  outage  probability  which  is  tight  at  high  signal-to-noise  ratio  (SNR)  and  reveals  the  asymp¬ 
totic  performance  of  the  DAF  relaying  scheme.  Then,  based  on  the  tight  outage  probability  approximation, 
we  were  able  to  determine  an  asymptotically  optimum  power  allocation  scheme  for  DAF  relaying.  Numerical 
calculations  and  simulation  results  have  been  carried  out  to  validate  and  illustrate  our  theoretical  findings. 

2.4  Overview  of  Joint  Power  Optimization  for  Multi-Source  Multi-Destination  Relay  Net¬ 
works 

Recent  work  on  information-theoretic  aspects  of  cooperative  relaying  as  well  as  recent  proposals  of  practical 
cooperative  relaying  protocols  ([50]— [57]  and  references  therein)  suggest  that  cooperative  relaying  may  lead  to 
significant  improvements  in  detection  reliability  at  destinations  and  overall  system  performance.  In  cooperative 
relaying,  a  user/node  may  serve  as  a  relay  and  assist  others  by  forwarding  their  signals  to  destinations,  thus 
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enhancing  detection  reliability  at  the  destination.  Various  relaying  strategies  have  been  studied  in  the  literature 
for  relay  to  forward  signals.  For  example,  relay  may  decode  the  received  signal  and  forward  the  decoded 
information  to  a  destination,  or  it  may  simply  amplify  the  received  signal  and  forward  it  to  the  destination. 

More  recently,  there  is  increasing  interest  in  investigating  the  advantages  of  relaying  in  multi-source  multi¬ 
destination  networks  [58]— [65],  [71] — [73].  The  simplest  multi-source  multi-destination  relay  network  is  mod¬ 
eled  as  an  interference  relay  channel  (IRC)  [58]  where  a  relay  helps  two  independent  source-destination  pairs  by 
using  different  relaying  strategies  such  as  dc  code  -  an  d-  fo  rw  aid .  amplify-and-forward,  or  comprcss-and- forward. 
Past  literature  on  multi-source  multi-destination  relay  networks  focused  primarily  on  information-theoretic 
studies  including  achievable  rate  regions  or  bounds  of  capacity  region  [58]— [65] .  For  example,  in  [58]  a  rate 
splitting  technique  is  used  to  study  the  problem  of  achievable  rate  region  for  a  Gaussian  IRC  channel,  where 
each  message  is  split  into  a  common  message  which  is  decodable  at  all  destinations  and  a  private  message  which 
is  decodable  only  at  the  intended  destination  [61].  Since  the  receivers  arc  able  to  decode  paid  of  the  interference 
messages,  the  effect  of  interference  is  reduced  and  the  overall  communication  rate  is  therefore  increased.  The 
achievable  rate  region  of  [58]  was  further  improved  in  [59]  by  considering  both  intended  message  and  inter¬ 
ference  forwarding  at  the  relay  (optimal  relay  strategies  were  studied  under  the  assumption  that  the  relay  is 
connected  to  each  source  and  each  destination  via  orthogonal  and  finite  capacity  links).  The  capacity  region  of 
the  interference  channel  with  a  single -relay  was  investigated  in  [62,  63],  where  it  was  shown  that  forwarding 
the  intended  message  of  just  one  source,  the  achievable  rates  for  both  source-destination  pairs  can  be  improved. 
By  assuming  that  the  relay  knows  the  source  message  a  priori,  a  relaying  strategy  was  proposed  in  [64,  65] 
where  generalized  beamforming  with  dirty  paper  coding  was  considered  for  a  two-source  two-destination  relay 
network  to  further  improve  the  capacity  region  of  the  network. 

Power  control  is  an  important  technique  to  improve  the  performance  of  multi-source  multi-destination  re¬ 
lay  networks.  While  past  literature  offers  a  significant  amount  of  work  on  power  allocation  for  single-source 
single-destination  relay  networks  (see,  for  example  [67]-[70]  and  references  therein),  there  arc  rather  limited 
studies  on  power  optimization  for  multi-source  multi-destination  relay  networks.  In  [71],  power  allocation  was 
optimized  by  exhaustive  search  for  a  two-source  two-destination  relay  network  where  a  half-duplex  decode- 
and-forward  relay  was  considered.  Unfortunately  the  exhaustive  search  is  not  scalable  and  leads  to  prohibitive 
optimization  complexity  for  networks  with  larger  number  of  source-destination  pairs.  In  [72,  73],  a  power 
allocation  scheme  was  proposed  for  a  multi-source  multi-destination  relay  network  based  on  geometric  pro¬ 
gramming  (the  scheme  assumes  that  signals  from  different  sources  are  sent  through  orthogonal  channels  and 
the  direct  transmission  link  is  not  involved  in  detection  at  destinations). 

In  this  project,  we  designed  and  optimized  a  general  multi-source  multi-destination  relay  scheme  with  K 
sources  ( K  can  be  large)  for  cognitive  airborne  networks.  The  network  allows  simultaneous  multi-source  trans¬ 
missions  through  non-orthogonal,  in  general,  channels.  We  considered  two  power  optimization  strategies:  (i) 
Minimization  of  the  total  power  consumption  of  all  sources  and  relay  under  the  constraint  that  the  signal  to 
interference  plus  noise  ratio  (SINR)  requirement  of  each  source-destination  pair  is  satisfied,  and  (ii)  Maximiza¬ 
tion  of  the  minimum  SINR  among  all  source-destination  pairs  subject  to  any  given  total  power  budget.  Thanks 
to  an  asymptotically  tight  approximation  of  the  SINR  that  we  developed,  we  were  able  to  reformulate  the  orig¬ 
inal  optimization  problems,  which  involve  K  power  variables,  to  single-variable  optimization  problems.  Then, 
the  resulting  optimization  problems  can  be  easily  solved  by  a  simple  numerical  search  of  the  single  variable. 
The  proposed  optimization  schemes  are  scalable  and  lead  to  power  assignment  algorithms  that  exhibit  the  same 
optimization  complexity  for  any  number  (K)  of  source-destination  pairs  in  the  network.  Moreover,  for  the 
special  case  of  transmission  over  orthogonal  channels,  we  were  able  to  further  simplify  the  single-variable  op¬ 
timizations  and  obtain  analytical  solutions  for  a  symmetric  system.  Extensive  numerical  studies  show  that  the 
proposed  power  assignment  is  almost  identical  to  the  exhaustive  search  method,  and  the  optimum  power  as¬ 
signment  schemes  can  significantly  improve  the  performance  of  multi-source  multi-destination  relay  networks 
with  equal  power  assignment. 
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2.5  Overview  of  Cognitive  Code-Division  Channelization 


Recent  experimental  studies  [77]  demonstrated  that  much  of  the  licensed  radio  spectrum  experiences  low  uti¬ 
lization.  Cognitive  radio  (CR)  [78]  emerges  as  a  promising  technology  to  improve  spectrum  utilization  by 
allowing  secondary  users/networks  to  share  spectrum  licensed  by  primary  users.  As  licensees,  the  primary 
users  arc  to  have  guaranteed  access  to  the  spectrum  [79].  Therefore,  the  underlying  challenge  of  CR  technol¬ 
ogy  is  to  ensure  the  Quality-of-Service  (QoS)  requirements  of  the  primary  users  and,  simultaneously,  maximize 
in  a  best-effort  context  the  QoS  of  the  secondary  users  [80]-[83]. 

Herein,  we  arc  particularly  interested  in  cognitive  radio  built  around  a  code-division  multiple-access  (CDMA) 
primary  system.  In  contrast  to  frequency  or  time  division  operation  where  cognitive  secondary  users  may 
transmit  opportunistically  in  sensed  spectrum  holes/void  only,  cognitive  code-division  users  may  in  principle 
operate  in  parallel  in  frequency  and  time  to  a  primary  system  as  long  as  the  induced  spread-spectrum  inter¬ 
ference  remains  below  a  pre-defined  acceptable  threshold1.  Power  control  for  cognitive  code-division  systems 
was  considered  in  [85]  under  an  “interference  temperature”  constraint  (total  secondary  user  disturbance  power 
over  primary  band).  No  optimization  was  carried  out  with  respect  to  the  code  channels  (signatures)  of  the 
secondary  users  in  [85].  In  contrast,  in  [86]  a  secondary  code  assignment  scheme  was  presented  to  minimize 
the  mean-square  crosscorrelation  of  the  secondary  code  with  the  primary  received  signal.  Extension  to  multiple 
secondary  users  was  also  considered  in  the  form  of  iterative  secondary  code  set  construction.  Under  interference 
minimizing  code  assignments,  bit  rate  and  spreading  factor  adjustments  for  a  secondary  CDMA  system  were 
considered  in  [87].  Interesting  work  outside  the  framework  of  CDMA  CR  in  the  form  of  joint  beamforming 
and  power  allocation  algorithms  was  reported  in  [88],  [89],  while  auction  mechanisms  for  power  control  were 
presented  in  [90]. 

In  this  project,  we  considered  the  problem  of  cognitive  code-division  channelization  for  cognitive  airborne 
networks.  The  goal  is  to  improve  the  efficiency  of  spectrum  utilization  and  coexisting  of  primary  and  secondary 
users  in  heterogeneous  cognitive  airborne  networks,  by  jointly  assigning  power  and  code-channel  allocation  for 
a  secondary  transmitter/receiver  pair  coexisting  with  a  primary  code-division  multiple-access  (CDMA)  sys¬ 
tem.  In  pursuit  of  a  computationally  manageable  and  performance-wise  appealing  suboptimal  solution,  we 
first  converted  the  amplitude/code-vector  optimization  problem  to  an  equivalent  matrix  optimization  problem 
under  a  rank- 1  constraint.  Disregarding  (relaxing  in  formal  language)  the  rank- 1  constraint  makes  the  prob¬ 
lem  amenable  to  an  easy  polynomial-cost  semidefinite  programming  solution.  When  luckily,  a  rank- 1  matrix 
happens  to  be  returned,  optimal  secondary-line  design  is  achieved.  For  the  common  case  of  a  higher  rank, 
we  developed  an  iterative  linearized  polynomial-cost  convex  optimizer  with  much  appealing  (yet  suboptimal) 
amplitude/code-vector  design  solutions  after  a  few  iterations.  Extensive  numerical  studies  have  validated  our 
theoretical  developments  and  the  proposed  iterative  algorithm.  The  proposed  scheme  almost  doubled  the  oc¬ 
currences  of  secondary  transmission  compared  to  some  early  work.  The  joint  power  and  sequence  optimization 
executed  by  the  proposed  scheme  results  in  superior  SINR  performance  for  the  secondary  receiver  compared 
to  prior  work.  The  proposed  scheme  has  great  potential  to  improve  the  efficiency  of  coexisting  of  primary  and 
secondary  users  in  heterogeneous  cognitive  airborne  networks. 


3  Models,  Assumptions,  Methods  and  Procedures 

3.1  Cooperative  Decode-and-Forward  ARQ  Relaying 

In  this  subsection,  first  we  describe  briefly  the  DF  cooperative  ARQ  relay  scheme  and  model  the  fast  (packet- 
rate)  fading  channels.  Second,  we  develop  two  useful  lemmas  which  arc  the  bases  of  our  analytical  approach, 
and  calculate  the  outage  probability  and  derive  the  asymptotically  tight  outage  probability  for  the  DF  cooper¬ 
ative  ARQ  relay  scheme.  Third,  based  on  the  tight  approximation  of  the  outage  probability,  we  determine  the 
asymptotically  optimum  power  allocation  for  the  DF  cooperative  ARQ  relay  scheme.  Numerical  and  simulation 
studies  arc  presented  and  some  conclusions  arc  drawn  at  the  end  of  this  subsection. 

'While  early  standardization  and  regulation  discussions  have  begun  [84],  no  conclusive  “interference  temperature”  rules  and  agree¬ 
ments  have  been  reached  yet. 


Relay 


Source 


Destination 


Figure  1:  Illustration  of  the  cooperative  ARQ  relay  scheme  with  one  source,  one  relay  and  one  destination. 


3.1.1  System  Model 

We  consider  a  cooperative  ARQ  relay  scheme  with  one  source,  one  relay  and  one  destination  as  illustrated  in 
Fig.  1.  The  DF  cooperative  ARQ  relay  scheme  works  as  follows.  First,  a  data  packet  of  b  bits  are  encoded 
into  a  codeword  of  length  AT,  where  L  is  the  maximum  number  of  ARQ  retransmission  rounds  allowed  in  the 
protocol  and  T  is  the  number  of  channel  uses  in  a  single  ARQ  retransmission  round.  Then,  the  codeword  is 
divided  into  L  different  blocks  each  of  length  T.  In  each  ARQ  retransmission  round  a  block  of  the  message  is 
sent,  so  the  transmission  rate  is  R  =  b/T.  When  the  source  transmits  a  block  of  the  message  to  the  destination, 
it  is  also  received  by  the  relay.  The  destination  indicates  success  or  failure  of  receiving  the  message  by  feeding 
back  a  single  bit  of  acknowledgement  (ACK)  or  negtive-acknowledgement  (NACK).  The  feedback  is  assumed 
to  be  detected  reliably  at  the  source  and  at  the  relay.  If  an  ACK  is  received  or  the  retransmission  reaches  the 
maximum  number  of  rounds,  the  source  stops  transmitting  the  current  message  and  starts  transmitting  a  new 
message.  If  a  NACK  is  received  and  the  retransmission  has  not  reached  the  maximum  number  of  rounds,  the 
source  sends  another  block  of  the  same  message.  If  the  relay  decodes  successfully  before  the  destination  is 
able  to,  the  relay  starts  cooperating  with  the  source  by  transmitting  corresponding  blocks  of  the  message  to 
the  destination  by  using  a  space-time  transmission  [23],  for  example,  using  the  Alamouti  scheme  [28].  The 
destination  combines  the  received  signal  in  current  round  and  those  in  previous  rounds  to  jointly  decode  the 
data  packet.  After  L  ARQ  retransmission  rounds,  if  the  destination  still  cannot  decode  the  data  packet,  an 
outage  is  declared  which  means  that  the  mutual  information  of  the  DF  cooperative  ARQ  relay  channel  is  below 
the  transmission  rate. 

The  DF  cooperative  ARQ  relay  scheme  can  be  modeled  as  follows.  The  received  signal  yr.m  at  the  relay  at 
the  m-th  (1  <  m  <  L)  ARQ  retransmission  round  can  be  modeled  as 

yr,m  =  \J Pshsr,m%s  T  (1) 

where  Ps  is  the  transmitted  power  of  the  source  signal  xs,  hsr,m  is  the  coefficient  of  the  source-relay  channel  at 
the  m-th  ARQ  retransmission  round,  and  ryvm  is  the  additive  noise.  If  the  relay  is  not  involved  in  forwarding, 
the  received  signal  y,/rr,  at  the  destination  at  the  m-th  ARQ  retransmission  round  is 

Ud,m  =  \J Pshsd,m%s  ~F  Vd,rm  (2) 

where  hs^m  is  the  source-destination  channel  coefficient  at  the  m-th  ARQ  retransmission  round.  If  the  relay 
receives  the  data  packet  from  the  source  successfully,  it  helps  in  forwarding  the  packet  to  the  destination  using 
the  Alamouti  scheme.  Specifically,  each  block  of  the  data  packet  can  be  partitioned  into  two  parts  as  xs  = 
[x\s.  i  xs^\,  then  the  relay  forwards  a  corresponding  block  xr  =  [  —x*  2  x*  -J.  The  received  signal  y,/,m  at  the 
destination  at  the  m-th  ARQ  retransmission  round  can  be  written  as 

Vd,m  =  \/ Ps hsd,m%s  T  \/ Prhrd,m%r  T  (3) 
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where  Pr  is  the  transmitted  power  at  the  relay  and  hrd.m  is  the  channel  coefficient  from  the  relay  to  the  desti¬ 
nation  at  the  m-th  ARQ  retransmission  round.  At  the  destination,  the  message  block  xs  can  be  recovered  based 
on  the  orthogonal  structure  of  the  Alamouti  code  [23,  28].  The  channel  coefficients  hsd,m ,  hsr,m  and  hrd,m 
are  modeled  as  independent,  zero-mean  complex  Gaussian  random  variables  with  variances  cr2sd ,  af.  and  cr2d, 
respectively.  We  consider  a  fast  fading  scenario,  i.e.  the  channels  remain  fixed  within  one  ARQ  retransmission 
round,  but  change  independently  from  one  round  to  another  (packet-rate  fading).  The  channel  state  information 
is  assumed  to  be  known  at  the  receiver  and  unknown  at  the  transmitter.  The  noise  7/r  m  and  r}d.m  are  modeled 
as  zero-mean  complex  Gaussian  random  variables  with  variance  Mq. 

3.1.2  Outage  Probability  Analysis 

In  this  subsection,  first  we  develop  two  lemmas  which  paly  key  roles  in  our  outage  probability  analysis.  Sec¬ 
ond,  we  derive  the  outage  probability  of  the  direct  ARQ  transmission  scheme.  Finally,  we  derive  the  outage 
probability  of  the  DF  cooperative  ARQ  relay  scheme  under  packet-rate  fading  conditions. 

First  we  develop  two  lemmas  and  some  notations.  The  two  lemmas  will  paly  key  roles  in  analyzing  the 
outage  probability  of  the  DF  cooperative  ARQ  relay  scheme. 

Lemma  1  and  uaiv..)SM  are  two  independent  random  variables  satisfying  the  following  properties 


M 

lim 

Si — >00 

n 

■  Pr[us  ir 

■;SM  <t]  =  a 

■m 

1  <i<M 

i=l 

M 

lim 

Si — >oo 

n 

42 

•  Pr[vSu.. 

:,SM  <  t]  =  b  ■ 

9(t), 

l<i<M 

i=  1 

where  d\ ,  d.2,  a  and  b  are  constants,  f(t )  and  g(t )  are  monotonically  increasing  functions,  and  f'[t )  is  inte- 
grable,  then 

M  ft 

II  si1+d2  '  Pr  ,-,»m  +  ®«i . <t]=ab-  /  g{x)f'(t  -  x)dx.  (4) 

Proof  :  For  any  partition  of  the  interval  [0,  t],  denoted  as  U  =  {tto,  u\, ...,  uj}  with  uo  =  0  and  uj  =  t, 
we  can  obtain  upper  and  lower  bounds  of  the  event  {wSl,...,SM  +  vSl,...,sM  <  t}  as  follows: 

J 

{US1,...,SM  +  VS1,...,SM  A  0  —  1  —  US1,...,SM  S  Uj}  3  {Ai,...,Sm  A  t  —  Uj—  l}, 

3= 1 
J 

{US1,...,SM  +  VS1,...,SM  A  t}  D  J  { Z4 J — 1  A  Uslt...JsM  <  Uj}  n  {'Us1,...IsM  <  t  —  Uj}. 

7=1 

The  upper  and  lower  bounds  arc  considered  as  a  union  of  rectangles  between  Uj-\  and  uj  for  0  <  j  <  J.  First, 
let  us  focus  on  the  upper  bound.  The  probability  of  the  subset  {uj-i  <  <  'Uj}n{vSlr..)SM  <  t—Uj- 1} 

can  be  calculated  as 


Pr  [uj-i  <  us1,...,sM  <  Ui,vsl,...,SM  <  t  -  Uj- 1] 


=  {Pr  [ws1)...,,m  <  ui\  ~  Pr  [ 


u 


<  Uj- i]}Pr  [vai)...,SM  <  t  -  Uj- 1]  . 


(5) 


Let  Si  — >  oc  for  1  <  i  <  M,  according  to  the  assumption  in  the  lemma,  we  have 


M 


gjim,  II  si1+d2  '  Pr  [^i-1  -  U»u-,SM  <  uPvsi,-,sM  <  t  -  uj-i]  =  ab  •  {/(«j)  -  f{uj-i)}g{t  ~  Uj- 1). 

(6) 


l<i<Mi=l 
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Since  f(t )  and  g(t)  arc  monotonically  increasing  functions,  we  obtain  an  upper  bound  of  (4)  as  follows 

M  J 

SUP u  aJ™o  II  si1+d2  '  Pr  +  vsi,..,SM  <t\<ab-J2g{t  -  «i-i){/(«i)  -  /(«j-i)}-  (7) 

l<i<M  i=l  j=l 

Similarly,  we  can  obtain  a  lower  bound  as 

M  L 

infc/  aj™,  II  si1+d2  '  Pr  . ««  +  <  *]  >  ab  ■  g(t  -  Uj){f(uj)  -  f{uj- 1)}.  (8) 

KKM  j=l  j= 1 

The  above  upper  and  lower  bounds  are  good  for  any  partition  U  =  {uq,  u\  ,  ...,*tj}  over  the  interval  [0,  t\. 
Since  f'(t)  is  integrable,  so  let  J  ->  oc,  the  the  summations  in  (7)  and  (8)  converge  to  the  same  integral 
Jo  9(x)f'(t  ~  x)dx.  Therefore  we  have  the  result  in  (4).  □ 

We  note  that  the  special  case  of  Lemma  1  with  M  =  1  was  presented  in  [29] .  Lemma  1  will  be  used  to 
approximate  the  outage  probability  of  the  ARQ  schemes  at  high  SNR  scenario.  In  the  following,  we  would  like 
to  define  a  special  function  Fn((3 1, ....  Pn\  t)  which  will  be  used  to  characterized  the  outage  probability  of  the 
DF  cooperative  ARQ  relay  scheme.  For  any  integer  n  >  2  and  non-zero  constants  P\ .  fa,  •  •  •  ,  Pn,  define 

rt  fXn  rx  2 

Fn(P  1, Pn,  t)  =  /  /  ■  2^+^+-+^dxldx2  •  •  •  dxn . 

7o  Jo  Jo 

A  closed-form  expression  for  calculating  the  special  function  Fn(fa, ....  pn:  t )  can  be  obtained.  We  have  the 
following  result. 


Lemma  2  For  any  integer  n  >  2  and  non-zero  constants  fa,  fa,  •••  ,  fa,  the  function  Fn(fa, fa;  t)  can  be 
calculated  as  follows 


Fn{fa,  fa’,t)  —  y  ' 

G{0,1} 


(_l)n+5i+-+5„_i^n2)-r 

.n:Ui[E£iW(«)ft] 


(9) 


where  the  variables  5i,  62,  ■■■,  5n- 1  £  {0, 1},  8  =  {5i,  82, ...,  5n_i},  the  coefficients  { *m,/(5 )  :  1  <  m  < 
n,  1  <  l  <  rn\  are  specified  as 


*1,1(5)  =  *2,2(5)  =  •  •  •  =  *n,n(5)  =  1, 


and,  for  any  m  =  2, 3,  .fan, 


—  5m_i  •  im—\  i (5),  /  —  1,2,..., *7i .  1 . 

Proof  :  We  use  induction  to  prove  the  result  for  any  integer  n  >  2.  When  n  =  2,  it  is  easy  to  see  that 
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2^Xl+^dxidx2 
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(>n2)  1  r2oj.+ft)<  _  p  _  (ln2>  1  fa 
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-  1 


(10) 
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i.e.  the  closed-form  expression  in  (9)  is  valid  for  n  =  2.  Next,  we  assume  that  the  result  in  (9)  is  good  for  any 
k  >2.  Then,  for  n  =  k  +  1, 


Fk+i(Pi,—,Pk+i',t)  = 


t  rxk+i  rx 2 


1 0  Jo 

ft  rxk+  i  rx  2 


0  JO 


fX  2 

/  2l3lxl+P2X2~{  h/3fc+ia;fe+1dx’idx2  •  •  •  dx’fc+i 

Jo 

rx  2 

/  2/3ia:i+/32a;2+"'+/3fca;fc  dxidx’2  •  •  •  dxfc  2/3fc+lXfe+1  dxfc+i 

Jo 


Fk((h  vj/^fc  5  *£fc+i ) 


=  f  Fk(Pi,-,Pk',  xk+l)2h+lXk+1  <\xk+1. 

Jo 

According  to  the  induction  assumption  that  the  result  in  (9)  is  good  for  any  k  >  2,  we  have 

fXk+l  fXk  fX  2 

Ffc(/3i,...,/?fc;a*+i)  =  /  /  •  •  •/  2^1+^2+-+^dxidx2  •  •  •  dxfc 

do  do  do 


(11) 


_1  - |-5fe-i  (}n2)_fc 


(-1) 


6(0,1} 


^2[Xu=i  *fc,i(^)A]®fc+i  _  _  (i2) 


Substituting  (12)  into  (1 1),  we  have 


Fk+ l(Pl,  —,Pk+ 1;  t) 

(_l)fe+<5i+-+4-i(in2)-fe  r* 

~  s,h-,  nL,  WiukJm  Jo 
6(0,1) 

- t-<5fc — i  (}n2)  — (k^1) 


v  _ 

5i,...,4_i  Um=l  [X)/=l*m,i(^)A]  •  [  X^=l  4,i(4A  +  Jk+ 1 

6(0,1} 

^ _ 2^fe+<5i-t — t-<5fc-i  (ln2)_(fc+1l 


^2lX*=i  *fc,i(^)A+/3fc+i]®fc+i  _  2^fc+l3:fe+1^dxfe+i 

2 [5D f=i  Ai;(<5)/3i+/3fc+i]t  _  ^ 


,..X,  iiL.[iz.m«aF/w 

6(0,1) 


fc+l1  _  ^ 


=  V-  (-l)fc+1+^+ -+^(ln2)-^+1)  n 

IKS.  [2Z.WWA]  1  ’’ 

6(0,1} 

where  Sk  E  {0, 1},  ik+i,k+i(S)  =  1  and  ik+ij(S)  =  Sk  ■  ik,i(8),l  =  1, 2, k.  Also,  we  have 


(13) 


k+ 1 

y^4+i,K^)A 

/=i 


Zf=i  h,i(s)Pi  +  4+i,  if  4  =  l; 

Pk+ 1)  if  4  =  o. 


(14) 


Therefore,  the  closed-form  expression  in  (9)  is  valid  for  n  =  k  +  1.  By  induction,  we  conclude  that  the  result 
in  (9)  is  true  for  all  n  >  2.  □ 

For  comparison  purpose  we  specify  the  outage  probability  for  the  direct  ARQ  transmission  scheme.  In 
the  direct  ARQ  transmission  scheme,  the  destination  receives  information  from  the  source  directly,  without 
involving  the  relay.  The  mutual  information  between  the  source  and  the  destination  in  the  m-th  round  of  the 
direct  ARQ  transmission  scheme  is 


ps 

Isd,m.  lo§2  (id-  ,  -  |d.s(/., 

Ao 


(15) 
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The  total  mutual  information  after  L  ARQ  rounds  is  Ijf  =  YJm= t  -Tsd,m •  Thus,  the  outage  probability  of  the 
direct  ARQ  scheme  after  L  ARQ  rounds  is 


P 


out,L  _ 


=  Pr  [1%  <  R  ] 


(16) 


A  closed-form  expression  of  (16)  is  not  tractable.  An  approximation  of  the  outage  probability  can  be  obtained 
for  high-SNR  scenario  as  follows  [29] 


P 


>out,L 


9l(R) 


•A/q 


where  gji.-)  is  defined  as 


9n(t)  =  /  gn-i(x)f'(t  -  x)dx,  n  >  1, 

J  o 


(17) 


(18) 


with  go(t)  =  1  and  /(f)  =  2*  —  1.  Note  that  (17)  is  a  direct  result  of  Lemma  1  with  M  =  1. 

The  calculation  of  the  coefficient  pl(7?)  involves  L  recursive  integrals.  In  the  following,  we  develop  a  closed- 
form  expression  for  the  function  gn(t)  for  any  n  >  1.  Since  gn(t)  =  gn-\(x)  f  (f  —  x)dx  and  /'(f)  =  2*ln2, 

we  have 


P%n  —  1  r%2 

9n{t)  =  /  /  •••/  gi(xi)f'(x2  -  Xi)f'(x3  -  X2)  ■  ■  -/'(Xn-i  -  Xn_2) 

7o7o  7o 

/'(f  -  xn-i)  dx\dx2  ■  •  •  dx;n_ i 
rx2 

■  (2X1  -  l)(ln2)n_12:E2"Xl2a:3“a:2  •  •  •  2a:"-1 

Jo 

t  P'Pn—  1  P-^2 

■  (1  -  2_Xl)dxidx’2  •  •  •  dxn_i 

7o 


J'  P'P'n  —  l  /'3'2 


-l-xn-22t-Xn-ldXldX2  ■  ■  ■  dxn_l 


I0J0 


=  2f(ln2  )r 


=  2 


(ln2) 


oJo 

n—  1 


(n-2) 


y  [ xn_2(l  —  2~t+x)  dx. 
’  7o 


Since  f*  xn  2  dx  =  ^yfn  1  and  ([30]) 


n—1 


xn~22x  dx  =  —2t  ^ 


(n-2)! 


_ j/n—m—l  ,  2)! 

lo  (-ln2)m(n  -  m  -  1)!  (-htf)"-1 ' 

therefore,  a  closed-form  expression  of  gn(t)  can  be  obtained  as  follows 

;(f  •ln2)n-1 


n—1 


9n(t)  —  2 


2‘E 


(n  —  1)! 
(-1)’ 


t  ^  (-l)m(f  •ln2)n-m-1  ,  Nr 

+  2'E  ~  ~  +  (-1) 


m=l 


(n  —  m  —  1)! 


^  (m  —  1)! 


(f  •ln2)m"1  +  (-l)n, 


(19) 


which  can  be  calculated  efficiently. 

In  the  following,  we  derive  the  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme  under  packet-rate 
fading  conditions.  In  the  DF  cooperative  ARQ  relay  scheme,  if  the  relay  decodes  the  message  from  the  source 
correctly,  say,  at  the  A:-th  round,  then  at  the  ( k  +  l)-th  round,  the  relay  stalls  forwarding  appropriate  ARQ 
blocks  to  the  destination.  Let  {Tr  =  k}  denote  the  event  of  successful  message  decoding  by  the  relay  at  the 
A’-th  round  and  subsequent  ARQ  block  forwarding  at  the  (A;  +  l)-th  round.  Let  P™Lk  denote  the  probability  that 
the  destination  decodes  the  message  unsuccessfully  after  L  ARQ  retransmission  rounds  if  the  event  {Tr  =  k} 
occurs.  Then,  the  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme  after  L  ARQ  retransmission 
rounds  can  be  written  as 

L 

pout,L  =  po^lk  .  pr  [Tr  =  k]  .  (20) 

k= t 
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Note  that  the  channels  change  independently  over  each  ARQ  retransmission  round  in  a  fast  fading  scenario,  so 
the  mutual  information  of  fading  channels  can  be  viewed  as  a  sum  of  independent  random  variables. 

First,  we  calculate  the  probability  of  the  event  {Tr  =  k }  ,  i.e.  Pr  [Tr  =  k\.  We  note  that  the  mutual  informa¬ 
tion  between  the  source  and  the  relay  in  the  m-th  ARQ  round  is 


Isr,m  —  lo§2 


Ps_ 

Wo 


I  h 


sr,m 


(21) 


The  probability  that  the  relay  decodes  the  message  successfully  at  the  first  round  (Tr  =  1)  is 

Pr  [Tr  =  1]  =  Pr  [ISrA  >  R]=  exp  ( -2  2  •  ^ 


ct2  P 

°  sr  rs 


For  any  Tr  =  k,  k  =  2, 3, L  —  1,  we  have 


(22) 


Pr  [Tr  =  /,;]  =  Pi¬ 


le-  1 


^  '  Isr,m  ^  Ri  ^  '  Isr,m.  E  R 


\_m.=  l 
k- 1 


=  Pr 


^  ^  Isr,m  ^  R 


m=l 


-Pr 


\_m=l 


^  ^  Isr,m  ^  R 


gk-i(R) 


CT2  P 
u  srr  S 


[_m=  1 

-  9k(R) 


A/q  A'-1  _x  (  ATo 


CT2  P 

°  srrs 


(23) 


where  gk-i(-)  and  fjk(-)  are  specified  in  (19)  in  the  previous  subsection.  The  approximation  in  (23)  is  obtained 
by  applying  Lemma  1  with  M  =  1.  Finally,  if  Tr  =  L,  we  have 


Pr  [Tr 


L\=  Pr 


'  L—l 

E 

m=  1 


<  R 


'9l-\{R) 


L- 1 


(24) 


Next,  we  calculate  the  conditional  outage  probability  P™Lk  when  the  relay  decodes  correctly  at  the  fc-th 
round  and  stalls  forwarding  at  the  ( k  +  1 ) -th  round.  When  the  relay  cooperates  with  the  source  by  jointly 
sending  a  message  block  via  the  Alamouti  scheme,  the  mutual  information  of  the  cooperative  channels  in  the 
m-th  ARQ  round  is  given  by  [23] 


Is rd,m  —  lo§2  ( 

Thus,  with  L  ARQ  rounds,  the  total  mutual  information  is 


V1  + 


2  4- 
m  | 


Pr 

A/n 


I  hrd,r 


r 


tot 

d,Tr=k 


_  J  Em=l  Isd,m  +  Sm=fc+1  Isrd^rri')  1  ^  k  L\ 

I  Em=l  Isd,mi  k  =  L. 


(25) 


(26) 


We  note  that  if  Tr  =  L,  the  relay  has  no  chance  to  cooperate  since  the  source  starts  sending  a  new  packet.  The 
conditional  outage  probability  can  be  evaluated  as 


P^ik  =  Pr  [I^r=k  <  R]  . 


(27) 


When  Tr  =  L,  the  conditional  outage  probability  is  reduced  to  the  direct  ARQ  scenario  and  it  is  given  by 


PtOUt  _ 

Tr=L  ~ 


Pr  [Id°h=i 


(28) 


In  the  following,  we  calculate  the  conditional  outage  probability  (27)  for  any  Tr  =  k,  k  =  1,2, ...,  L  —  l.  For 
simplicity  in  presentation,  we  introduce  the  following  notation 


=  flog2  (1  +  St|Wd,m|2)  , 

\log2  (l  +  Si\hsd,m\2  +  S2\hrd,m\2) 


1  <  m  <  k; 
k  +  1  <  m  <  L, 


(29) 
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(30) 


where  si  =  Psj A/q  and  S2  =  Pt/Mq.  Then,  the  total  mutual  information  can  be  written  as 


jt.ot  _ 
1d,Tr=k  ~  /  , 


Urn- 


m=  1 


Note  that  for  any  1  <  m  <  k,  \hsdjm\2  is  an  exponential  random  variable  with  parameter  a sd  ,  so 


lim  si  •  Pr  [um  <  t]  =  lim  si  •  Pr 

Si — >oo  Si  — >oo 


I  hsd 


I2  < 

,m\  - 


2*  -  1 
si 


=  ^-(2t-l),  m  =  1,2 


,k. 


a 


(31) 


sd 


Since  urn ,  1  <  m  <  k,  arc  independent  random  variables,  by  applying  Lemma  1  with  M  =  1  recursively,  we 
have 

k  1  , 

=  (PY 


lim  Si  •  Pr 

Si— XXD 


J2um<t 


\_m=l 


x  I 


9k{t), 


(32) 


where  gk(t)  is  given  in  (19).  For  any  m,  k  +  1  <  m  <  L,  um  involves  the  sum  of  two  independent  exponential 
random  variables  hsdim | 2  and  hrd,Tn 1 2  with  parameters  a~d  and  a~d,  respectively,  and  the  distribution  of  um 
can  be  specified  as 

Pr  [um  <t]  =  Pr  [si\hsdtm\2  +  s2\hrdtm\2  <  2t  —  l] 

1  2*  — 1 ' 


exp 


1_  + 

V  ^ 

1  ~  ^sd-^rd0  ^  _ 


/  1  2‘-lN| 

va  pyn  ( 

1  2*-l\ 

v  / 

S2^d-S1(T^  P  l 

^  J 

if 

^rd  ^sd 

if  4L  X  4jL, 

CTrd  CTid 


for  any  m  =  A:  +  1, L.  Thus,  for  any  m  =  k  +  1, L,  we  have 

slinroSiS2-Pr[um<f]  =  — ^ 

l<i<2  asdard 


(24  -  l)2 


(33) 


(34) 


Let  go(A)  =  1  and  p(t)  =  (2t  —  l)2,  then  p'(t )  =  2(22*  —  24)ln2.  Since  um,  k  +  1  <  m  <  L,  are  independent 
to  each  other,  by  applying  Lemma  1  with  M  =  2  recursively,  we  can  show  that  for  any  n  =  1,  2, L  —  k, 

k+n 

y.  um<t  -~ 


lim  (si-S2)ra  •  Pr 

Si — >oo  x  7 

l<i<  2 


\_m=k+l 


Qn{t)  j 


(35) 


in  which 


qn(t)  =  /  qn-i(x)p(t  -  x)dx,  n  =  1,2, L  —  k. 

Jo 

Based  on  Lemma  2,  a  closed-form  expression  for  qn(t)  can  be  obtained  as  follows: 

rt  rxn-i  rx 2 

qn(t)  =  //  •••/  Ql(.Tl)p,(x2  -  Xl)p'(x3  ~  X2)  ■  ■  -p'ixn-l  ~  Xn-2) 

J0J0  Jo 

xp(t  —  xn_i)dxidx’2  •  •  •  dxn_i 

■2  n~ 2 
(2X1  -  l)2(21n2)n_1  f2Xm+1~Xrn  -  l'' 


(36) 


rt  r. x„~i  rx 2 

I0J0  Jo 


m= 1 


x^2i-a;n_i  _  f)2t_a:idx’1d.'C2  •  •  -dxn-i 

(— 21n2)n_1  E  - |-On-12(l+a™-l)i 


6(0,1} 


rt  rxn- 1  fX2 

I0J0 


n—  1 


/•x2  "■  x 

■  (2X1  -  l)2  •  2-"1*1  TT  2(am-1_am)a;mdx’idx2  •  •  -dsn-i 

m=2 

(— 21n2)n_1  ^  (-l)Ql+-+“"-12(1+a"-^{Fn_i(l  -ai,02,...,l3n-i\t) 

ai,...,an_i 

G{0,1} 

2Fn—  1  (  CTl ,  (32 1  ■■■1  fin—  1  j  A)  T  Pn—  1  (  1  Ob >  /?2 >  ■  fin—  1  i  A)  }  j 


(37) 
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where  /%  =  «i  —  o>2,  P3  =  «2  —  «3,  /3n-i  =  «n-2  —  an-i.  and  the  function  Fn_ i(  •  ;  i)  is  defined  in  (9) 

and  it  has  a  closed-form  expression  specified  in  Lemma  2.  Note  that  in  case  for  some  i,  Pi  is  zero  which  does 
not  satisfied  the  non-zero  condition  in  Lemma  2,  we  may  calculate  Fn- 1(  •  ;  t)  as  follows.  Since  the  function 
Fn((3 1, ....  Pn:  t)  defined  in  (9)  is  continuous  in  terms  of  each  variable  Pt  and  the  closed-form  expression  in 
Lemma  2  is  also  continuous  in  terms  of  Pi,  so  we  can  apply  Lemma  2  with  Pi  =  £,  where  e,  is  sufficiently 
small  (i.e.  sy  — ►  0)  to  accommodate  the  calculation  of  Fn_\(  ■  ;  t). 

According  to  the  result  in  (35)  with  n  =  L  —  k,  we  have 


lim  (•S1S2)1'  k  ■  Pr 

Si— >00  v 
l<i<2 


‘  L 

y  um<t 

_m=k-\- 1 


2asdard 


L-k 

QL-k(t), 


(38) 


where  qL-k(t )  can  be  calculated  specifically  based  on  (37).  Combining  (32)  and  (38),  and  applying  Lemma  1, 
we  obtain 


lim  sf  Sn  k  •  Pr 

Si~ >00  ±  z 

l<i<2 


^  "  Isd,m  4~  ^  "  Isrd,m  F 


R 


m=  1 


m=k+ 1 


(39) 


where 


&fc(*)  =  /  9k{x)qL_k(t  -  x)dx. 


(40) 


Since  .sy  =  Ps/Nq  and  S2  =  Pr/-N 0,  so  for  any  Tr  =  k,  k  =  1,2, ...,  L  —  1,  the  conditional  probability  (27) 
can  be  asymptotically  approximated  as 


Pt%  =  Pr 


^  Isd,v 


L  771=1 


+ 


L 

E 

m=/c+l 


Isrd.m  ^ 


6fc(i2)  f  -A/o 

2L-k 


<dPs 


•A/p 


L-fc 


(41) 


Finally,  combining  (22)-(24)  and  (41),  we  can  obtain  the  outage  probability  for  the  DF  cooperative  ARQ 
relay  scheme  as  follows 


P‘ 


out,L 


L—l 

E 

fc=i 


bk(R ) 

2L-k 


gk-i(R)  -  gk(R) 


+  9L{R)gL-i(R)  ( 


<dps 


A/q 

O’sr-Ps 

A/~o 


A/q 

L-l 


A/~o 


L-fc 


A/q 

0‘srA’s 


fc-1 


(42) 


where  gk(R)  and  bk(R)  are  specified  in  (19)  and  (40),  respectively.  Furthermore,  we  note  that  the  term 
9k(R)  y°P  would  be  much  smaller  than  gk-i(R)  at  high  SNR  so  the  asymptotic  outage  probability  in 
(42)  can  be  further  simplified  as 


,out,L  yA  bk(R)gk-i(R)  (  A/o  \L  (  A/o  \L“V  A/o  \fc_1 


(43) 


where  6l(7?)  =  giPJl).  The  asymptotic  outage  probability  is  tight  at  high  SNR  which  will  be  shown  in 
simulation  results. 

Based  on  the  above  asymptotic  outage  probability,  we  observe  that  the  term  (  )  in  (43)  contributes  a 

V  ® sd  s  / 

diversity  order  L  in  the  asymptotic  outage  performance,  which  is  due  to  the  fast  fading  nature  of  the  channels. 


The  term 


Mo 


L-k 


Mo 

tFPs 


fc-1 


contributes  an  overall  diversity  order  (L  —  k)  +  (k  —  1)  =  L  —  l  which  is  due 


to  the  cooperative  relaying.  Thus,  the  asymptotic  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme 
has  an  overall  diversity  order  2L  —  1.  The  observation  is  more  obvious  in  case  of  an  equal  power  allocation 
scenario,  i.e.,  Ps  =  Pr  =  P.  As  a  comparison,  we  recall  that  the  diversity  order  of  the  direct  ARQ  transmission 
scheme  is  only  L,  which  is  much  less  than  that  of  the  DF  cooperative  ARQ  relay  scheme.  The  diversity  order 
of  the  direct  ARQ  scheme  comes  from  the  fast  fading  nature  of  the  source-destination  channel. 


16 


3.1.3  Optimum  Power  Allocation  for  the  DF  Cooperative  ARQ  Relay  Scheme 

In  this  subsection,  we  derive  an  asymptotic  optimum  power  allocation  for  the  DF  cooperative  ARQ  relay  scheme 
based  on  the  tight  approximation  of  the  outage  probability  we  presented  in  the  previous  section.  Without  loss 

of  generality,  we  denote  the  total  transmission  power  Ps  +  Pr  =  2 P.  For  any  given  total  transmission  power 
2 P,  we  try  to  determine  optimum  power  Ps  used  at  the  source  and  power  Pr  at  the  relay  in  order  to  minimize 
the  asymptotic  outage  probability. 

Let  us  denote  A  as  the  ratio  of  the  transmit  power  Ps  to  the  total  transmission  power,  i.e.  A  =  jp,  then 
0  <  A  <  1  and  Pr  =  (1  —  A)2 P.  The  asymptotic  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme 
can  be  written  as 


P' 


out,L 


at 


(2(7 


2  rr2  )L 

sd°rd> 


A/q 
2 P 


2L-1  L 

^hWgk-^R) 

k= l 


Ha 


cr2r  J  AL+fc_1(l  -  X)L~k‘ 


(44) 


We  try  to  find  the  optimum  power  ratio  A  (0  <  A  <  1)  such  that  the  asymptotic  outage  probability  is  minimized. 

A  /  9 rr2  \  k 

Let  Ak(R)  =  b^g^R)  and 


p/n  A  v-  Ak{R) _ 

(  (  AL+fe-i(l  _  X)L~k' 

then  the  optimization  problem  can  be  formulated  as  follows 

min  G(A) 

A 


s.t.  0  <  A  <  1. 


By  taking  derivative  of  G( A)  with  respect  to  A,  we  have 


dG{  A) 
dX 


J2mr) 


k= 1 


+ 


L-k 


-(L  +  fc-1) 

XL+k(l  —  X)L~k  '  AL+fc_1(l  —  A)L_fc+1  / 


Let  c)Cg^  =  0,  an  optimum  power  ratio  can  be  found  by  solving  the  following  equation 


Y^Ak(R)\-(L  +  k-  1) 


k=  1 


1  -  A 


+  (L  —  k) 


1  -  A 


k- 1' 


=  0. 


(45) 


(46) 


(47) 


(48) 


The  equation  can  be  easily  solved  by  using  the  Newton  method.  Based  on  the  equation,  we  have  two  observa¬ 
tions.  First,  since  Ak(R)  is  positive,  so  1  in  the  equation  must  be  less  than  1,  otherwise  the  left-hand  side 
of  (48)  is  negative.  It  implies  that  A  >  i.e.  Ps  >  P  and  Pr  <  P,  which  means  we  should  allocate  more 
power  at  the  source  and  less  power  at  the  relay.  It  also  shows  that  the  equal  power  allocation  scheme  assigning 
power  equally  between  the  source  and  the  relay  is  not  optimum  in  general.  Second,  we  observe  that  for  a  given 
transmission  rate  R,  the  parameters  Ak(R)  in  the  equation  (48)  depend  only  on  <j2sr  and  a2d  which  are  the 
variances  of  the  source-relay  and  relay-destination  channel  links,  respectively.  Thus,  the  asymptotic  optimum 
power  ratio  A  depends  only  on  the  the  variances  of  the  source-relay  and  relay-destination  channels,  not  on  the 
source-destination  channel  link.  A  similar  observation  was  reported  in  [7]  where  the  optimum  power  allocation 
between  the  source  and  the  relay  was  determined  based  on  the  analysis  of  the  symbol-error-rate  performance. 

When  L  =  2,  a  closed-from  expression  of  the  optimum  power  ratio  can  be  obtained  as  follows.  In  this  case, 
the  equation  (48)  is  reduced  as 

3 A2(R)  +  2Ai (R)  "  MR)  =  0,  (49) 
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so  the  optimum  power  ratio  is 


A 


1  + 
2  + 


V^3 


A2(R) 

Ai(R) 


\/1  +  3 


A2(R)  ' 
Ai(R) 


Thus  the  coiTesponding  optimum  power  allocation  at  the  source  and  at  the  relay  is  given  by 


(50) 


Ps 


Pr 


1  +  1/1  +  3 


A2{R) 

Ai(R) 


P, 


l  +  \A++{f 

1 


1  +  5V/l  +  3SW 


P. 


(51) 

(52) 


In  (51)  and  (52),  the  parameter  =  7^-  So  the  optimum  power  allocation  dependents  only 

on  the  ratio  of  the  variances  of  the  source-relay  and  relay-destination  channels.  This  result  is  consistent  with 
the  observation  based  on  the  general  optimization  equation  in  (48).  Also,  from  (51)  and  (52),  we  can  see  that 
P  <  Ps  <  2 P  and  0  <  Pr  <  P,  i.e.  we  should  put  more  power  at  the  source  and  less  power  at  the  relay 
to  optimize  the  overall  performance  at  the  destination.  Furthermore,  from  (51)  and  (52)  we  observe  that  if  the 
relay  is  located  close  to  the  source,  i.e.,  a2r  +>  a2d,  Ps  goes  to  |P  and  Pr  goes  to  |P.  On  the  contrary,  if 
the  relay  is  located  close  to  the  destination,  i.e.,  a2r  <C  a2d,  Ps  goes  to  2 P  and  Pr  goes  to  0,  which  means  we 
should  allocate  most  power  at  the  source  in  this  case.  It  is  reasonable  to  allocate  most  of  the  total  power  2 P  to 
the  source  since  the  role  of  the  relay  helping  in  forwarding  is  minor  in  this  case. 


3.1.4  Simulation  Results  and  Discussion 

In  this  section,  we  present  numerical  and  simulation  studies  for  the  DF  cooperative  ARQ  relay  scheme  to 
validate  our  theoretical  analysis  in  the  previous  sections.  We  also  compare  the  performance  of  the  DF  co¬ 
operative  ARQ  scheme  with  that  of  the  direct  ARQ  scheme.  In  all  studies,  the  variance  of  the  channel 
hij  {(*,  j )  G  (s,  d),  ( s ,  r),  (r,  d)}  is  assumed  to  be  cr|-  =  dP\  where  dij  is  the  distance  between  two  nodes 
and  //  is  the  path  loss  exponent  which  is  assumed  to  he  //  =  3  in  a  typical  fading  environment.  We  assume  that 
the  source-destination  distance  is  dsci  =  10  m  and  the  relay  is  located  in  the  midpoint  between  the  source  and 
the  destination.  We  consider  a  target  transmission  rate  of  R  =  2  bits/s/Hz. 

Figs.  2,  3  and  4  present  the  simulation  studies  for  the  DF  cooperative  ARQ  relay  scheme  when  the  maximum 
number  of  ARQ  retransmission  rounds  is  L  =  2,  3  and  4,  respectively.  In  these  simulations,  we  allocated 
power  equally  at  the  source  and  at  the  relay,  i.e.,  Ps  =  Pr  =  P.  The  simulation  results  show  that  the  theoretical 
approximation  of  the  outage  probability  of  the  DF  cooperative  ARQ  relay  scheme  is  loose  at  low  SNR  and  tight 
at  high  SNR.  For  example,  in  case  of  L  =  2  in  Fig.  2,  the  analytical  approximation  matches  with  the  simulated 
curve  at  an  outage  performance  of  around  10  3.  Moreover,  the  larger  the  number  of  ARQ  retransmission 
rounds,  the  higher  the  diversity  order  of  the  DF  cooperative  ARQ  relay  scheme.  This  observation  is  consistent 
with  the  theoretical  result  that  the  diversity  order  of  the  DF  cooperative  ARQ  relay  scheme  increases  in  terms 
of  the  number  of  ARQ  retransmission  rounds. 

For  comparison.  Figs.  2-4  also  include  the  performance  of  the  direct  ARQ  scheme.  We  can  see  that  the  DF 
cooperative  ARQ  relay  scheme  significantly  outperforms  the  direct  ARQ  scheme.  At  an  outage  performance  of 
ICR4,  the  performance  of  the  DF  cooperative  ARQ  relay  scheme  is  about  8dB  better  than  that  of  the  direct  ARQ 
scheme.  For  the  same  maximum  number  of  retransmission  rounds  P,  the  performance  of  the  DF  cooperative 
ARQ  relay  scheme  shows  a  higher  diversity  order  than  that  of  the  direct  ARQ  scheme.  This  observation  is 
consistent  with  our  theoretical  developments  showing  that  the  DF  cooperative  ARQ  relay  scheme  has  diversity 
order  2L  —  1,  while  the  direct  ARQ  scheme  has  diversity  order  only  L. 

In  Fig.  5,  we  show  the  optimum  power  ratio  A  for  the  DF  cooperative  ARQ  relay  scheme  with  L  =  2,  3 
and  4,  respectively.  We  plot  the  optimization  function  G( A)  in  terms  of  A  (0  <  A  <  1).  We  assume  that  the 
quality  of  the  source-relay  link  is  the  same  as  that  of  the  relay-destination  link,  i.e.,  a2r  =  a2d.  We  observe 
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P/JV0(dB) 

Figure  2:  Outage  probability  of  the  direct  and  DF  cooperative  ARQ  schemes  (L=2). 


P/JV0(dB) 

Figure  3:  Outage  probability  of  the  direct  and  DF  cooperative  ARQ  schemes  (L=3). 
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Figure  4:  Outage  probability  of  the  direct  and  DF  cooperative  ARQ  schemes  (L=4). 


Figure  5:  Optimum  power  ratio  A  for  the  DF  cooperative  ARQ  scheme.  When  L=2,  3  and  4,  the 
asymptotic  power  allocation  is  A  =  0.8203,  A  =  0.7969  and  A  =  0.7838,  respectively. 
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P/N0(dB) 

Figure  6:  Outage  probability  of  the  DF  cooperative  ARQ  scheme  with  equal  and  optimum  power 
allocations  (L=2). 


P/N0(dB) 

Figure  7 :  Outage  probability  of  the  DF  cooperative  ARQ  scheme  with  equal  and  optimum  power 
allocations  (L=3). 
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Figure  8:  Outage  probability  of  the  DF  cooperative  ARQ  scheme  with  equal  and  optimum  power 
allocations  (L=4). 


that  the  optimum  power  ratio  A  is  around  0.8  for  the  three  cases.  More  precisely,  from  the  numerical  results, 
the  optimum  power  ratios  arc  A  =  0.8203,  A  =  0.7969  and  A  =  0.7838  for  L  =  2,  3  and  4,  respectively. 
It  appears  that  the  optimum  power  ratio  decreases  gradually  when  the  maximum  number  of  retransmission 
rounds  is  increased.  Furthermore,  the  optimum  power  ratio  is  much  larger  than  1/2  which  is  consistent  with 
our  analysis  that  we  should  allocate  more  power  at  the  source  and  less  at  the  relay. 

We  compare  the  performances  of  the  the  DF  cooperative  ARQ  relay  scheme  with  the  optimum  power  alloca¬ 
tion  as  well  as  the  equal  power  allocation  in  Figs.  6-8  L  =  2,3  and  4,  respectively.  In  case  of  L  =  2  in  Fig. 
6,  the  optimum  power  allocation  is  PS/2P  =  0.8203  and  Pr/2P  =  0.1797  according  to  the  numerical  results 
in  Fig.  5.  From  the  figure,  we  observe  that  the  performance  of  the  DF  cooperative  ARQ  relay  scheme  with  the 
optimum  power  allocation  is  about  1.5dB  better  than  that  of  the  scheme  with  the  equal  power  allocation.  In 
cases  of  L  =  3  in  Fig.  7,  the  optimum  power  allocation  is  PS/2P  =  0.7969  and  Pr/2P  =  0.2031,  and  the 
optimum  power  allocation  shows  a  performance  improvement  of  1 .25dB  compare  to  the  equal  power  allocation. 
In  case  of  L  =  4  in  Fig.  8,  the  optimum  power  allocations  is  PS/2P  =  0.7838  and  Pr/2P  =  0.2162.  We  can 
see  that  the  DF  cooperative  ARQ  relay  scheme  with  the  optimum  power  allocation  also  shows  a  1.25dB  gain 
compared  to  that  with  the  equal  power  allocation.  In  Figs.  6-8,  we  also  plot  the  approximation  of  the  outage 
probability  for  the  DF  cooperative  ARQ  relay  scheme  with  the  optimum  power  allocation  as  well  as  the  equal 
power  allocation.  We  can  see  that  the  approximation  of  the  outage  probability  matches  with  the  simulation 
curves  at  high  SNR  in  all  scenarios,  which  further  validate  our  theoretical  analysis. 

3.2  Optimal  Power  Assignment  for  Hybrid-ARQ  Rayleigh  Fading  Links 

In  this  subsection,  first  we  review  briefly  the  H-ARQ  transmission  scheme  and  formulate  the  power  assignment 
optimization  problem.  Second,  we  find  the  optimal  power  assignment  strategy  for  the  H-ARQ  protocol  and 
present  an  exact  recursive  calculation  algorithm.  Third,  we  develop  a  simple  approximation  of  the  optimal 
power  assignment  sequence  and  compare  it  with  the  exact  calculation  result.  By  the  end  of  this  subsection, 
numerical  and  simulation  studies  are  carried  out  to  compare  the  performance  of  the  equal  and  optimal  power 
assignment  strategies. 
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Figure  9:  Illustration  of  a  hybrid- ARQ  protocol  with  transmission  power  Pi  in  the  Zth  (re- 
)transmission  round,  1  <  l  <  L. 


3.2.1  System  Model  and  Problem  Formulation 


We  consider  an  H-ARQ  transmission  protocol  implemented  between  a  source  node  and  a  destination  node  as 
illustrated  in  Fig.  9.  Assume  that  L  is  the  number  of  retransmission  rounds  allowed  in  the  H-ARQ  protocol. 
The  H-ARQ  transmission  scheme  operates  as  follows.  First,  the  source  transmits  an  information  packet  to  the 
destination  and  the  destination  indicates  success  or  failure  of  receiving  the  packet  by  feeding  back  a  single 
bit  of  acknowledge  (ACK)  or  negative-acknowledgement  (NACK),  respectively.  The  feedback  channel  is  as¬ 
sumed  error-free.  Then,  if  a  NACK  is  received  by  the  source  and  the  maximum  number  of  retransmissions 
L  is  not  reached,  the  source  retransmits  the  packet  at  a  potentially  different  transmission  power  to  be  deter¬ 
mined/optimized.  If  an  ACK  is  received  by  the  source  or  the  maximum  retransmission  number  L  is  reached, 
the  source  begins  transmission  of  a  new  information  packet.  In  each  retransmission  round,  the  destination  at¬ 
tempts  to  decode  an  information  packet  by  combining  received  signals  from  all  previous  transmission  rounds 
by  the  standard  maximal -ratio-combining  (MRC)  technique  [37].  If  the  destination  still  cannot  decode  an  infor¬ 
mation  packet  after  L  (re)transmission  rounds,  then  an  outage  is  declared  which  means  that  the  signal-to-noise 
ratio  (SNR)  of  the  combined  received  signals  at  the  destination  is  below  a  required  SNR. 

The  H-ARQ  transmission  scheme  can  be  modeled  as  follows.  With  L  maximum  retransmission  rounds 
allowed  in  the  H-ARQ  protocol,  the  based-band  received  signal  ysd,i  at  the  destination  at  the  Zth  transmission 
round  can  be  written  as 

Vsd,i  —  C  Pi  hsd  xs  T  T]sd,h  Z  —  1;  2,  •  ,  L ,  (53) 


where  xs  is  the  transmitted  information  symbol  from  the  source.  Pi  is  the  transmission  power  used  by  the 
source  at  the  Zth  transmission  round,  hsd  is  the  source-destination  channel  coefficient,  and  rjsdti  is  additive  noise 
at  the  Zth  round.  The  channel  coefficient  hsd  is  modeled  as  zero-mean  complex  Gaussian  random  variable  with 
variance  a2sd.  The  channel  is  assumed  to  be  quasi-static ,  i.e.  the  channel  does  not  change  during  retransmissions 
of  the  same  information  packet  and  it  may  change  independently  when  a  new  information  packet  is  transmitted. 
The  source-destination  channel  coefficient  is  assumed  to  be  known  at  the  receiver  side,  but  unknown  at  the 
transmitter  side.  The  additive  noise  contribution  fjsd.l  is  modeled  as  a  zero-mean  complex  Gaussian  random 
variable  with  variance  Mq. 

At  the  destination  side,  the  receiving  node  combines  the  received  signals  from  all  previous  retransmission 
rounds  and  jointly  decodes  the  information  packet  based  on  the  MRC  combining  technique  [37].  Note  that  the 
MRC  combining  is  applied  over  base-band  symbol-level  signals  in  (53)  before  decoding  an  entire  information 
packet.  With  the  assumption  that  the  channel  does  not  change  in  retransmissions  of  the  same  information 
packet,  the  SNR  of  the  combined  signal  at  the  destination  at  the  Zth  (1  <  Z  <  L)  retransmission  round  can  be 
given  as  [37],  [38] 


' 7sd,l 


Etifll  hsd\2\xs\2 
A/q 


(54) 


Without  loss  of  generality,  let  us  assume  the  average  power  of  the  transmitted  information  symbol  is  1,  then 
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we  have  7,,//  =  ^~‘i=1^hsd^  •  Since  /;,sr/  follows  a  Rayleigh  distribution  with  mean  zero  and  variance  a2sd, 
so  for  any  targeted  SNR  70,  the  probability  of  the  event  that  the  destination  cannot  decode  correctly  after  l 
transmission  rounds  can  be  calculated  as 

t'qA/'o 

V°uU  =  Pr  [7 sd,i  <  7o]  =  1  -  e  A .  (55) 

Set  pout’°  =  1.  Then,  the  probability  that  the  H-ARQ  protocol  stops  successfully  at  the  /th,  1  <  l  <  L, 
transmission  round  is  pouiJ  1  —  pnuL1\  which  means  the  destination  cannot  decode  correctly  at  the  (1  —  l)th 
round,  but  succeeds  at  the  It h  round. 

Our  goal  is  to  find  an  optimal  power  assignment  sequence  P  =  [Pi,  P2,  Pj\  for  the  H-ARQ  protocol  such 
that  under  a  targeted  outage  probability  po,  the  average  total  transmission  power  for  the  protocol  to  deliver  an 
information  packet  is  minimized.  Since  the  probability  that  the  protocol  succeeds  exactly  at  the  /th  (1  <  l  < 
L  —  1)  round  is  pout’1- 1  —  pnvLl  and  the  coiTesponding  total  transmission  power  is  P\  +  P2  +  •  •  •  +  P/,  so  the 
average  total  transmission  power  of  the  H-ARQ  protocol  can  be  expressed  as 


P  = 


£  {pout'l~l 
1=1 


p 


out, l 


1=1 


i= 1 


(56) 


Note  that  the  last  term  in  (56)  is  due  to  the  fact  that  the  protocol  stops  retransmissions  after  the  Pth  round  no 
matter  whether  decoding  at  the  Pth  round  is  successful  or  not.  For  the  H-ARQ  protocol  with  a  targeted  outage 
probability  po,  the  problem  of  finding  optimal  power  assignment  can  be  formulated  as  follows: 


min  P  with  respect  to  Pi ,  P2,  ■  ■  •  ,  Pl  >  0 
subject  to  pout’L  <  p0 


(57) 


where  P  is  specified  in  (56). 

3.2.2  Optimal  Transmission  Power  Assignment 

In  this  section,  we  investigate  the  optimal  power  assignment  strategy  for  the  H-ARQ  protocol  to  minimize  the 
average  total  transmission  power.  We  obtain  a  set  of  equations  that  describe  the  optimal  transmission  power 
values,  and  then  develop  a  recursive  algorithm  to  exactly  calculate  the  optimal  transmission  power  level  for 
each  retransmission  round. 

The  average  total  transmission  power  in  (56)  can  be  rewritten  by  switching  the  summation  order  (between 
the  indices  l  and  i)  as  follows 


p  =  ±p, 

i=  1 

where  we  first  consider  the  summation  by  enumerating  the  index  i  from  1  to  P,  then  consider  the  summation 
index  l  (i  <  l  <  L  —  1).  Since  for  each  i,  —  pout )l)  =  pout^~l  —  pOut,L-i ,  so  average  total 

transmission  power  can  be  represented  as 


L- 1 

E 

L  l=i 


~out.l — 1 


p 


-pout^  +p 1 


out.L—  1 


(58) 


L 

P  =  Pi  +  £  Pi  pout ’l~1.  (59) 

1=2 

Moreover,  the  constraint  in  (57)  means  that  with  a  targeted  SNR  70,  the  outage  probability  of  the  H-ARQ 
protocol  with  P  retransmissions  should  not  be  larger  than  the  specified  outage  probability  value  77,  i.e. 

_  70-Vh 

pout,L  =  1_e  PdEf=i7  <  pQ'  (60) 
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Denote  Pq  =  — s  then  the  constraint  is  equivalent  to 

u  sa  11 1—pQ 

L 

1=1 

and  the  optimization  problem  in  (57)  can  be  further  specified  as 


L 

minp1;... ,pl> o  P  =  P1+y^Pi 

1=2 


L 

subject  to  ^  Pi  >  Po- 

i=i 


(61) 


(62) 


Next,  we  relax  temporarily  the  non- negative  condition  on  Pi,  l  =  1, 2,  •  •  •  ,  L,  and  consider  the  sum-power 
constraint  in  (62)  with  equality  in  order  to  consider  a  Lagrange  multiplier  method  to  solve  the  optimization 
problem.  We  will  prove  later  that  the  obtained  solution  is  indeed  optimal  under  the  constraint  in  (62)  and  it 
satisfies  the  non-negative  condition.  Let  us  form  a  Lagrangian  objective  function  as 


L 

70-V*0 

r  L  -| 

£(P,A)  =  ftL  +  ^ft 
1=2 

1  (T2  1  P- 

1  —  g  sd  ^  i=l  4 

+  A 

T.o-p„ 

L  t=t  J 

(63) 


Taking  the  derivative  of  C  (P,  A)  with  respect  to  A  and  setting  it  equal  to  zero,  we  have  the  power  constraint  as 
1  Pi  —  Po  =  0.  The  derivatives  of  C  (P,  A)  with  respect  to  ft,  are 


dC 

dP i 

dC 

dPk 

dC 

~dP~L 
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ToA^n 


1=2  asd(^2i= 1  ^*) 


e  a?dEi=l  Pi  +  A, 
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1  /r2  ^  p. 
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PiloMo 


l=k+ 1  (  Ei=  1  ft) 


_  7pA/h 

e  °sdft=i  pi  -)-  A, 


k  =  2,  3, 


ft 


_ -I'pA/p 

2  ^  —  Id 

1  —  e  °sd^-i=l  2  -f-  A. 


Based  on  =  0  and  =  0,  we  have 


dC  dC 


dP\  dP2 


1  - 


E270A/0 


sd 


70  Wp 


e  CT*dpi  =  0, 


which  implies 


P2  = 
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sd 
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For  any  k  =  3, 4,  •  •  •  ,  L,  according  to  ap  -  —  ,,  ap 


7&M)  ’ 

C>c  =  0  and  Mfr  =  0,  we  have 
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ftft-i  dPk 
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which  means 


7o2Vo 


ft  = 


^vESft)  J 

pfc-l7oAo 


2  V'  /c  —  1  n 

e  %d^i=i  pi  =  0, 


1  _  e  °2sd(T.ki=l  pi){p.ki=t  Pi) 


(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 
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for  any  k  =  3, 4,  •  •  •  ,L.  We  can  easily  verify  that  the  Lagrangian  solutions  p>.  P>,  ■  ■  ■  ,Pl  in  (68)  and  (70)  arc 
positive. 

In  the  following,  we  would  like  to  show  that  the  average  total  transmission  power  P  cannot  be  further  min¬ 
imized  with  strict  inequality  in  (62).  If  there  exists  a  power  sequence  Pf,  P|,  •  •  •  ,  P£  such  that  P*  +  P|  + 
■■■  +  n  >  P)  and  the  average  total  transmission  power  P  is  minimized,  then  let  us  consider  another  power 
sequence 


pk  =  rpk;  k  =  1,2,- ••  ,L, 


in  which  r  is  an  arbitrary  number  satisfying 

r  > 


Pn 


P*  +  P*  - b  P£ 


(71) 


(72) 


We  can  see  that  the  new  power  sequence  P\,  P2,  ■  ■  ■  ,  Pl  satisfies  the  power  constraint  in  (62).  With  the  new 
power  sequence,  the  resulting  average  total  transmission  power  is 


A 


70-Vq 


f(r)  = 

1= 2 

L 

rp;  +  YJrPi 


2  yl-1  p, 
1  —  g  sd  ^i=  1  * 


1=2 


7qATq 


1  a2  Y1  —  1  rP* 

\  —  g  asd  ^  i=l  r  4 


(73) 


which  is  a  function  of  r.  Taking  derivative  of  f(r)  with  respect  to  r,  we  have 

L  L 


9f(r) 

dr 


T,pi~T,pi 


1=1 


1=2 


1  + 


7oM) 


rPi 


-rn-wn 

G2  1  rP* 

g  sd^i=l  t  . 


(74) 


Since  e  X(1  +  x)  <  1  for  any  positive  x,  so  in  (74), 
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EEi  rPt 
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sd 


T ()-'Vt) 
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sdL-‘i=  1  1  1 . 


(75) 


Thus,  we  have 


d  f(r) 
dr 


L 


>EP‘ 


L 


pi  >  0, 


(76) 


which  means  that  f(r)  is  an  increasing  function  for  any  r  >  p*+P*+...+P* ,  and  the  minimum  of  f(r)  is 
achieved  when  r  =  P*+P^...+P*  ■  It  implies  that  the  average  total  transmission  power  resulting  from  the  new 
power  sequence  with  r  =  p*+P*+...+P*  (<  1)  is  less  than  that  based  on  the  power  sequence  P*,  •  •  •  ,  Pj. 

This  is  contradictory  to  the  assumption  that  the  power  sequence  P± ,  P% ,  ■  ■■  ,  P£  minimizes  the  average  total 
transmission  power  P.  Therefore,  the  minimum  average  total  transmission  power  P  can  be  achieved  at  the 
boundary  of  the  constraint  (with  equality)  in  (62).  We  note  that  with  r  =  P*+P*p...+P* ,  the  new  power 
sequence  satisfies 

-Pl  +  p2  +  •  •  •  +  pL  =  Pqj  (77) 


which  is  the  boundary  of  the  constraint  in  (62). 

We  note  that  in  general  a  Lagrangian  solution  may  not  guarantee  global  optimality,  i.e.  it  may  lead  to  a  local 
minima  or  maxima.  Fortunately,  the  Lagrangian  solution  in  (68)  and  (70)  leads  to  a  global  minima  as  explained 
as  follows.  From  the  Lagrangian  solution  in  (68)  and  (70)  and  the  total  power  constraint  Pi+P2+-  •  -+Pl  =  Pi, 
we  can  see  that  there  is  only  one  unique  power  sequence  Pi,  P2,  ■ ' '  ■  pr.  that  results  from  the  Lagrangian 
solution.  So,  the  unique  power  sequence  guarantees  the  global  optimality  which,  however,  can  be  either  global 
minima  or  maxima.  We  further  exam  that  with  a  trivial  power  assignment  Pi  =  Pq,  P2  =  P3  =  •  •  •  =  Pl, 


26 


Table  1:  Algorithm  to  determine  the  optimal  power  assignment  sequence  I\  ( 1  <  k  <  L ) 


Step  1 
Step  2 


Input  7o,po:  v'sri-No-  Calculate  Pq  = 
Set  lower  =  0  and  upper  =  Pq- 


7qA/q 

j2,ln  ^ — ] — 

sd  1—  PQ 


Step  3  :  Let  P\  =  (lower  +  upper)/ 2,  and  calculate 
70A/0 


_2  p2 

p2  =  Adpi 


70-W) 

for  any  k  =  3, 4,  •  •  •  ,  L. 


pfc- l'TO-A) 


l_e  pi) 


4  :  Check  if  abs(P\  +  P2  +  •  •  •  +  Pl  —  Po)  <  e(=  0.001),  then  stop 
and  output  power  sequence  Pf ,  P2 ,  ■  •  •  ,Pl ;  otherwise, 
if  P\  +  P2  +  ■  ■  ■  +  Pl  —  Po  <  0,  set  loiver  =  Pi; 

z/  Pi  +  P2  H - h  Pl  -  Po  >  0,  set  upper  =  Pi; 

and  go  to  Step  3. 


the  resulting  average  total  transmission  power  is  P  =  Po,  which  is  larger  than  the  average  total  transmission 
power  resulting  from  the  power  sequence  associated  with  the  Lagrangian  solution.  Therefore,  the  unique  power 
sequence  from  the  Lagrangian  solution  guarantees  the  global  minima.  We  summarize  the  above  discussion  in 
the  following  theorem. 


Theorem  1  In  the  H-ARQ  transmission  protocol,  to  minimize  the  average  total  transmission  power,  the  optimal 
transmission  power  If  at  the  kth,  1  <  k  <  L,  transmission  round  must  satisfy  the  following 


P2  = 


rr2  P2 
asdr 1 

7oA/q  ’ 


y-UY^Zl  Pi 

7oM) 


Pk  = 


pk-  mWo 


1  —  e  <Jsd 


(E  til  Pi)(^i=t  Pr) 


for  k  =  3, 4,  •  •  •  ,  L,  and 


Pi  +  P2  H - 1-  Pl  =  P( 


0, 


(78) 


(79) 


(80) 


where  Po  =  2'^~0i — >  7o  P  die  required  SNR  for  correct  decoding,  Nq  is  the  additive  white  noise  variance, 

u  sdin  1—  pq 

cr2sd  is  the  Rayleigh  fading  variance,  and  po  is  the  targeted  H-ARQ  outage  probability. 


□ 


From  Theorem  1,  we  can  see  that  the  optimal  transmission  power  sequence  is  uniquely  determined  by  the  set 
of  equations  (78)— (80).  The  optimal  transmission  power  level  for  each  (re)transmission  round  can  be  calculated 
recursively.  According  to  (78)  and  (79),  for  any  k  =  2, 3,  •  •  •  ,  L,  the  optimal  transmission  power  value  Pl  can 
be  calculated  based  on  P\ .  P2.  ■  ■  •  ,  P/._  | .  So  for  any  given  power  Pi,  all  other  transmission  power  P/c ,  k  = 
2, 3,  •  •  •  ,  L,  can  be  subsequently  determined.  The  optimal  initial  power  Pi  can  be  numerically  found  based 
on  (80)  by  the  Newton  method.  A  complete  algorithm  to  recursively  determine  the  optimal  power  assignment 
sequence  Pl,  k  =  1, 2,  •  •  •  ,  L,  is  detailed  in  Table  I. 

When  L  =  2,  we  have  a  closed-form  solution  for  the  optimal  power  sequence.  In  this  case,  P\  +  If  =  p. 
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and  P‘2 


rr2  P2 

crsd±'l 

7oWo  ' 


By  solving  the  two  equations,  the  optimal  transmission  power  P\  and  P>  arc  given  by 


Pi 


P-2 


2P0 

1  ,  .  /l  ,  4asdP0  ’ 

1  +  V1  +  TTW 

70-A/q  (  L  4°sdp' P 

4^  7oAA0 


(81) 


(82) 


From  (81)  and  (82),  we  can  see  that  if  "a2  ^  2°  ’  then  Pi  >  which  implies  that  Pi  >  P2,  i.e.  the  power 

sd 

assigned  in  the  first  transmission  round  should  be  larger  than  that  for  the  second  retransmission  round.  The 
condition  '"'Y0  >  51  means 

(TZ  ,  ^ 

sd 


70A/0 


1 

>  - 


70A/0 


(7 


sd 


2  ">T 4 


which  is  true  when  po  >  1  —  e  2  k  0.3935.  On  the  other  hand,  if  ^  (i.e.  po  <  1  —  e  2),  then 

sd 

P\  <  which  means  Pi  should  be  less  than  P2  (an  opposite  power  assignment  strategy  compared  to  that  of 
Pi  >  P2).  Especially,  when  the  targeted  outage  probability  is  po  =  1  —  e_2,  the  optimal  power  assignment 
is  Pi  =  P2,  i.e.  an  equal  power  assignment,  no  matter  what  arc  the  required  SNR  70,  the  noise  variance  A/o 
and  the  channel  variance  a2d.  From  the  case  of  L  =  2,  we  can  see  that  the  optimal  power  can  be  assigned 
either  in  an  increasing,  decreasing  or  equal  way  depending  on  the  targeted  outage  probability  performance  of 
the  H-ARQ  protocol.  This  is  different  from  the  case  in  [35]  where  the  optimal  transmission  power  must  be 
increasing  in  every  retransmission. 

For  the  general  case  of  L  >  2,  numerical  results  (shown  in  Figs.  2-4  in  Section  4)  reveal  that  the  optimal 
power  assignment  sequence  can  be  neither  increasing  nor  decreasing.  Actually,  from  the  theorem  we  can  see 
that  when  Pi  the  optimal  transmission  power  P2  is  less  than  Pi  according  to  (78).  On  the  other  hand, 

asd 

when  Pi  >  ^  the  optimal  transmission  power  P2  is  larger  than  Pi.  This  phenomenon  is  fundamentally 

® sd 

different  from  the  case  in  [35]  where  the  optimal  transmission  power  sequence  is  always  increasing. 

We  note  that  if  Pi  >  qlcn  qie  optimal  power  assignment  sequence  in  Theorem  1  is  monotonically 

® sd 

increasing,  i.e.  P\  <  P2  <  •  •  •  <  Pl-  From  (78),  it  is  easy  to  see  that  P2  >  Pi  in  this  case.  For  any 
k  =  3, 4,  •  •  •  ,  L,  since  1  —  e~x  >  x  —  \x2  for  any  x  >  02,  so  from  (79)  we  have 


Pk  > 


olAYtzlPi) 

70A/0 

Pk- 1 

^k—2 


Pk- 1  70-A/q 
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^dEtp.r  (Etpi 


2  /  k- 2  '2 


2^  (Eti  Pi 
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7oM) 


^  7oA/o 
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old^Pi)  ^  ^ 


uzi  p, 

<  1,  it  is  easy  to  see  that 


Pk- 1- 


(S3) 


Pk- 1 

Ek— 2  p 

i= 1  % 


Pk-  170A/0 


>  0. 


(84) 


2In  this  footnote,  we  would  like  to  prove  that  G(x)  =  (1— e~x)— (x—  |x2)  >  0  for  any  a:  >  0.  We  can  see  that  G' {x)  =  e~x+x—  1 
and  G”(*)  =  —e~x  +  1.  Since  G'(0)  =  0  and  G"(x)  >  0  for  any  x  >  0,  so  G'(x)  >  0  for  any  x  >  0,  i.e.  G(x)  is  monotonically 
increasing  for  x  >  0.  With  G(0)  =  0,  we  conclude  that  G(x )  >  0  for  any  x  >  0. 
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Combining  (83)  and  (84),  we  conclude  that  Pf.-  >  Pk-i  for  any  k  =  3,4,  ••  •  ,  L.  Thus,  the  optimal  power 
assignment  sequence  in  Theorem  1  is  monotonically  increasing  in  this  case.  Actually,  in  this  case  the  optimal 
power  assignment  sequence  increases  as  a  function  of  Pi  roughly  in  a  polynomial  way,  which  is  shown  in  the 
next  section. 


3.2.3  Approximation  of  the  Optimal  Power  Sequence 

To  reduce  calculation  complexity  in  the  optimal  power  assignment,  we  present  in  this  section  a  simple  and  tight 
approximation  for  the  optimal  transmission  power  sequence.  The  tight  approximation  allows  us  to  get  more 
insight  understanding  of  the  optimal  power  assignment  strategy  for  the  H-ARQ  protocol. 

Since  1  —  e~x  ~  x  for  small  x,  so  for  any  k  =  3, 4,  •  •  •  ,  L,  the  optimal  transmission  power  P}-  in  (79)  can  be 
approximated  as 
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The  approximation  in  (85)  is  tight  when  ,, 
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that  when  k  =  3, 
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is  small  and  it  is  true  in  general.  We  can  verify 
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70A0  1 

which  is  strictly  less  than  1,  and  it  becomes  smaller  when  the  power  P\  is  larger.  When  k  >  3, 


Pk- iToA/q 
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which  is  small  when  any  of  Pi,  Pi,  ■  ■  ■  ,  Pf;-2  is  large. 

a2  P2 

When  k  =  3,  substituting  P2  =  into  (86),  we  have 
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P3  ~  P2  T  "77“  = 
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When  k  =  4,  substituting  P2  = 


_  gsVi2 


70A0 


and  the  above  approximation  of  P3  into  (86),  we  can  approximate  P4  as 
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Assume  that  for  any  k  <  ko ( >  2),  it  is  true  that 
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then  for  k  =  ko  +  1,  we  have 
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Table  2:  Algorithm  to  determine  the  approximation  of  the  optimal  power  sequence  P/,.  (1  <  k  <  L) 


Step  1  :  Input  7o,pp,  <7^,  A/p.  Calculate  Pp  =  ■ 

sa  1 1—  pq 

Step  2  :  Set  lower  =  0  and  upper  =  Pp. 


Step  3  :  Let  P\  =  (lower  +  upper)/ 2,  and  calculate 
temp  =Pi(l+  -  Pq. 

Check  if  abs(Temp )  <  e(=  0.001),  then  output 
the  optimal  power  P\  and  go  to  Step  4;  otherwise 
if  temp  <  0,  set  lower  =  P\\ 
if  temp  >  0,  set  upper  =  Pp, 
and  repeat  Step  3. 


Step  4  :  Calculate  the  optimal  transmission  power  P&  as  follows  : 

k  =  2,3,- 


Pk 


a2  P 2 
asd\l 

70A/0 


1  4-  ^PiY 
'  to A/o  ) 


,  L. 


i.e.  the  result  in  (89)  is  also  true  for  k  = 
2, 3,  •  •  •  ,  L,  we  have 


Pk 


ho  +  1.  Thus,  by  induction  we  can  conclude  that  for  any  k  = 


ojjP?  (x  |  -2sdPAk~\ 
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(93) 


Based  on  the  approximation  and  the  sum-power  constraint  in  (80),  we  have  a  constraint  on  the  optimal  power 
Pi  as  follows 


(94) 


or  equivalently 


Pi 


<dPi 

7oA/'0 


(95) 


The  left-hand  side  of  the  equation  (95)  is  an  increasing  function  in  terms  of  power  Pi,  so  there  is  a  unique 
solution  for  the  equation.  Thus,  the  optimal  power  Pi  can  be  easily  determined  based  on  the  equation  in  (95) 
by  using  the  Newton  method.  We  summarize  the  above  discussion  in  the  form  of  the  following  theorem. 


Theorem  2  In  the  H-ARQ  transmission  protocol,  the  optimal  transmission  power  at  each  round  can  be  ap¬ 
proximated  as 
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,  L  where  P\  is  determined  by  the  equation 
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where  Pq  = 
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□ 


From  Theorem  2,  we  observe  that  for  any  k  =  2,  3,  •  •  •  ,  L,  the  optimal  transmission  power  Pf:  can  be  ap¬ 
proximated  as  a  function  of  Pi.  The  optimal  transmission  power  Pi  can  be  directly  determined  by  the  equation 
(97),  then  all  other  optimal  transmission  power  values  /), ,  k  =  2,  3,  •  •  •  ,L  can  be  obtained  immediately  based 
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on  the  closed-form  expression  in  (96).  The  procedure  is  detailed  in  the  algorithm  in  Table  II.  We  can  see  that 
the  calculation  complexity  of  the  algorithm  in  Table  II  is  much  less  than  that  of  the  recursive  algorithm  in  Table 
I.  The  approximation  in  Theorem  2  provides  some  insight  understanding  that  the  optimal  transmission  power  in 
each  retransmission  round  varies  in  term  of  P\  (the  transmission  power  in  the  first  round)  in  a  polynomial  way. 

When  L  =  2,  the  constraint  in  (97)  is  reduced  to  Pi  f  1  +  \  =  Po-  By  solving  the  equation,  we  have 


the  optimal  power  value  for  the  first  transmission  round  as  Pi  = 


2Po 


1+ 


4  ct^jPq 
1  sd  u 

70  W) 


which  matches  with  the  exact 


power  value  given  in  (81).  When  L  >  2,  based  on  the  constraint  in  (97),  the  optimal  transmission  power  Pi 
can  be  bounded  as  follows.  Since  the  geometric  mean  is  not  greater  than  the  arithmetic  mean,  we  have 
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Therefore,  the  optimal  transmission  power  Pi  is  bounded  as  follows 
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in  which  the  upper  bound  is  tight  when  Po  is  large.  The  difference  between  the  lower  bound  and  the  upper 
bound  is  less  than 

asd 

In  Figs.  10  and  11,  we  show  comparisons  of  the  approximation  of  the  optimal  transmission  power  sequence 
by  Theorem  2  with  the  exact  optimized  power  sequence  by  Theorem  1.  In  these  two  figures,  we  assumed 
that  the  targeted  SNR  is  70  =  10  dB,  the  required  outage  performance  is  po  =  10  a2sd  =  1  and  Mq  =  1. 
The  maximum  number  of  transmission  rounds  is  L  =  3  in  Fig.  10,  and  L  =  5  in  Fig.  11.  We  can  see 
that  the  approximations  of  the  optimal  transmission  power  values  (solid  line  with  V)  match  very  well  with 
those  based  on  exactly  numerical  calculation  (solid  line  with  ’o’).  For  comparison,  we  also  include  in  the 
figures  the  transmission  power  level  of  the  equal-power  assignment  strategy.  We  observe  that  in  the  first  few 
(re)transmission  rounds,  the  optimal  power  assignment  strategy  assigns  significantly  less  transmission  power 
compared  to  the  equal-power  assignment  strategy. 

The  optimum  transmission  power  sequence  is  increasing  in  both  cases  in  Figs.  10  and  11.  Actually,  the 
optimum  transmission  power  sequence  can  be  neither  increasing  or  deceasing,  which  is  shown  in  Fig.  12.  In 
this  case,  the  maximum  number  of  retransmission  rounds  is  L  =  10,  the  targeted  SNR  is  70  =  10  dB  and  the 
required  outage  performance  is  po  =  10  '■  We  can  see  that  the  optimal  power  assignment  is  decreasing  in 
the  first  two  rounds  and  increasing  after  that.  Moreover,  we  observe  that  in  this  case  there  is  a  gap  between 
the  approximations  of  the  optimal  transmission  power  sequence  and  the  exactly  calculated  sequence.  Since 

when  k  =  3,  the  term  ,,  ■> —  =  - ~  0.63  which  is  not  small  enough  in  this  case, 

A )(£{«!  Pi)  i+irA  6 

70-A/0 

so  the  approximation  of  the  exponential  term  in  (85)  is  not  tight.  The  difference  between  the  approximated 
sequence  and  the  exactly  calculated  power  sequence  can  be  more  significant  when  the  the  maximum  number 
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Figure  10:  Transmission  power  sequence  of  the  optimal  power  assignment  strategy  with  L  =  3 
70  =  10  dB,  po  =  10-3. 


Figure  1 1 :  Transmission  power  sequence  of  the  optimal  power  assignment  strategy  with  L  =  5 
7o  =  10  dB,  po  =  10-3. 
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Figure  12:  Transmission  power  sequence  of  the  optimal  assignment  strategy  with  L  =  10,  70  = 
lOdB.po  =  10"1. 


of  retransmission  rounds  L  goes  to  infinity.  Fortunately,  this  is  not  the  case  in  practice  where  a  reasonable  L  is 
normally  less  than  10  due  to  delay  consideration. 

Finally,  based  on  the  approximation  of  the  optimal  transmission  power  sequence,  the  average  total  transmis¬ 
sion  power  of  the  H-ARQ  protocol  accounting  for  requested  retransmissions  can  be  approximated  as  follows 
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(103) 


The  approximation  of  the  average  total  transmission  power  is  a  function  of  Pi.  Moreover,  when  Pi  >  , 

O sd 

then  using  the  approximation  1  —  e~x  ~  x  for  small  x,  the  average  total  transmission  power  across  requested 
retransmissions  can  be  further  approximated  as 
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(104) 


This  approximation  is  tight  when  L  is  small,  since  Pi  is  normally  larger  than  7o^°  in  this  case.  It  shows  that 

® sd 

the  average  total  transmission  power  is  roughly  the  product  of  the  transmission  power  in  the  first  round  and  the 
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number  of  retransmission  rounds  allowed  in  the  H-ARQ  protocol.  If  we  approximate  Pi  as 
based  on  (101),  then  the  average  total  transmission  power  with  the  optimal  power  assignment  sequence  can  be 
approximated  as 
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3.2.4  Performance  Comparisons  between  the  Equal  and  Optimal  Power  Assignments 


In  this  subsection,  we  compare  the  power  efficiency  of  the  H-ARQ  protocols  with  the  optimal  power  assignment 
strategy  derived  in  this  work  and  the  conventional  equal-power  assignment  approach.  In  numerical  calculation, 
we  assume  that  the  variance  of  the  channel  hsd  is  a2sd  =  1  and  the  noise  variance  is  Mq  =  1. 

For  a  targeted  outage  probability  po,  according  to  (61),  the  equal-power  assignment  approach  should  also 
follow  the  power  constraint 
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Thus,  the  equal-power  assignment  is  Pj  =  P$/L  for  each  l  =  1, 2,  •  ■  •  ,  L.  The  corresponding  average  total 
transmission  power  is 
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When  P0  is  large,  the  average  total  transmission  power  of  the  equal-power  assignment  strategy  can  be  approxi¬ 
mated  as 
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Therefore,  the  power  efficiency  of  the  optimal  power  assignment  strategy  compared  to  the  equal-power  assign¬ 
ment  approach  can  be  quantified  by  the  following  ratio 
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According  to  Inyz^  ~  po,  the  ratio  can  be  approximated  as 

Pegu  ,-y  Pp_  (  1  ,  1 

Popt  L  y  Lp0  ^  l  -  1 J 


(109) 


(HO) 


For  small  targeted  outage  probability  po,  the  ratio  can  be  further  approximated  as  ~  P(^2  .  We  can  see  that 
the  smaller  the  targeted  outage  probability  po,  the  more  power  saving  the  optimal  power  assignment  strategy 
compared  to  the  equal-power  assignment  strategy. 

For  different  targeted  SNR  70  (from  0  dB  to  40  dB),  we  compare  the  average  total  transmission  power  of 
the  optimal  power  assignment  strategy  and  the  equal-power  assignment  scheme  in  Figs.  13,  14,  and  15  for  the 
cases  of  L  =  2,  L  =  3,  and  L  =  5,  respectively.  The  required  outage  performance  of  the  H-ARQ  protocol  is 
set  at  po  =  10~3.  When  L  =  2,  from  Fig.  13  we  observe  that  the  optimal  power  assignment  saves  about  9 
dB  in  average  total  transmission  power  compared  to  the  equal-power  H-ARQ.  When  L  =  3,  we  can  see  from 
Fig.  14  that  the  optimal  power  assignment  shows  about  10  dB  gain  compared  to  the  equal-power  assignment 
scheme.  When  L  =  5,  Fig.  15  shows  that  the  optimal  power  assignment  strategy  significantly  outperforms  the 
equal -power  assignment  scheme  with  a  performance  improvement  of  about  1 1  dB.  Moreover,  it  is  interesting  to 
observe  that  in  each  figure,  the  performance  gain  of  the  optimal  power  assignment  strategy  is  almost  constant 
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Figure  13:  Comparisons  of  the  average  total  transmission  power  between  the  equal  and  optimal 
power  assignment  strategies  with  different  targeted  SNRs.  L  =  2,  po  =  HP  3. 


Figure  14:  Comparisons  of  the  average  total  transmission  power  of  the  equal  and  optimal  power 
assignment  strategies  with  different  targeted  SNRs.  L  =  3,  po  =  10”3. 
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Figure  15:  Comparisons  of  the  average  total  transmission  power  of  the  equal  and  optimal  power 
assignment  strategies  with  different  targeted  SNRs.  L  =  5,  po  =  1CP3. 


Figure  16:  Comparisons  of  the  average  total  transmission  power  of  the  equal  and  optimal  power 
assignment  strategies  with  different  targeted  outage  probabilities.  L  =  2,  70  =  10  dB. 
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Figure  17:  Comparisons  of  the  average  total  transmission  power  of  the  equal  and  optimal  power 
assignment  strategies  with  different  targeted  outage  probabilities.  L  =  3,  70  =  10  dB. 


Figure  18:  Comparisons  of  the  average  total  transmission  power  of  the  equal  and  optimal  power 
assignment  strategies  with  different  targeted  outage  probabilities.  L  =  5,  70  =  10  dB. 
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p 

for  different  targeted  SNR  70  (from  0  dB  to  40  dB).  This  is  consistent  with  the  theoretical  approximation 

^  Opt 

in  (110)  which  does  not  rely  on  the  targeted  SNR  70.  When  L  =  2  and  po  =  10-3,  the  ratio  in  (110)  is 
=  8.99  dB  (the  observed  power  saving  in  Fig.  13  is  9  dB).  When  L  =  3  and  /y,  =  10-3,  the  ratio  in  (1 10) 

p 

is  =  10.48  dB  (the  observed  power  saving  in  Fig.  14  is  10  dB). 

We  also  compare  the  average  total  transmission  power  required  in  the  two  power  assignment  strategies  with 
different  targeted  outage  probability  values.  We  assume  the  required  SNR  is  70  =  10  dB.  Figs.  16,  17  and  18 
present  comparison  results  for  the  cases  of  L  =  2,  L  =  3,  and  L  =  5,  respectively.  From  the  three  figures, 
we  can  see  that  for  an  outage  performance  of  po  =  10-4,  the  power  savings  of  the  optimal  power  assignment 
strategy  compared  to  the  equal-power  assignment  scheme  arc  15  dB  when  L  =  2,  17  dB  when  L  =  3,  and 
19  dB  when  L  =  5.  The  lower  the  required  outage  probability,  the  more  important  optimization  of  the  power 
sequence  becomes.  Moreover,  we  also  observe  that  with  the  same  targeted  outage  performance,  the  larger  the 
number  of  retransmission  rounds  allowed  in  the  H-ARQ  protocol,  the  larger  the  performance  gain  between  the 
optimal  power  assignment  scheme  and  the  equal-power  assignment  scheme. 

We  also  show  the  average  total  transmission  power  based  on  the  approximated  optimal  power  sequence  in  the 
six  figures.  We  can  see  that  the  average  total  transmission  power  based  on  the  approximated  power  sequence 
matches  matches  tightly  with  that  from  the  exact  optimal  power  sequence  in  each  case.  Note  that  the  exact 
calculation  of  the  optimal  transmission  power  sequence  is  based  on  Theorem  1  and  the  approximated  power 
sequence  comes  from  Theorem  2. 

3.3  Differential  Amplify-and-Forward  Relaying:  Performance  Analysis  and  Power  Optimiza¬ 
tion 

In  this  subsection,  first  we  briefly  review  the  DAF  relaying  scheme.  Second,  we  present  our  exact  outage 
probability  result  and  the  asymptotically  tight  approximation.  Third,  the  asymptotically  optimum  power  allo¬ 
cation  scheme  is  presented.  At  the  end  of  this  subsection,  numerical  and  simulation  results  are  given  and  a  few 
conclusions  are  made. 

3.3.1  System  Model 

We  consider  a  system  model  which  consists  of  one  source,  one  relay,  and  one  destination  using  a  DAF  relay¬ 
ing  protocol,  as  shown  in  Figure  1.  In  DAF  transmissions,  each  information  symbol  vrn  is  first  differentially 
encoded  at  the  source  as 


xT  =  vmxT  1  (111) 

where  xT  denoted  the  signal  to  be  transmitted  from  the  source  at  time  r,  vrn  =  ej4'"' ,  and  {(pm}^  '  is  a  set  °f 
M  information  phases  [44,  45].  We  consider  a  cooperative  strategy  with  two  phases,  in  which  users  transmit 
their  information  through  orthogonal  channels  by  frequency  division  or  time  division  multiplexing. 

In  phase  1,  the  source  broadcasts  the  signal  xT  to  the  destination  and  the  relay.  Then  the  received  signals  at 
the  destination  and  at  the  relay  can  be  modeled,  respectively,  as 

VTs,d  =  \/P~ihs)dxT  +  Vld  (! 12) 

yj. r  =  \Z~P\hs,rxT  +  p^r  (1 13) 

where  Pi  is  the  transmitted  power  at  the  source,  hs  ^  and  hs,r  are  the  channel  coefficients  from  the  source  to 
the  destination  and  to  the  relay  respectively,  7/J d  and  pTsr  are  additive  noise,  respectively. 

In  Phase  2,  the  relay  amplifies  the  received  signal  and  forwards  it  to  the  destination  with  an  amplification 
factor  a.  The  received  signal  at  the  destination  can  be  written  as 

Vr,d  =  \fP2hr,dayTsr  +  rild  (114) 
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Figure  19:  An  illustrative  cooperative  system  model. 


where  P2  is  the  transmitted  power  at  the  relay,  hrJi  is  the  channel  coefficient  from  the  relay  to  the  destination, 
and  r]^d  is  additive  noise. 

The  channel  coefficients  hStd,  hSiT  and  hrJi  arc  modeled  as  independent  zero-mean  complex  Gaussian  ran¬ 
dom  variables  with  variances  5^ d,  hj:  r  and  5%d,  respectively.  The  noise  r/(rf,  77  Jr  and  rj^d  are  modeled  as 
independent  complex  Gaussian  random  variables  with  zero  mean  and  variance  Nq.  Note  that  /i, | 2 ,  (i.  j)  E 
{(s,d),  (s,r),  (r,  d) }  are  exponential  random  variables  with  PDF  p\h.  12(f)  =  ^4- exp( — -A- ) .  To  normalize 

the  transmitted  power  at  the  relay,  we  consider 
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From  ( 1 13)-(1 15),  the  received  signal  at  the  destination  in  Phase  2  can  be  rewritten  as 


Vr,d  =  hr,dXT  +  fj^a 


(116) 


where 


hr,d  — 


PlP2 
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' s,r 


Vr’d  \  Pl8g.r  +  Nq 


K,dVTs,r  +  Vl,d 


(117) 

(118) 


and  f/^d  is  a  zero-mean  complex  Gaussian  random  variable  with  variance  ( p^r'+No  +  1)  Ajj. 

For  differential  demodulation  at  the  receiver  side,  the  fading  coefficients  are  assumed  to  be  approximately  the 
same  over  two  consecutive  symbol  periods.  Then,  according  to  (111)  and  (112),  we  may  rewrite  the  received 
signal  at  the  destination  in  Phase  1  as 


yTs,d  =  yTs/vm  +  nTS)d  (ii9) 

where  nTs  d  =  ifsd  —  Vmrfsd  •  Similarly,  with  (111)  and  (116),  we  have 

Vr,d  =  yTr~dlvm  +  K,d  (120) 

where  nTrd  =  fj^d~  v™tfrd-  Finally,  at  the  destination,  we  combine  the  received  signal  from  the  source  in 
Phase  1  and  that  from  the  relay  in  Phase  2  to  jointly  detect  the  transmitted  information.  Based  on  the  maximum 
ratio  combining  (MRC)  detection  [46],  the  combined  output  at  the  destination  is 


T— 1\ 


=  °i  (yB,d 


fyTs,d  +  a2{yTridlTyTr,d 


(121) 
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where  the  coefficients  ai  and  a2  arc  given  by  [42,  37] 


0l  N0 


a2  = 


PlSlr  +  N0 


(122) 


N0(PiSlr  +  P2\hr,d\2  +  N0) 

which  maximize  the  SNR  of  the  combined  output.  The  transmitted  information  symbol  is  detected  as  follows 

vm  =  argmaxrn=0,i,...,M-iRe{v^nyT}  (123) 


Note  that  the  optimum  combining  coefficient  in  (122)  needs  the  instantaneous  amplitude  square  of  the  channel 
hrtd-  In  practice,  this  may  be  replaced  by  the  statistics  of  the  channel  such  as  the  variance  of  the  channel,  which 
results  in  the  following  combining  coefficients  [42] 


at 


a2 


1 

No 

Pl6lr  +  No 

No(Pi6lr  +  P2S2d  +  No) 


(124) 


It  was  observed  in  [42]  that  the  performances  of  the  two  combining  coefficients  (122)  and  (124)  arc  close.  Since 
the  analysis  of  the  outage  probability  based  on  the  combining  coefficient  (124)  is  intractable  which  may  not 
provide  any  insight  understanding,  we  consider  in  the  following  the  analysis  based  on  the  combining  coefficient 
(122).  The  analysis  will  provide  a  benchmark  for  the  performance  of  the  DAF  relaying  scheme  and  reveal  some 
insight  understanding  of  the  DAF  relaying. 


3.3.2  Outage  Analysis  for  The  DAF  Relaying 

In  the  DAF  relaying  scheme,  the  instantaneous  SNR  7  of  the  combining  output  at  the  destination  based  on  the 
MRC  combining  coefficients  (122)  can  be  represented  as  [42,  37] 

7  =  7i  +  72,  (125) 


where 
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An  outage  is  declared  if  the  instantaneous  SNR  7  falls  below  a  certain  specified  threshold  7 th,  so  the  outage 
probability  is  Pout  =  Pr{ 0  <  7  <  77,. }  [47].  If  we  denote  the  CDF  of  the  SNR  7  as  then  the  outage 

probability  of  the  DAF  relaying  is  Pout  =  77.  (77,  ). 

We  develop  the  CDF  of  the  SNR  7  as  follows.  Note  that  71  and  72  arc  two  independent  random  variables. 
Hence  the  CDF  of  7  is 
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(128) 


Let  (3s4  =  Trrh,(3s,r  =  jrp~  and  fir/1  =  jM-,  then  p7l(7i),  the  PDF  of  71,  in  (128)  is  given  by  p7l(7i)  = 

*-us,d  1  ^ur,d 

fSS)d^~^s’d'11  ■  Furthermore,  if  we  denote  'ys^r  =  P\\hS}r\2 /Nq  and  7^  =  p2\hr,d\2 / Nh  then  the  CDF  of  72  can 
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be  determined  as 


F7I(x)=Pr{ 0  <  7;TT  <  x} 
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Substitute  (129)  with  argument  x  —  71  into  (128),  we  have 

F7(a;)  =  J  Ps,de  ^s,djld'yi 

-J[jr,de-0’’Ax-'11 }  1 1+^  ^  d^de-P^d^ . 


(129) 


(130) 


In  (130),  the  first  paid  is  f*  f3s,d^~^s’dlld'y\  =  1  —  e~@s’dX.  The  second  paid  can  be  further  calculated  by 
switching  the  integration  order,  i.e.  integrating  in  terms  of  the  variable  71  first  and  then  for  the  variable  t,  as 
follows: 
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By  combining  the  above  two  parts,  we  obtain  the  CDF  of  the  SNR  7.  Therefore,  the  outage  probability  of  the 
DAF  relaying  can  be  specifically  given  by 


Poutilth)  =  1  -  e  &s’dlth  -  [ 
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(131) 


Compaidng  to  the  result  in  [43],  in  which  the  outage  probability  was  presented  by  a  triple  integral,  our  result 
involves  only  a  single  integral,  which  can  be  easier  to  calculate  numerically. 

For  a  special  case  of  (3S^  =  fts.r,  i.e.  the  source-destination  and  source-relay  links  have  the  same  channel 
variance,  the  outage  probability  in  (131)  can  be  simplified  as 
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which  can  be  further  represented  as 


Pout{lth)=l  -  e~^h  -  e-^h— -A 


(1  +  (3s,r)f3r,d 
+e-/3s,r7 th  2j3s^thK2(2^J (1  +  Ps,r)Pr,dlth), 


(132) 


(133) 


where  K2{-)  is  the  2nd  order  modified  Bessel  function  of  the  second  kind  [48,  (3.478.4)]. 

Although  the  exact  expression  in  (131)  is  much  simpler  than  the  result  in  [43],  it  still  involves  an  integral. 
To  get  further  insight  understanding  of  the  DAF  relaying  scheme,  we  try  to  approximate  the  outage  proba¬ 
bility  in  the  following.  With  the  aid  of  an  approximation  of  the  modified  Bessel  function,  we  can  obtain  an 
asymptotically  tight  approximation  for  the  outage  probability  at  high  SNR.  We  have  the  following  result. 
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Theorem  3  If  all  channel  links  hs td,  hSjT  and  hr^d  are  available,  i.e.  5S:d  /  0,  5StV  /  0,  and  5r:d  /  0,  then 

for  sufficiently  high  SNR,  the  outage  probability  Pout  of  the  DAF  relaying  can  be  approximated  as 
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Proof:  First,  we  try  to  approximate  the  CDF  of  72  in  (129)  as  follows.  Denote  t  =  0r,d,lr,d  and  y  = 

(1  +  0s,r)0r,dX,  then  the  CDF  of  72  (129)  can  be  rewritten  as 


For  simplicity,  denote  g(y) 
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Jo 
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where  Kq(-)  is  the  zero-order  modified  Bessel  function  of  the  second  kind.  Note  that  for  sufficiently  high  SNR, 
y  =  (1  +  0s,r)0r,dX  is  small,  which  may  go  to  zero.  For  small  y,  there  exists  a  positive  value  yo  (0  <  yo  <  y) 
such  that 


g(y)  ~  g(yo)  +  g'(yo)y, 


where  for 


small  yo,  g(y0 )  and  g'(y0)  c^n  be  approve 

g'{yo)  =  - 


imated  as  [49] 

2Ko{2y/yf)  «  ln(y0). 


Thus,  for  sufficiently  small  vaiiable  y,  we  have 

g(y)  ~  1  +  ln(y0)y  ~  1  +  ln(y)y.  (135) 

With  the  approximation  in  (135)  and  e  ^s’rX  ~  1  —  0s,rX,  we  have  an  approximation  for  the  CDF  of  72  in  (129) 
as  follows: 


Fl2  (x)  ~  0s, rx  -  0r,dxln(0r,dx) .  (136) 

Substituting  the  approximation  (136)  into  (128),  and  with  cT^sMn  ~  1  —  0SJf)\  for  sufficiently  small  argument 
0s, dll,  the  CDF  of  7,  Fcy(x),  can  be  approximated  as 

rx 

Fcy{x)  «  /  [0s,r(x  ~  7l)  -  0r,d{x  -  71  )ln(0r}d(x  -  71))] 

Jo 

X  [0s, d(f  -  0s,dll)\dll 

=  h+h  +  h  +  h,  (137) 


where 


h 

h 

h 


[x  1 

0s,d0s,r  (X  —  ll)dji  =  ~0s,d0s,rX2, 

[X  1 

~0s,d0s,r  Jv  (X  ~  7l)7l dll  =  -Q0s,d0s,rX3, 

-0s,d0r,d  /  (x  -  li)ln{/3r,d{x  -  7i))d7i 

Jo 

y^0s,d0r,dX  2  0s,d0r,dX  ln[(3rdX'), 

I4  =  0sd0r,d  (x  -  ll)llln(Pr,d{x  -  7l))d7l 
Jo 

5  1 

=  -7^0s,d0r,dX3  +  - 0ld0r:dX3ln(0r,dX )■ 


(138) 
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Substituting  7i,  I3  and  I4  into  (137),  we  have 

T^y(x)  —  2  T  -^(3s,d{3r,d'E  ^ (3s,dfir,d% 

-\p2s,d,Ps, rx3-  ^0ldPr,dX3+  \f32s  '  d(3r,dxsln{(3r4x).  (139) 

We  observe  that  when  fts,d-  Ts,r  and  /3r^  go  to  zero,  the  terms  involving  d2  dBs,r  and  (3$  df3r,d  g°  to  zero  faster 
than  the  terms  with  (3s^(3s,r  and  /3s,dPr,d-  As  a  result,  for  high  SNR,  the  CDF  of  7  can  be  approximated  as 

Fry{x)  ~  ~  Ps,df3s,r%  T  ~^f3s,dftrld%  ^  Ps,df3r,d%  l,ft'(Pr,d%')  5  (140) 

which  implies  the  approximation  of  the  outage  probability  in  the  theorem. 


3.3.3  Optimum  Power  Assignment 


Based  on  the  closed- form  tight  approximation  of  the  outage  probability  in  (134),  we  can  determine  an  asymptot¬ 
ically  optimum  power  allocation  scheme  for  the  DAF  cooperation  scheme.  For  a  given  total  transmitted  power 
budget  of  P\  +  P2  =  P,  we  tty  to  minimize  the  outage  approximation  to  find  an  optimum  power  allocation 
ratio. 

Let  Pi  =  rP.  P2  =  (1  —  r)P,  and  denote  B  =  \  —  ln('p'^" ),  then  the  outage  approximation  (134)  can  be 

~  r,c£ 

rewritten  as 


2P2S2  \r2S2r  '  r(l-  r)b2 


Pont  - 


1  ^  B  +  ln{  1  —  r) 


r.d 


To  find  an  optimum  power  ratio  r,  it  is  equivalent  to  minimize 

A  1  B  +  ln(l-r) 
r2b2r  r{  1  -  r)b2rd 

By  taking  derivative  of  (142)  with  respect  to  r,  we  have 

df(r)  _  2  B(  1  -  2r) 

dr  r3$s,r  ^r/2(^-r)2 

1  ln{  1  —  r)(l  —  2  r) 


S2dr( 1  -  r)2  52dr2(  1  -  r)2 


r,d 


(141) 


(142) 


(143) 


Set  the  derivative  (143)  equal  to  0,  we  have 

2(i  -r)2  +  %^(1  -2r)r 

°r,d 

b2  b2 

+  -py~r2  4-  1  —  r)(l  —  2r)?’  =  0.  (144) 

°r,d  K,d 

For  simplicity,  denote  u(r )  =  ln{  1  —  r)(l  —  2r)r,  then  it(|)  =  0,u'(A)  =  ln{2)  and  u"{\)  =  Aln{2)  +  4. 
Therefore,  the  first  three  terms  of  the  Taylor  expansion  to  approximate  u(r)  are 

ln{l  -r)(l  -2r)r«u(^)  +  u'(^)(r-  +  ^u"(^)(r-  ^)2 

=  (2Zn(2)  +  2)r2  -  (ln( 2)  +  2)r  +  (145) 

By  plugging  (145)  into  (144),  we  obtain  the  following  equation 

C^r2  +  C\r  +  Co  =  0,  (146) 

in  which  Co  =  -  +  2,  Ci  =  S^L{—ln{2)  —  2  +  B)  —  4,  and  C2  =  wr(2(n(2)  +  3  —  2B)  +  2.  Solving  the 

r,d  r,d  ur,d 

equation,  we  obtain  the  following  power  allocation  result. 
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Table  3:  Comparing  the  optimum  power  ratio  by  exhaustive  search  and  approximation, 

assuming  {S2sd,  82sr,  82rd}  =  {1,1,1} 


P  in  dB 

Exhaustive  search  r 

Approximated  r 

25 

0.6746 

0.6662 

30 

0.6296 

0.6272 

35 

0.6052 

0.6041 

40 

0.5891 

0.5886 

Table  4:  Comparing  the  optimum  power  ratio  by  exhaustive  search  and  approximation, 

assuming  {52d,  b2sr,  82rd}  =  {1, 1, 100} 


P  in  dB 

iSs.dV  Ss,r>  ^r.d}  =  {1,1,100} 

Exhaustive  search  r 

Approximated  r 

25 

0.8690 

0.8607 

30 

0.8580 

0.8513 

35 

0.8485 

0.8429 

40 

0.8401 

0.8354 

Theorem  4  If  all  channel  links  hsd,  hS)T  and  hr^  are  available,  i.e.  8S)d  0,  8s^r  0,  and  8T)d  0,  then 

for  sufficiently  high  SNR,  the  optimum  power  allocation  ratio  r  =  P\/ P  for  the  DAF  relaying  is 


(Cr,d—  2)  S2r—  452d  +-\J,C2d—2)  ^r+  (8Crtd— 12)  8gr82d 
(4CM-6)52r-4b;2d 


(147) 


where  Cr  d  =  —  ln(2p2Jth)  +  i.  Furthermore,  when  SNR  goes  to  infinity,  Cr  d  is  sufficiently  large  and  then  the 

^ °r,d  ’ 

ratio  (147)  can  be  further  approximated  as: 


r 


1  + 


Kd 

51,rCr,d 


4 


(148) 


We  list  in  Tables  3  and  4  the  comparisons  of  the  optimum  power  ratio  r  based  on  the  closed-form  approx¬ 
imation  presented  in  Theorem  4  and  the  ratio  obtained  by  exhaustive  search  based  on  the  exact  expression 
in  (142).  The  comparison  results  are  shown  for  various  total  transmitted  power  budget  P,  and  for  variances 
{b2d,  8f  r ,  82 d }  =  {1, 1, 1}  and  {1, 1, 100},  respectively.  From  the  tables,  we  can  see  that  the  power  ratios 
calculated  based  on  Theorem  4  arc  very  close  to  those  by  exhaustive  search. 

From  Theorem  4,  we  observe  that  that  the  for  high  SNR,  the  asymptotically  optimum  power  ratio  does  not 
depend  on  the  source-destination  link,  and  it  depends  only  on  the  variances  of  the  source -relay  and  relay- 
destination  links.  Moreover,  as  the  total  transmitted  power  P  goes  to  infinity,  the  optimum  ratio  r  goes  to  0.5, 
which  implies  that  an  equal  power  allocation  scheme  (Pi  =  P>  =  P/2)  is  optimum  in  this  case. 


3.3.4  Simulation  Results  and  Discussions 

To  verify  our  analytical  results,  we  carried  out  some  computer  simulations.  In  our  simulations,  we  assume  that 
the  variance  of  the  noise  is  1  (i.e.  Nq  =  1),  and  the  outage  SNR  threshold  is  yth  =  10  dB. 

Figures  20-22  depict  the  simulated  outage  probability  performance  curves,  the  outage  probability  approxi¬ 
mation  based  on  Theorem  1,  and  numerical  calculation  based  on  the  close-form  expression  in  (131)  for  three 
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-  New  exact  formula 

-  Closed-form  approximatio 

-  Simulation 


P/N0  (dB) 


igure  22:  Outage  probability  of  the  DAF  relaying  with  variances  {8Z  d,  8Z  r .  5zd} 


DAF  relaying  with  { 5^  8^  sf  }={1  .1  ,1  } 


0.1  0.2  0.3  0.4 


3:  Numerical  search  of  the  optimum  power  ratio  r  for  the  DAF  relaying  with  vari 
1).  and  different  transmit  nower  P=20dB.  30dB  and  40dB. 


Figure  24:  Performance  comparison  of  the  DAF  relaying  with  the  equal  and  the  optimum  power  allocation, 
with  variances  <5*r,  S2d}  =  {1, 1, 1}. 


Figure  25 :  Numerical  search  of  the  optimum  power  ratio  r  for  the  DAF  relaying  with  variances  {52  d,  52  r,52  d} 
=  {1,1, 100},  and  different  transmit  power  P=20dB,  30dB  and  40dB. 
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Figure  26:  Performance  comparison  of  the  DAF  relaying  with  the  equal  and  the  optimum  power  allocation, 
with  variances  {Sgd,Sgr,Sfd}  =  {1, 1, 100}. 


scenarios: 


=  {1,14},  {Sld,Slr,62rJ  =  { 1 , 1 , 100},  and}^,^  ^}  =  {M00,1} 


Js,di 


Jr,d 


s,d>  us,r  i  ur,d i 


re¬ 


spectively.  We  equally  allocated  power  for  the  source  and  the  relay  in  the  three  scenarios.  We  can  see  that 
the  numerical  results  based  on  the  exact  expression  in  (131)  matches  with  the  simulation  results.  The  outage 
probability  approximation  is  loose  at  low  SNR  and  it  merges  with  the  simulation  curve  around  10"  2  which  is 
tight  in  all  the  three  scenarios. 

Figure  23  shows  the  outage  probability  performance  versus  power  ratio  r  =  P\/P  for  different  total  trans¬ 
mitted  power  budgets,  in  which  the  channel  variances  are  f)}  r,  iry/}  =  {1, 1, 1}.  In  this  case,  the  optimum 
power  ratio  is  close  to  1/2.  The  performance  of  the  optimum  power  allocation  scheme  is  close  to  that  of  the 
equal  power  allocation  scheme  in  this  case,  which  is  shown  in  Figure  24. 

Figure  25  shows  the  outage  probability  performance  versus  power  ratio  r  =  P\/P  for  the  case  of  82sr,  <5/  d} 
=  {1, 1, 100}.  In  this  case,  the  optimum  power  ratio  is  around  0.85.  Simulation  results  in  Figure  26  show  that 
the  performance  of  the  optimum  power  allocation  scheme  is  about  2  dB  better  than  that  of  the  equal  power 
allocation  scheme  in  this  case. 


3.4  Joint  Power  Optimization  for  Multi-Source  and  Multi-Destination  Relay  Network 

In  this  subsection,  first  we  introduce  briefly  the  system  model  of  a  multi-source  multi-destination  relay  network 
where  transmissions  occur  over  non-orthogonal,  in  general,  channels.  Second,  we  determine  the  maximum 
ratio  combining  of  the  received  signals  at  each  intended  destination  and  exploit  the  resulting  SINR.  Third, 
we  determine  the  optimum  power  assignment  for  the  sources  and  the  relay  that  minimizes  the  total  power 
consumption  under  the  condition  that  the  SINR  requirement  of  each  source-destination  pair  is  satisfied.  Fourth, 
we  determine  the  optimum  power  assignment  that  maximizes  the  minimum  SINR  among  all  source -destination 
pairs  subject  to  any  given  total  power  budget.  Finally,  numerical  studies  arc  provided  and  some  discussion  arc 
made  at  the  end. 

The  following  notation  is  used  throughout  this  subsection.  Bold  letters  in  uppercase  and  lowercase  denote 
matrices  and  vectors,  respectively.  (•)*,  (■)  r  and  (■)n  represent  the  conjugate,  the  transpose  and  the  Hermitian 
transpose  operation,  respectively.  |  •  |  and  1 1  •  1 1  represent  Euclidean  norm  of  a  complex  number  and  a  vector. 
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Figure  27 :  Multi-source  multi-destination  relay  network. 

respectively.  1^  is  an  L  x  L  identity  matrix.  diag(h \ ,  /12,  •  •  •  ,  hj,)  is  an  L  x  L  diagonal  matrix  with  diagonal 
elements  hi,  /12,  ■  ■  •  ,  hr,.  Aj;  denotes  a  sub-matrix  of  A  obtained  by  deleting  the  kth  column  and  klh  row  of 
A.  If  afc  represents  the  kth  column  of  the  matrix  A,  then  a,,  denotes  the  vector  obtained  after  removing  the  kth 
entry  from  a/-. 


3.4.1  System  Model 

For  illustration  puiposes  and  simplicity  in  presentation,  we  consider  a  single  relay  code  division  multiplexing 
system  with  K  sources  and  K  destinations  as  shown  in  Fig.  27,  where  transmissions  occur  over  non-orthogonal, 
in  general,  channels.  Our  developments  can  be  generalized  to  multiple-relay  systems  and  other  multiplexing 
schemes  in  frequency  and/or  time.  Let  Sk  denote  the  kth  source  and  Dk  the  corresponding  destination,  k  = 
1, 2,  •  •  •  ,  K,  and  let  R  denote  the  relay.  The  relay  forwards  simultaneously  the  signals  received  from  all  sources. 
Let  bf.  denote  the  transmitted  information  symbol  of  the  source  Sk  with  unity  energy,  i.e.  E { \  by-  | 2  }  =  1,  Vfc. 
The  signal  sent  by  the  source  Sk  can  be  expressed  as 


Sfc  —  c  fchfe, 

where  C/c  =  (cl,  c[2')  •  •  •  ,  cj7'  1 ) 7  is  the  code/signature  of  the  source  Sk,  which  is  a  unit-energy  column  vector 

with  length  L.  The  codes/channels  of  different  sources  are,  in  general,  correlated.  Let  pkj  =  c/  c,  denote  the 

cross-correlation  between  codes/channels  k  and  j,  where  pkj  E  [0, 1)  for  k  ^  j,  and  pkk  =  1.  Let  R  =  (pkj) 
denote  the  K  x  K  cross-correlation  matrix,  i.e. 

R  =  (ci,c2,...,c^)t(ci,c2,...,c^).  (149) 


We  consider  the  following  two-phase  amplify-and-forward  relay  strategy  with  L  time  slots  in  each  phase.  In 
Phase  1,  each  source  Sk  transmits  the  signal  s k  with  transmitted  power  P/,:.  Then,  the  received  signals  at  the 
destination  Dk  and  at  the  relay  R  during  the  7-th  (1  <  i  <  L)  time  slot  can  be  modeled,  respectively,  as 


K 


y's,dk  =  J2^hsldkclbl  +  ns4k'1-i-L’ 
1=1 
K 

^  <  L. 

v  *  o  1,1  L  v  0,1  >  — 

1  =  1 


Vs,  r  =  lhsurcll)bl  +ns}r,1<i 


(150) 

(151) 


In  Phase  2,  the  relay  amplifies  the  received  signals  and  forwards  them  to  the  destination  with  an  amplification 
factor  a  and  transmission  power  Pr.  The  received  signal  at  the  destination  Dk  during  the  7-th  (1  <  %  <  L)  time 
slot  can  be  written  as 


y/Prh® 


(0 

r,dkuys,r 


+  n 


(0 

r,dk  > 


i  =  1,2.---  ,L. 


(152) 
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In  (150)-(152),  d  .  h^r,  and  ,  ( l ,  k  =  1, K)  are  the  coefficients  of  the  channels  between  the  source 
Si  and  the  destination  Dk,  between  the  source  Si  and  the  relay  R,  and  between  the  relay  R  and  the  destination 


Dp.,  respectively,  during  the  i-th  (1  <  i  <  L)  time  slot.  n2\  and  nyi,d  represent  noise  at  the  destination  Dp. 


r,dk 


(i) 

during  the  i-th  time  slot  of  Phase  1  and  Phase  2,  correspondingly,  while  ns,r  represents  noise  at  the  relay  R 

(i)  ( i ) 

during  the  i-th  time  slot.  The  channels  hs’ d) ,  K’r,  and  KX  are  assumed  to  be  independent  Gaussian  random 
variables  with  zero-mean  and  variances  cr2  d  ,<7^  r,  and  a2d  ,  respectively.  All  noise  terms  are  assumed  to 
be  independent  Gaussian  random  variables  with  zero-mean  and  variances  a2.  Without  loss  of  generality,  we 
assume  cr2  =  1. 


The  channel  coefficients  in  matrix  format  can  be  given  as  H Sudk  =  'O-  H 


si,r 


diag(h^r,  ■  ■  ■  ,  hifj.)  and  Hr,dk  =  d,iag(h^d  ,  ■  ■  ■  ,  h/^).  Then,  the  received  signals  can  be  expressed  as 
follows: 

K 


ys,dk 

/*  ^  \/Pl.y^sp,dkclJJl  +  ns,dk  j 

1=1 

K 

(153) 

ys,r  = 

^  "  \/ 7)HS(rc/ bp  T  n.  ;- . 

Z=1 

\j Rratl,,dhys  r  T  nr  cpk , 

(154) 

y  >'.dk  ~ 

(155) 

(L)  \T 

Vs,dk)  ’  y s,r 

-  f7/(1)  7/(2)  -  -  -  r,(L)lTandv  -  G/(1)  7/(2)  .. 

—  {ys,r ,  Vs,r ,  ,  Vs,r  )  ana  y r  dk  —  (yrdfc ,  yr^k  j 

’  ^ r,dk  >  ■ 

where  y =  (y^ ,  y$fc ,  • 

In  ( 153)— (1 55),  the  noise  vectors  ns,dk  •  n.s.r  and  n r.dk  have  elements  that  are  independent  Gaussian  random 
variables  with  zero  mean  and  unit  variance.  The  amplification  factor  in  (155)  is  specified  as 


«2  = 


1 


1 


(156) 


where  /3sur  =  c^.E{H^rHSi)r}Q.  By  substituting  (154)  and  (156)  into  (155),  we  obtain 

K 

JV,c(fc  =  ^  ^  \/Pl^si,r^lbl  +  \/ Pratlr,dk Psvr  T  ^r,dk  ■ 


(157) 


i=i 


The  destination  Dp.  combines  the  signals  received  from  the  sources  in  Phase  1  and  the  signals  received  from 
the  relay  in  Phase  2  to  jointly  detect  the  information  symbol  transmitted  by  the  source  Sk.  The  combined 
received  signal  from  Phase  1  and  Phase  2  at  destination  Dp,  can  be  expressed  in  vector  form  as  follows: 


yk  t  ^ 

yr,dk 


I< 


+  nk, 

1=1, tyk 


where 


H,)fc  = 


ypiKsi,dk  \ 

y/ i ) /'  nil iS .  rH r,dk  / 


is  a  2 L  x  L  virtual  channel  matrix  from  the  source  Si  to  the  destination  Dp,,  and 


n  k 


ns,dk 

\fPrOtW.r,dkns,r  T  ^r.dk 


is  an  equivalent  noise  vector  of  length  2 L.  We  note  that  H kk  is  the  channel  matrix  associated  with  the  desired 
source  Sp„  while  H//;.,  (I  /  k)  arc  the  channel  matrices  of  the  interfering  sources.  Based  on  maximum  ratio 
combining  (MRC)  detection  [?],  the  transmitted  signal  from  the  source  Sp,  is  generally  detected  as, 

bk  =  argminbk&A\wk  oyk  -  bk\2, 
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where  A  is  the  transmitted  signal  set.  For  example. 


bk  =  sign(Re{  w^0yj), 


for  BPSK  and 


h  =  sign{Re{  w^yj)  +  sign(Im{w£0yk})j, 


for  QPSK  in  which  j  =  y— 1.  The  combining  weight  vector  wk  0  of  size  2  L  is  chosen  to  maximize  the  SINR 
at  the  destination  Dk,  which  is  given  by 


SINR(wfc) 


_ E{\w^Uk)kckbk\2} _ 

£{|wf  (EiLi&k  H  i,kdbi  +  nfc)l2}  ’ 


(158) 


i.e. 


Wfc,0  =  argm,cixWk  SINR(wfc). 

Note  that  in  (158),  the  expected  value  in  numerator  is  taken  over  the  random  variable  bk,  while  the  expected 
value  in  denominator  is  taken  over  the  random  variables  bi,l  ^  k  and  all  independent  noise  terms  in  n/,.. 


3.4.2  System  Performance  Analysis 

To  determine  the  maximum  SINR  weight  vector  w/,..0  and  the  corresponding  SINR  at  the  destination  Dk,  we 
first  define 


I\ 

Ufc  =  +  Tfc, 

1=1, tyk 


Tk  =  E{nknj? } 


I L  0 

o  pra2 nr4kn*dk  +  iL, 


Then  by  taking  the  expectation  over  the  signals  bk  in  numerator,  bf.  l  ^  k  and  noises  in  denominator,  the  SINR 
in  (158)  with  any  given  combining  weight  vector  vvy.  can  be  written  as 


SINR(wfc) 


w^Ufcwfc 


(159) 


It  is  easy  to  check  that  U/,.  is  Hermitian  and  it  can  be  represented  in  terms  of  its  eigenvalues  Ai,  A2,  •••,  A k  and 
their  associated  eigenvectors  q1;  q2, q^  as  follows: 


K 

u  k  =  XMk<lk, 

k= 1 

where  \k  >  0 ,\/k  and  qj^q;  =  0,  Mk,  l,k  /  l.  Let 

=  1". 

fc=i 

1  11 

then  U i  is  Hermitian  and  U/,.  =  U;:  U? .  Moreover, 


k=  1 
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which  is  also  Hermitian.  According  to  Schwartz  inequality,  we  have 


Iw^H^I2  =  |w^U|U7Hfc,fccfc|2 

w^Ufcwfc  w^U|U|wfc 

<  HwfUll|2||U?HMcfc||2 
HwfUflP 
=  ||U“^Hfc,fccfc||2, 

1  _  1 

where  the  equality  holds  when  w^U|  =  (Ufc  Thus,  the  maximum  SINR  weight  vector  Wfc;0  is 

given  by 

Wk,o  —  U^, 


The  corresponding  maximum  SINR  at  the  destination  D};  with  the  optimum  weight  vector  w/.,0  is  equal  to 

SINRfc  =  H^U^H*,***.  (160) 

In  order  to  optimally  allocate  power  to  all  sources,  we  further  exploit  SINR/,.  in  (160)  as  follows.  We  define 
Cfc  =  diag( Ci, c^-i,  Cfc+i, Ck),  which  is  a  L(K  —  1)  x  (K  —  1)  block  diagonal  matrix  formed  by  placing 

all  code  vectors  except  c&  in  the  diagonal  positions,  and  =  [Hi  fc, Hfc_i  which  is  a 

2 L  x  L(K  —  1)  interference  channel  matrix.  Using  the  above  notation,  SINR/,  can  be  expressed  as 

SINRfc  =  c^H^fc(HfcCfcCfc  Hfc  +  rfe)-1Hfc)fcCfc.  (161) 


According  to  the  Woodbury  matrix  inversion  lemma  [75],  we  have 

(Hfccfccf  Hf  +  Tfc)"1  =  T"1  -  r^HfcC^I/c-t  +  cf  Hf  T"1. 

Let  us  now  define  the  following  matrix  =  ( fmn ),  where 

,,,  , - (  a2  Pr\hr  dJ2  h*  rhSr, 

fmn  =  PmnV PmPn  ^sm,dfc^Sn,<4  H  p  |/2  “  |2  _|_  J 


(162) 


and  denote  by  f(/7  the  kth  column  vector  of  the  matrix  ¥'ik'! .  Then,  after  some  algebraic  calculations,  we  can 
see  that 


Ak) 

J  kk 

=  r^H^Cfc 

0 

(163) 

-p(k) 

<  k 

(k) 

where  fj~ ’  contains  the  channels  and  the  cross-coiTelation  between  the  intended  source  ,S'/-  and  the  interfering 

(k) 

sources,  while  F-  contains  the  channels  and  the  cross-correlation  among  interfering  sources.  Based  on  (161)- 
(163),  we  can  represent  SINR/,.  as 


SINRfc 


r(k) 

■I  kk 


k)H 

k 


(Ifr-i  + 


(164) 


where  the  superscript  (-)®  indicates  the  coiTesponding  destination  D/:. 
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(i) 

Note  that  when  all  the  channels  arc  quasi-static  over  an  information  symbol  period  L,  i.e.  when  hg  d  = 

hSn,dk,  h(slr  =  hSn,r  and  h^dk  =  hr/lk  Vi  =  1,  2,  •  ■  •  ,  L  and  Vn,  k  =  1,  2,  •  •  •  ,  K,  then  /f,J  in  (162)  can  be 
reduced  to 


,dfc^Sn,<4 


a  Pr\hr,dk\  kSmrhSn}r\ 
a2Pr\hr4k\2  +  1  )■ 


(165) 


When  there  arc  only  two  source-destination  pairs  with  one  relay,  i.e.  K  =  2,  the  SINR/,.  in  (164)  with  quasi¬ 
static  channels  is  specified  as 


'  kk 

=  Pk 


'33 

ii  2  ,  a2Pr\hr,dk\2\hSk,r\ 
1  Sk4kl  a2Pr\hr>dk\2  +  l 


P^PkPj 


,*  ,  ,  a2^|^A-l2^ 

nsk,dkn^Ak  +  a2pr\hi 


lr,dk  I 


+  1 


a2Pr|h.ridJ2P,|/).Sjir|2 

CX2  Pr  |  hr^dk  | 2  +  1 


where  p  =  pi 2  =  P21  and  the  subscript  j  =  2  if  k  =  1  while  j  =  1  if  k  =  2.  Furthermore,  when  the 
code/channel  vectors  C&  are  orthogonal  to  each  other,  i.e.  pkj  =  c^c j  =  0,  for  j  ^  k ,  then  ’  =  0  for 
quasi-static  channels.  Hence  the  SINR  at  destination  Dk  is  given  by 


SINRfc 


,.(k) 
J  kk 


Pk\hSk,dk  I 


a2Pr\hr4k\2Pk\hSk,r\2 
a2Pr\hrydk  | 2  +  1 


(166) 


which  is  the  sum  of  the  SINRs  of  the  direct  link  and  the  relay  link. 


3.4.3  Optimum  Power  Assignment  Under  SINR  Constraints  for  All  Source-Destination  Pairs 

In  this  subsection,  we  determine  the  optimum  power  assignment  for  the  sources  and  the  relay  that  minimizes  the 
total  power  consumption  under  the  condition  that  the  SINR  requirement  of  each  source-destination  pair  is  satis¬ 
fied.  First,  we  consider  the  power  optimization  for  the  relay  network  in  a  general  setting  where  codes/channels 
of  different  sources  may  have  arbitrary  correlation.  Then,  we  discuss  a  simplified  power  optimization  scheme 
for  a  special  case  where  codes/channels  of  different  sources  arc  orthogonal,  and  provide  an  intuitive  interpreta¬ 
tion  for  the  proposed  scheme. 

Let  us  assume  that  the  SINR  requirement  for  the  source-destination  pair  (Sk,  Pk)  is  7 k,  k  =  1, 2,  •  •  •  ,  K. 
Then,  the  problem  of  optimizing  power  to  minimize  the  total  power  consumption  and  satisfy  all  source- 
destination  SINR  requirements  can  be  formulated  as: 


min p1,...,pK-pr  J2k=i  Pk  +  Pr, 

s  .  t .  SINRfc  >  7 k,  1  <  k  <  K, 


(167) 


where  the  transmission  power  terms  Pi,  P2,  ■  ■  ■  ,Pk  and  Pr  arc  all  non-negative. 

Let  us  define  a  normalized  power  factor  at  the  relay  as  follows 

x  =  a2Pr,  (168) 


where  a  is  the  amplification  factor  specified  in  (156).  The  parameter  x  will  play  a  key  role  in  the  optimization 
procedure.  Let  us  also  denote  by  U :  k '  =  (glnl)  a  matrix  with  elements 


^k'j  (  7  4s  1 

Smn  =  Pmn\P'srn  ,d 


+ 


x\hr  tdk 


1 2h*  rhs 

1  5 1  0 


x\hry 


dk  I 


+  1 


— ),  (m,n=  1,2,---  ,K). 


(169) 
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for  any  m,  n  =  1, 2,  •  •  •  ,K.  Then,  from  (162),  we  can  represent  each  entry  in  Fik>  by 

f(k)  /  p  p  ( k ) 

J  mn  V  1  fn1  n[3rnn' 


It  is  straightforward  to  verify  that 


det(F^)  = 


det(F\ 


(  If  Pi)d<Gl 

1=1, tyk 


From  (164),  we  know  that  \k-i  +  F^  is  invertible.  So,  based  on  the  Schur  complement  formula3,  we  have 
det{  lK+F^) 


=  det 


1  ,  f(k)  Uk)H 

1  “T  Jkk 


Ik- t  +  Ff 


=  det( IK-!  +  Ff )  •  [l  +  fj£  -  f^)H(lK-i  +  F^  j^ff' 
=  det( IK-i  +  F^})  •  (1  +  SINR*), 

where  the  last  equality  follows  from  the  expression  of  SINR/,.  in  (164).  Thus,  we  have 


1  +  SINRfc  = 


det{lK+F^) 
det(IK-i  +  Ff; 


We  note  that  for  moderate  or  high  SINR,  we  can  approximate  1  +  SINRfc  ~  SINR/.  and  1+ ~  f^k\  \/k,l  = 
1,  •  •  •  ,  K  [55] [74].  So,  based  on  (173),  we  can  approximate  SINR/.  as  follows: 


SINRfc 


det{ F<*>)  _  Pkdet( G(fc>) 
det(F^)  ~  det{ Gf)  ' 


Therefore,  the  optimization  problem  in  (167)  can  be  written  as 


lin Pi,-  ,PK;  PrY.k= 1  Pk  +  pr 


Pkdet(  GW) 
det(G P) 


>7 k,  1  <  k  <  K. 


Let  V  denote  the  feasible  set  of  the  optimization  problem  in  (175),  i.e. 

V  =  {Pi,  ■  •  •  ,  PK,  Pr  |  SINRfc  >  7fc,  VI  <  k  <  K}. 


We  may  partition  V  into  disjoint  subsets  such  as 


V=  |J  14 


3  If  matrix  D  is  invertible,  then  [75] 


det  ^  ®  )  =det(D)-det(A-BD^1C). 

\  C  D 
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where 


Vx  =  {Pi,--  -  ,  Pk ,  Pr  |  SINRfc  >  7fc,  VI  <  k  <  K,  a2  Pr  =  x  } 


for  any  x  >  0. 

We  note  that  for  any  given  value  of  the  parameter  x  in  (168),  the  transmission  power  at  the  relay,  Pr,  can  be 
determined  from  (156)  as 


Pr  —  9  —  '  r  'N  Pkdsk .r  xL. 

or  z — ' 


K 


(176) 


k= 1 


Thus,  for  any  given  x  >  0,  the  optimization  problem  in  (175)  over  the  feasible  set  Vx  becomes 

minp1;...  ;Pk  J2k=i(xPsk,r  +  l)pk  +  xL, 
lkdet(G[P) 


(177) 


s  .  t .  Pk> 

We  observe  that  in  (177),  for  any  given  x  >  0, 

j  in  (177)  is 

lkdet{G(-^ 


det(  G(fc)) 
' Ykdet(G[k) ) 


- ,  1  <  k  <  K. 


is  a  constant  which  is  independent  of  Pk  (k  = 


det(  G^) 

1, 2,  •  •  •  ,  K ).  Hence  the  minimal  value  in  (177)  is  obtained  when  all  constraints  hold  with  equality,  i.e. 


Pk  = 


k  =  1,2,  •••  ,K. 


det{G^) 

Then,  the  corresponding  minimal  total  power  of  (177)  is 

A  det{ 

v(x)  = 


(178) 


/\  X - X  UlOl/  VJ7  ) 

w  =  S7l(l's“'  +  1,SV5)+Ii- 


(179) 


which  is  a  function  of  x  >  0.  Let  us  denote  v*  as  the  minimal  value  of  the  objective  function  over  (175)  in  the 
feasible  set  V.  Then,  we  can  see  that 


v*  =  minx>  ov(x).  (180) 

The  above  discussion  shows  that  we  arc  able  to  convert  the  optimization  problem  in  (167)  over  a  multi¬ 
dimension  space  to  the  minimization  problem  in  (180),  which  depends  only  on  one  variable  x  >  0,  i.e.  over 
a  one-dimension  space.  The  minimization  of  v(x)  in  (179)  can  be  easily  solved  by  a  numerical  search  for  the 
optimal  value  of  the  parameter  x  >  0.  With  the  optimal  value  x*  that  minimizes  the  function  v(x)  in  (179),  we 
can  obtain  the  corresponding  optimal  power  P^  and  P*  based  on  (178)  and  (176),  respectively. 

In  the  previous  section,  we  solved  the  optimization  problem  for  a  general  multi-source  multi-destination 
relay  network  with  arbitrary  correlation  among  user  codes.  In  this  subsection,  we  arc  able  to  further  simplify 
the  optimization  when  the  signatures  of  different  sources  arc  orthogonal  and  the  fading  channels  arc  quasi-static 
during  the  transmission  period  of  a  signature  code. 

In  particular,  for  multi-source  multi-destination  relay  networks  with  orthogonal  transmissions,  the  cross¬ 
correlation  matrix  R  in  (149)  is  an  identity  matrix,  so  both  and  in  (174)  arc  diagonal  matrices,  and 

where  is  specified  in  (169).  Under  the  assumption  of  quasi-static  channels,  we  can  write  =  hSk>r, 
}^ldk  =  hSk  A,  and  =  hrtdk,  Vi  =  1,2, ...,  L,  and  \/k  =  1, 2, ...,  K.  We  note  that  (3Sk)r  =  cr2k  r  where 
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<jg  r  is  the  variance  of  the  source-relay  channel  hSk^r,  i.e.  a2k  =  E{\hSk)T\2}.  Thus,  by  substituting  (181) 
into  (179),  we  have 


where 


v(x)  =  Y'  7  k(XC7Sk,r  +  ^ 

V  )  ^  \h  12  ,  4 

k=l  \'hk,dk  I  T"  x\hrAk\2+ 1 


xL 


A' 


fc=l 


x  +  dfc 


ak 


bk 


Ck 


dk 


l^sfc,r|2  /-,  _  \b-sk,dk\  ask,r\ 

Ai\hr4k\2[  A\hr4k\2  > 

L  *Lr 

Klk  A  ’ 

1  \K,r\2  ask,r 

A2\hr,dk\2  ’ 

\bsk4k  | 2 

A\K,dkr 


(182) 


(183) 


and  A  =  \hSk4k\2  +  |/iSfe,r|2. 

We  note  that  for  any  k  =  1, 2,  •  •  ■  ,  K,  if  a/.  >  0,  each  term  Ck  +  bkX  +  x°fd  in  (182)  is  convex  with  respect 
to  x  >  0,  and  it  can  be  minimized  by 


x  =  max(  0,  —  dk  + 


If  cik  <  0,  Ck  +  bkX  +  x°£dk  is  increasing  in  [0,  -Too),  it  implies  that  the  minimum  point  occurs  at  x  =  0.  Let  us 
denote 


where  for  any  A:  =  1,2,---  ,  K 


Xk 


Xmin  —  mi7l{x\ ,  X2,  ... ,  Xk)  > 
Xmax  —  max(xi,X2,  ■■■,  Xk), 


max(  0,  —  dk  + 
0 


if  ak  >  0, 
if  ak  <  0. 


(184) 

(185) 


(186) 


Then,  each  term  Ck  +  bkX  +  x^kd  in  (182)  is  decreasing  in  the  range  [0,  xrmn]  and  increasing  in  the  range 
\xmax ?  +oo),  which  implies  that  the  optimal  solution  x*  that  minimizes  the  function  vix)  in  (182)  is  bounded 
as 

Xmin  X  T  Xmox  •  (187) 

Thus,  to  find  the  optimal  solution  x*,  we  only  need  to  search  within  the  range  \xrnm ,  xmax\.  We  note  that 
xmin  =  0  if  there  exists  k  G  {1,  2,  •  •  •  ,  K}  such  that  ak  <  0. 

From  (186),  we  can  see  that  the  necessary  condition  for  Xk  >  0  is 


which  implies  that 


\hSk,r\2\hr4k  | 


a 


2 

sk,r 


>  \hSk,dk\ 


(188) 
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In  a  non-fading  or  slow-fading  scenario,  where  the  coherence  time  of  the  channel  is  much  longer  than  the  delay 
requirement  of  the  channel.  In  this  situation,  the  channel  gain  remains  roughly  constant  over  all  the  time  [76], 
then  we  can  safely  assume  that  \hij\2  ~  afj,  where  nf  j  is  the  variance  of  the  channel  ht.j.  Then  (188)  is 
equivalent  to 


G 


2 

r,dk 


>  G 


2 

dk  ' 


(189) 


If  none  of  the  inequalities  in  (189)  is  true,  i.e.  G2dk  <  G2k  dk  for  all  k  =  1,  2,  •  •  ■  ,  K,  then  xmax  =  0,  which 
implies  that  the  optimal  solution  x*  =  0.  So,  the  corresponding  optimum  power  at  the  relay  P*  =  0,  which 
means  that  the  relay  is  not  needed.  In  other  words,  if  each  relay-destination  channel  link  is  weaker  than  the 
intended  source-destination  channel  link,  then  we  should  not  use  relay. 

In  addition,  if  a&  <  0,  \/k  =  1,  2,  •  •  •  ,  K,  then  each  term  c/;;  +  bkX  +  ^  in  (182)  is  an  increasing  function 
of  x  G  [0,  Too),  therefore  the  optimal  solution  is  x*  =  0.  Given  that  <  0  is  equivalent  to 


\hSk,dk\  ask,r  —  {\hsk,dk\  +  \hsk,r\  )\hr,dk\  j  (190) 

it  implies  that  if  all  channels  \hSkidk\2,  (k  =  1,  2,  •  ■  •  ,  K )  that  correspond  to  direct  links  arc  strong  enough  or 
all  channels  of  the  relay  links  \hr^dk  |2,  (k  =  1.  2,  •  •  •  ,  K)  arc  weak,  then  avoiding  use  of  the  relay  can  leads  to 
power  savings.  In  other  words,  if  there  is  at  least  one  relay-destination  link  that  is  better  than  the  corresponding 
source-destination  link,  the  relay  should  forward  its  received  signal  to  the  destinations  in  order  to  help  reduce 
the  overall  power  consumption  while  maintaining  at  the  same  time  the  target  performance  level. 

Finally,  for  a  symmetric  system  where  a2si  di  =  g\^  =  ■■■  =  g2k4k1g2u7,  =  of2>r  =  ■  ■  ■  =  o%K>r  and 
(J2dl  =  <x2 d2  =  ' ' '  =  ar  dK>  we  are  able  to  obtain  an  analytical  solution  for  the  minimization  of  v(x).  In  this 
case,  all  dk  arc  equal,  say  d  =  dk,Vk.  Then 

v(x)  =  c  +  bx-\ - ——n  (191) 

x  +  d 


where  c  =  ff/k=i  7 fccfc,  b  =  ffjk=x  kkbk-  a  =  ffk=\  7 Consequently,  the  value  x*  that  minimizes  v(x)  in 
(191)  is 


x  = 


~d+Vk  if  a  >  °> 

0  if  a  <  0. 


(192) 


Substituting  (192)  into  (178)  and  (176),  we  obtain  the  optimal  power  P £  and  P* ,  respectively. 

We  note  that  for  a  symmetric  system,  all  «/,.  arc  equal.  When  the  SINR  requirement  is  high  enough  at 

(^2 

each  destination,  L/(K 7*.)  is  relatively  small  compared  to  “f  ,  so  all  6*7 s  arc  approximately  the  same  ( k  = 
1, 2,  •  •  •  ,  K).  As  a  consequence,  the  ratio  a/6  in  (192)  is  independent  of  7*.  Hence,  the  optimal  value  x*  in 
(192)  is  independent  of  the  SINR  thresholds  7*.  in  this  case. 

In  the  following,  for  comparison  purposes,  we  discuss  the  equal  power  allocation  scheme  where  all  the 
sources  and  the  relay  are  allocated  the  same  power  P.  In  this  case,  the  corresponding  parameter  x  in  (168)  is 


x  = 


1 


£ 


K 

k=  1 


<7 


2 

skX 


(193) 


To  find  the  optimal  power  P  that  minimizes  the  total  power  of  the  system  under  the  constraints  that  the  SINR 
requirements  of  all  source-destination  pairs  arc  satisfied,  we  can  follow  the  procedure  described  in  the  previous 
section,  i.e. 


min  p  PY,k=iixcr2sk,r+l)  +  xL, 
lkdet(G[P) 


S.t.  P> 


det{  G(fc)) 


k  =  1,2,- ■■  ,K. 


(194) 
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We  note  that  the  objective  function  in  (194)  is  increasing  in  terms  of  increasing  P  >  0.  Thus,  the  optimal 
solution  P*  of  (194)  is  given  as 


p* 


rWlQ'%k=l,...,K 


det(G^)  J 


(195) 


Since  the  channel  quality  is  not  the  same,  in  general,  for  all  links,  it  is  implied  that  the  terms  ^  t[G(h)  ’  ^  = 
1, 2,  •  •  •  ,  AT  in  (194)  arc  not  equal  in  general,  which  means  that  the  equality  in  (194)  might  not  hold  for  all  k. 
As  a  result,  the  equal  power  allocation  strategy  generally  spends  more  power  than  what  is  needed. 


3.4.4  Optimum  Power  Assignment  Under  a  Total  Power  Budget  Constraint 

In  this  subsection,  we  will  apply  the  methodology  that  we  developed  in  Section  4  to  solve  a  related  max-min 
SINR  based  power  optimization  problem  for  multi-source  multi-destination  relay  networks.  Specifically,  we 
will  design  a  power  assignment  scheme  that  maximizes  the  minimum  SINR  among  all  source-destination  pairs 
subject  to  any  given  total  power  budget.  We  may  say  that  such  a  scheme  introduces  a  type  of  fairness  among 
different  source-destination  pairs.  The  optimization  problem  can  be  formulated  as  follows: 


{maxp1)...  ,pK-pr  {SINRfc}  , 

S  •  t  .  Ylk= t  Pk  +  Pr  <  Ptotal 

where  Ptotal  is  die  given  power  budget  of  the  system. 

Let  us  denote 

z  =  mink=1)...)K  {SINRfc}  , 

then  (196)  is  equivalent  to 


(196) 


(197) 


(  roaxpj,...  ,PK;Pr 


z, 

Ztl  Pk  +  Pr  <  Ptotal , 


SINRfc  >z,k=  1,  ,  K. 


(198) 


According  to  (174),  we  may  approximate  SINR/.  as 


SINRfc 


Pkdet{  G(fci) 

det{Gf) 


(199) 


where  the  matrices  G^  and  arc  specified  in  (169).  If  we  define  x  =  a2Pr,  then  for  any  given  x  >  0,  the 
transmission  power  at  the  relay  Pr  is  determined  as 


I< 

Pr  =  y  =  X  \  '  Pkf3sk ,r  +  xL. 

az  z ' 

fc=i 


(200) 


As  a  consequence,  we  can  reformulate  the  problem  in  (198)  as 


maxp1,...  ,pKz, 

s  •  t -Zk= l(XPsk,r  +  1  )Pk  +  XL  <  Ptotal, 


Pk>  z- 


det(  GP) 
det( GW)  ’ 


1  <  k  <  K. 


(201) 
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For  any  x  >  0,  we  denote  the  maximal  value  of  the  objective  function  in  (201)  as  z{x).  Then,  z[x)  must  satisfy 


det(G\ 
Pk  >z{x)-  k 


(■ k)\ 


Mk  =  1,2,-  ••  ,K. 


det{G^)’ 

Substituting  the  above  constraints  into  the  total  power  constraint  in  (201),  we  have 

K 


(202) 


_  det(  Gp]) 

Y^.XPsk,r  +  1  )z{x)  ■  de^G(fc^  +  xL  ^  Ptotal  ■ 


(203) 


Thus,  for  any  fixed  x  >  0,  the  maximal  value  of  the  objective  function  in  (201)  is 


z(x) 


Ptotal  xL 


Ylk=l(XPsk,r  P  1)  ‘ 


det(G^)  ' 

det(  G(fe)) 


(204) 


We  can  see  that 
parameter  (x  >  0). 
the  sources  is 


maximizing  z(x)  can  be  implemented  by  a  numerical  search  of  single  (one-dimension) 
If  we  denote  the  optimal  parameter  of  (204)  as  x* ,  then  the  optimum  power  assignment  to 


z(x*) 


det(  Gf) 

det(GW)’ 


k  =  1,2,--  -  , K, 


(205) 


and  the  corresponding  power  assigned  to  the  relay  is 


K 


P*  =  x '■ 


:(E  P*kPsk,r  +  L)- 


k= 1 


(206) 


We  note  that  z(x)  in  (204)  should  be  non-negative  for  any  given  x,  which  implies  that  Ptotal  —  xL  >  0,  i.e. 
x  <  Pto£al .  Thus,  we  only  need  to  search  the  single  variable  x  over  the  interval  [0,  Pto*al  ]  to  obtain  the  optimal 
solution  x*  that  maximizes  the  function  z(x). 

From  (199)  and  (205),  we  can  see  that  for  any  k  =  1,  2,  •  •  •  ,  I\ ,  SINK/.  =  z(x* ) ,  which  implies  that  the 
optimum  power  assignment  achieved  through  the  max-min  SINR  based  optimization  in  (196)  leads  to  the  same 
SINR  values  for  all  source-destination  pairs.  We  recall  that  for  the  total  power  minimization  problem  in  (167), 
the  optimal  power  assignment  ensures  that  the  resulting  SINR  for  each  source -destination  pair  is  equal  to  the 
SINR  requirement  of  the  corresponding  source-destination  pair,  i.e.  SINR/,.  =  7/.,  k  =  1,  •  •  •  ,  K.  It  is  natural 
to  expect  that  when  all  source-destination  SINR  requirements  7*  in  (167)  are  the  same,  the  power  assignment 
strategy  for  the  max-min  SINR  based  optimization  and  that  for  the  total  power  minimization  based  optimization 
would  be  the  same.  This  intuitive  interpretation  can  be  seen  from  the  derivations  in  (179)  and  (204).  We  note 
that  for  reasonably  high  power  budget  Ptotah  the  term  xL  in  (204)  can  be  ignored,  and  maximization  of  z[x) 
in  (204)  is  equivalent  to  minimization  of  the  dominator  which  is  related  to  the  objective  function  in  (179).  If  all 
the  source-destination  pairs’  SINR  requirements  in  (179)  arc  the  same,  i.e.  7/,.,  k  =  1,  •  •  •  ,  K,  then  the  optimal 
solutions  x*  of  z(x)  in  (204)  is  also  optimal  in  minimizing  v(x)  in  (179). 

For  comparison  purposes,  in  the  following  we  illustrate  the  max-min  SINR  problem  under  equal  power 
allocation.  The  optimization  problem  in  (196)  with  an  equal  power  assignment  can  be  expressed  as 


f  maxp  minfe=i)...;x{SINRfc}, 
l  s.t.  (K  +  1)P  <  Ptotai, 


(207) 


where  all  the  sources  and  the  relay  arc  assigned  the  same  transmission  power  P.  With  the  SINR  approximation 
in  (199)  and  z  =  mmi<fc<i<-{SINRfc},  we  have 


maxp 

s.t. 


D  Ptotal 
r  -  K+ 1  ' 


P>  Z- 


det(  Gf>) 


-  ’  det( G(fc>)’ 


k  =  1,  -  -  -  ,  K. 


(208) 
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Table  5:  Optimum  power  assignment  that  minimizes  the  total  power  under  given  SINR  constraints  using  the 
proposed  method  and  the  exhaustive  search  method  for  an  asymmetric  setting  (K  =  2) 


(7i.72)(dB) 

-Pi(dB) 

^2(dB) 

Pr(dB) 

Proposed  method 

(10,10) 

13.14 

19.89 

20.54 

Exhaustive  search 

(10,10) 

13.22 

19.91 

20.41 

Proposed  method 

(10,20) 

11.59 

27.34 

29.74 

Exhaustive  search 

(10,20) 

11.76 

27.28 

29.70 

Table  6:  Optimum  power  assignment  that  minimizes  the  total  power  under  given  SINR  constraints  using  the 
proposed  method  and  the  exhaustive  search  method  for  a  symmetric  setting  ( K  =  2) 


(7i>72)(dB) 

P 1  (dB) 

P2(dB) 

Pr(dB) 

Proposed  method 

(10,10) 

14.56 

14.56 

15.86 

Exhaustive  search 

(10,10) 

14.91 

14.91 

15.31 

Proposed  method 

(10,20) 

13.25 

23.25 

22.44 

Exhaustive  search 

(10,20) 

13.32 

23.32 

22.34 

Thus,  the  optimal  solution  is  P 


and  the  maximal  value  of  the  worst-case  SINR  is  given  by 


z  = 


Pfotal 

K  +  1 


f  det(  G(fc))  1 

maxi<-k<-K\l^f))- 


(209) 


We  can  see  that  in  the  equal  power  assignment  scheme,  the  minimum  SINR  value  among  all  source -destination 
pairs  is  directly  determined  by  the  given  total  power  budget  value. 


3.4.5  Numerical  Results  and  Discussions 

In  this  subsection,  we  perform  numerical  studies  to  illustrate  the  proposed  optimum  power  assignment  algo¬ 
rithms.  In  our  studies,  we  consider  a  slow  fading  scenario  and  then  approximate  \hi.3\2  by  o2- .  The  channel 
gain  for  each  channel  link  is  assumed  to  follow  a  path  loss  model,  where  the  variance  of  channel  coefficient 
is  given  by  a2rj  =  Sf  j'  ( i,j  £  {s}..  df.,  r } )  with  8,^  as  the  distance  of  the  channel  link  and  A  as  the  path-loss 
exponent  (A  =  3  in  our  numerical  studies). 

In  the  first  set  of  numerical  studies,  we  illustrate  the  optimum  power  assignment  that  minimizes  the  to¬ 
tal  power  consumption  under  the  condition  that  the  SINR  requirement  of  each  source-destination  pair  is  sat¬ 
isfied.  First,  we  consider  a  system  with  two  source-destination  pairs  and  one  relay,  i.e.  K  =  2,  and  the 
cross-correlation  of  the  two  source  codes  is  p  =  0.25.  We  study  two  sets  of  SINR  requirements  for  the  two 
source-destination  pairs:  (i)  [71,72]  =  [10, 10]dB,  and  (ii)  [71,72]  =  [10,  20]dB. 

Fig.  28  plots  the  total  power  consumption  with  varying  parameter  x  >  0  for  an  asymmetric  system,  where 
the  values  of  the  distance  between  nodes  are  set  as  follows:  5Sl,d1  =  2,  SSud2  =  3,  Ssitr  =  5S2^r  =  8S2td1  = 
Sr  cp  =  1, 8S2  d2  =  3  and  8rjj2  =  2.  The  optimal  values  of  the  parameter  that  minimize  the  total  power 
consumption  are  x*  =  0.63  and  x*  =  1.52  for  [71, 72]  =  [10, 10] dB  and  [10, 20]dB,  respectively.  Based  on  the 
optimal  value  x*,  we  obtain  the  corresponding  optimal  power  assignment  P\ ,  If  and  Pr  according  to  (178)  and 
(176),  as  listed  in  Table  5.  In  this  table,  we  also  compare  the  optimal  power  values  obtained  by  our  proposed 
approximation  method  and  those  obtained  by  exhaustive  search  based  on  the  optimization  in  (167).  We  observe 
that  the  optimal  power  values  obtained  by  the  two  methods  arc  almost  indistinguishable.  In  Fig.  29,  we  plot 
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the  total  power  consumption  with  varying  parameter  x  >  0  for  a  symmetric  system,  where  the  values  of  the 
distance  between  nodes  are  set  as:  <5si)dl  =  <5 s2,d2  =  ^si,d2  =  ^s2,di  =  2  and  6Sl:r  =  5S2}r  =  6r_dl  =  5T)d2  =  1- 
The  optimal  values  of  the  parameter  are  x*  =  0.32  and  x*  =  0.59  for  [71,72]  =  [10, 10]dB  and  [10,  20]dB, 
respectively.  Based  on  the  optimal  value  x*,  we  obtain  the  corresponding  optimal  power  allocation  P\ ,  /A  and 
Pr  again  according  to  (178)  and  (176),  as  listed  in  Table  6.  We  observe  that  when  the  SINR  requirements  of 
the  two  source-destination  pairs  arc  equal,  the  power  assignments  to  the  two  sources  arc  the  same.  However, 
when  the  SINR  requirements  arc  different  ([71, 72]  =  [10,  20]dB),  the  equal  power  assignment  is  not  optimum 
anymore.  We  can  see  in  Table  6  that  the  optimal  power  values  obtained  by  our  proposed  approximation  method 
are  indistinguishable  from  those  obtained  by  the  exhaustive  search  based  on  the  optimization  in  (167). 

We  now  repeat  our  studies  for  a  system  with  three  source-destination  pairs  and  one  relay,  i.e.  K  =  3,  where 
the  signatures  of  the  three  sources  are  orthogonal  to  each  other.  We  examine  both  an  asymmetric  case  (Fig. 
30)  and  a  symmetric  case  (Fig.  31).  In  both  figures,  we  consider  three  sets  of  SINR  requirements  for  the  three 
source-destination  pairs,  namely:  (i)  [71,72,73]  =  [20, 20, 20]dB,  (ii)  [71,72,73]  =  [10,  20, 20]dB,  and  (iii) 
[71, 72, 73]  =  [10, 10,  20]dB.  In  Fig.  30,  we  consider  an  asymmetric  system  where  the  distance  values  are  set 
as  SSktr  =  1  ,(k  =  1,2,3 )A)dl  =  lA,d2  =  2,Srtda  =  3  and  <5si)dl  =  2 ,SS2>d2  =  3  and  SS3jd3  =  4.  In  this 
case,  the  optimal  values  of  the  parameter  that  minimize  the  total  power  consumption  are  x*  =  2.34,  x*  =  2.77 
and  x*  =  3.17  for  the  three  sets  of  SINR  requirements,  respectively.  Based  on  the  optimal  values,  the  optimum 
power  assignments  Pi,  P2,  P3  and  Pr  can  be  obtained  accordingly  based  on  (178)  and  (176).  In  Fig.  31,  we 
consider  a  symmetric  system  with  SSk>r  =  <)rj.h  =  1  ,(k  =  1,2, 3)  and  5Sk}dk  =  2,  (A:  =  1,  2,  3).  The  optimal 
values  of  the  parameter  arc  x*  =  0.77,  x*  =  0.76  and  x*  =  0.75  for  the  three  sets  of  SINR  requirements, 
respectively.  We  can  see  that  the  optimal  values  x*  in  the  symmetric  system  arc  almost  the  same  for  the  three 
sets  of  SINR  requirements.  This  result  is  consistent  with  the  theoretical  discussion  that  the  optimal  value  x*  is 
independent  of  the  SINR  requirement  7^  but  it  depends  on  the  the  channel  condition. 

In  Table  7,  we  show  the  power  efficiency  of  the  proposed  optimum  power  assignment  scheme  by  listing  the 
resulting  total  power  consumption,  compared  to  the  total  power  consumption  that  results  from  the  equal  power 
assignment  scheme.  We  consider  both  an  asymmetric  case  (the  system  setup  is  the  same  as  that  in  Fig.  30)  and 
a  symmetric  case  (the  system  setup  is  the  same  as  that  in  Fig.  31).  We  can  see  that  in  the  asymmetric  system, 
the  power  savings  of  the  optimum  power  assignment  scheme  is  4-6dB  for  the  three  sets  of  SINR  requirements. 
The  more  unbalanced  the  SINR  requirements  of  the  three  source-destination  pairs  are.  the  more  power  savings 
of  the  optimum  power  assignment  scheme  compared  to  the  equal  power  assignment  scheme  can  be  achieved.  In 
the  symmetric  system,  the  power  savings  of  the  optimum  power  assignment  scheme  is  0.5-4. 5dB  for  the  three 
sets  of  SINR  requirements.  Comparing  the  results  between  the  asymmetric  system  and  the  symmetric  system, 
we  observe  that  the  optimum  power  assignment  scheme  gains  more  power  savings  in  the  asymmetric  system 
than  in  the  symmetric  system.  Also,  we  observe  that  the  total  power  consumption  of  the  equal  power  assignment 
scheme  is  the  same  in  each  of  the  asymmetric  and  symmetric  systems,  which  is  consistent  with  our  previous 
discussion  that  the  equal  power  assignment  depends  only  on  the  most  challenging/weakest  source -destination 
pair. 

In  the  second  set  of  numerical  studies,  we  illustrate  the  optimum  power  assignment  algorithm  that  maximizes 
the  minimum  SINR  among  all  source -destination  pairs  under  a  given  total  power  budget.  We  consider  a  system 
with  two  source-destination  pairs  and  one  relay,  i.e.  K  =  2,  with  a  total  power  budget  Ptotai  =  30dB.  We 
assume  that  the  cross-correlation  of  the  two  source  codes  is  p  =  0.25.  We  consider  both  an  asymmetric  case 
(the  system  setup  is  the  same  as  that  in  Fig.  28)  and  a  symmetric  case  (the  system  setup  is  the  same  as  that  in 
Fig.  29). 

Fig.  32  shows  the  maximization  of  the  minimum  SINR  with  varying  parameter  x  >  0  for  both  the  asymmetric 
and  symmetric  systems,  respectively.  The  optimal  values  of  the  parameter  that  maximize  the  minimum  SINR 
of  all  source-destination  pairs  are  x*  =  1.16  for  the  asymmetric  system  and  x*  =  0.69  for  the  symmetric 
system.  For  both  the  asymmetric  and  symmetric  cases,  the  optimum  power  assignments  P\ ,  P2  and  Pr  are 
determined  based  on  (205)  and  (206)  with  the  corresponding  optimal  value  x*,  listed  in  Table  8.  In  Table  8,  we 
also  compare  the  optimal  power  values  obtained  by  the  proposed  approximation  method  and  those  obtained  by 
exhaustive  search  based  on  the  optimization  in  (198).  We  observe  that  the  optimal  power  values  of  the  power 
assignment  obtained  by  the  two  methods  are  indistinguishable.  We  note  that  for  the  symmetric  case,  the  sources 
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Total  power  vs.  variable  x,  K=2 


Figure  28:  Minimization  of  the  total  power  consumption  under  varying  parameter  x  for  an  asymmetric  multi¬ 
source  multi-destination  relay  network  with  given  SINR  constraints  (K  =  2). 


Total  power  vs.  variable  x,  K=2 


Figure  29:  Minimization  of  the  total  power  consumption  with  varying  parameter  x  for  a  symmetric  multi-source 
multi-destination  relay  network  under  given  SINR  constraints  ( K  =  2). 
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Total  power  vs.  variable  x,  K=3 


Figure  30:  Minimization  of  the  total  power  consumption  with  varying  parameter  x  for  an  asymmetric  multi¬ 
source  multi-destination  relay  network  under  given  SINR  constraints  ( K  =  3). 


Total  power  vs.  variable  x,  K=3 


Figure  3 1 :  Minimization  of  the  total  power  consumption  with  varying  parameter  x  for  a  symmetric  multi-source 
multi-destination  relay  network  under  given  SINR  constraints  ( K  =  3). 
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The  worst  SINR  vs.  variable  x,  K=2,  P,  4  ,=30dB 

total 


Figure  32:  Maximization  of  the  minimum  SINR  among  all  source-destination  pairs  under  given  total  power 
budget  Ptotai  =  30dB  (K  =  2). 


The  max-min  SINR  vs.  total  power  budget,  K=2,  non-orthogonal 


Figure  33:  Comparison  of  the  minimum  SINR  resulting  from  the  proposed  optimum  power  assignment  scheme 
and  the  equal  power  assignment  scheme  with  any  given  total  power  budget  ( K  =  2). 
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Table  7 :  Comparison  of  the  total  power  consumption  resulting  from  the  optimum  power  assignment  scheme 
and  the  equal  power  assignment  scheme  (K  =  3) 


71,72,73  (dB) 

Total  Power  Consumption  (dB) 

Optimum 

Equal 

Asymmetric 

System 

[20,20,20] 

37.15 

41.58 

[10,20,20] 

36.82 

41.58 

[10,10,20] 

35.63 

41.58 

Symmetric 

System 

[20,20,20] 

29.79 

30.30 

[10,20,20] 

28.25 

30.30 

[10,10,20] 

25.85 

30.30 

Table  8:  Optimum  power  assignment  that  maximizes  the  minimum  SINR  among  two  source-destination  pairs 
under  given  total  power  budget  Ptotai  =  30dB,  by  using  the  proposed  method  and  the  exhaustive  search  method 
(K  =  2) 


Pt(dB) 

P2(dB) 

Pr(dB) 

Asymmetric 

System 

Proposed  method 

19.22 

25.34 

27.59 

Exhaustive  search 

19.37 

25.45 

27.66 

Symmetric 

System 

Proposed  method 

24.49 

24.49 

26.41 

Exhaustive  search 

24.64 

24.64 

26.50 

are  allocated  equal  power  under  the  max-min  SINR  optimization  scheme,  while  the  relay  utilizes  transmission 
power  which  is  some  what  different  from  that  allocated  to  the  sources. 

Finally,  in  Fig.  33,  we  compare  the  minimum  SINR  resulting  from  the  proposed  optimum  power  assignment 
scheme  and  the  minimum  SINR  resulting  from  the  equal  power  assignment  scheme.  We  consider  both  an 
asymmetric  case  and  a  symmetric  case,  and  the  system  setup  is  the  same  as  that  in  Fig.  32.  We  can  see  that 
in  the  asymmetric  case,  the  minimum  SINR  of  the  two  source-destination  pairs  is  improved  significantly  (by  at 
least  2dB)  when  we  use  the  proposed  optimum  power  assignment  scheme  instead  of  the  equal  power  assignment 
scheme.  In  the  symmetric  case,  we  can  see  that  the  performance  of  the  equal  power  assignment  scheme  is  very 
close  to  that  of  the  optimum  power  assignment  scheme  (actually,  in  this  case  all  the  source-destination  pairs 
have  the  same  performance). 

3.5  Cognitive  Code-Division  Channelization 

In  this  subsection,  first  we  briefly  describe  the  CR  CDMA  system  model  and  our  formulation  of  the  optimization 
problem.  Then,  we  present  in  detail  our  proposed  power  and  code-channel  allocation  solution.  The  performance 
of  the  proposed  scheme  is  evaluated  through  simulations  and  a  few  concluding  remarks  arc  drawn  at  the  end. 

3.5.1  System  Model  and  Problem  Formulation 

We  consider  a  primary  CDMA  system  with  processing  gain  (code  sequence  length)  L,  K  primary  transmitters 
PTi,  i  =  1, 2, ....  A',  and  a  primary  receiver  PR  (for  example,  K  uplink  transmissions  by  users  PT, ,  i  = 
1, 2, . . . ,  K,  to  base  station  PR).  We  also  consider  a  potential  concurrent  secondary  code-division  link  in  the 
spectrum  band  of  the  primary  system  between  a  secondary  transmitter  ST  and  receiver  SR  (Fig.  34).  All 
signals,  primary  and  secondary,  arc  supposed  to  propagate  over  flat-fading  channels  and  experience  additive 
white  Gaussian  noise.  We  denote  by  ht,  qt,  i  =  1,2,...,  K,  the  path  coefficients  from  PR  to  PR  and  SR, 
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Figure  34:  Primary /secondary  CDMA  system  model  of  K  primary  transmitters  FT, ,  i  = 
1,  2, . . . ,  K,  a  primary  receiver  PR,  and  a  secondary  transmitter-receiver  pair  ST,  SR.  All  paths 
h\, ,  fiK,  qi, ... ,  qk,  hs,  qs  exhibit  independent  (quasi-static)  Rayleigh  fading. 


respectively.  The  path  coefficients  from  ST  to  PR,  SR  are  denoted  by  hs  and  qs,  respectively.  All  path 
coefficients  are  modeled  as  Rayleigh  distributed  random  variables  that  arc  independent  across  user  signals  and 
remain  constant  during  several  symbol  intervals  (quasi-static  fading). 

After  carrier  demodulation,  chip-matched  filtering  and  sampling  at  the  chip  rate  over  the  duration  of  a  symbol 
(bit)  period  of  L  chips,  the  received  signal  at  the  primary  receiver  PR  can  be  represented  as 

K 

r  —  \  '  \J Ej hj s j bj  \J Eshsssbs  T  n p,  (210) 

i=l 

while  the  secondary  signal  received  by  SR  is 

K 

y  —  ^  '  \J Pj q> s t bj  -f-  \J Esqsssbs  -t-  ns  (211) 

i= 1 

where  Tj  >  0,  b,  £  {±1},  and  s,  G  VJJ,  \ s,  |  =  1,  denote  bit  energy,  information  bit,  and  normalized  signature 
vector  of  primary  user  i,  i  =  1, 2, . . . ,  K,  respectively;  Es  >  0,  bs  G  {±1},  and  ss  e  ML,  ||ss||  =  1,  denote  the 
bit  energy,  information  bit,  and  normalized  signature  vector,  respectively,  of  the  secondary  transmitter  ST;  np 
and  ns  represent  additive  white  Gaussian  noise  (AWGN)  at  PR  and  SR,  correspondingly,  independent  from 
each  other  with  0  mean  and  autoco variance  matrix  cr2I. 

The  linear  filters  at  the  primary  and  secondary  receivers  that  exhibit  maximum  output  SINR  [91]  can  be  found 
to  be 


W maxSINR,i  —  CiRp  Sj,  i  —  1,2,...,  K, 

WmaxSIN  R,s  —  *42 

where  Rp  =  T'{rrr},  Rs  =  E{yyT},  a,  C2  >  0  (T'l-}  denotes  statistical  expectation  and  T  is  the  transpose 
operator).  The  output  SINR  at  PR  with  respect  to  the  signal  transmitted  by  PT,L,  SINRi,  is  given  below 
followed  by  the  output  SINR  at  SR,  SINRs, 


SINRi 


_ E{\™LaxSINR,i(VElhibisi)\2} _ 

E{\'WmaxSINR,i('^2k=l,k^i  +  V^shs^sbs  +  np)H} 


EihfsJ  Rp/^Sj, 


(212) 
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SINRs 


_ E{\WmaxSINR,s('^Es<lsbsSs)\2} _ 

E{\WmaxSINR,s(^2k=l  V^kQkskbk  +  ns)h} 


Esql  sTR 


s/s^ 


(213) 


where  Rp/j  and  Rs/s  are  the  “exclude  i”  and  “exclude  s ”  input  data  autocorrelation  matrices  at  PR  and  SR, 
respectively,  defined  by 


K 

Rp/i  =  Ekhlsksk  +  Esh2sssJ  +  CJ2I, 

k=l,k^i 

I< 

Rs/«  =  Ek(ilSkSk  +  (j21- 

k=  1 


In  our  cognitive  radio  setup,  the  secondary  transmitter  has  to  guarantee  the  SINR  QoS  of  all  primary  users. 
In  this  spirit,  our  objective  is  to  find  the  transmission  bit  energy  Es  and  the  real-valued  normalized  signature 
vector  ss  that  maximize  SINRs  under  the  constraints  that  SINRi,  i  =  1,2, ... ,  K,  are  all  above  a  certain 
threshold  a  >  0,  i.e.  we  would  like  to  identify  the  optimal  pair 

(Es,  Sg)^  =  arg  max  Es sj RjA  ss 

Es>0,ss£M.L 

subject  to  Eihf  sf  R“^  s i>  a,  i  =  1,2, ,  K,  (214) 

ss  —  1,  Eg  '  Emax 


where  Emax  denotes  the  maximum  available/allowable  bit  energy  for  the  secondary  user. 

The  optimization  task  of  maximizing  a  quadratic  objective  function  (R“^  is  positive  definite)  subject  to  the 
constraints  in  (214)  is,  unfortunately,  a  non-convex  NP-hard  (in  L)  optimization  problem  [92].  In  the  following 
section,  we  delve  into  the  details  of  the  problem  and  derive  a  novel  realizable  suboptimum  solution. 


3.5.2  Proposed  Cognitive  Secondary  Channel  Design 


Using  the  matrix  inversion  lemma  [93]  on  Rn/’,  we  can  express  the  key  quadratic  constraint  expression 

sFRp/is* in  (214) as 


stR_1s-  — 

s*  p/i - 


p/i 

sjR-'si 


1  —  Eih‘fsfRp1Si  ’ 


i  =  1,2,...,  A:, 


(215) 


where  we  recall  that  Rp  =  E{rrr\  is  the  autocorrelation  matrix  of  the  whole  input  to  the  primary  receiver 
PR.  Then,  the  PR  SINR  constraints  in  (214)  become 


rT-f-~k  _ 1  OL 

S'  N  £  E^  +  aEM 


A 


=  7 h  i  =  l,2,...,K, 


(216) 


and  the  optimization  problem  can  be  rewritten  as 


(Es,  ss)opt  =  arg  max  Es s^R^Sg 


ss>o,sseM-L 


subject  to  sfR  1Sj  >7 z  =  1,2,...,  K, 
—  1,  Eg  < "  Emax. 


(217) 


Using  the  matrix  inversion  lemma  on  R  1  this  time,  we  see  that 


EghSR  /sgS^R  } 

■0-1  _  -0-1  _  5  p/s  b  S  p/g 

p  p's  1  +  Egh^sjR-^Sg 


(218) 
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where  Rp/S  is  the  autocorrelation  matrix  of  the  input  to  the  primary  receiver  PR  excluding  the  secondary 
transmission. 


K  K 

R-p/s  =  E{ \/ EihiSibi  T  np)(^^  \/ P'/Ji/Sjb j  T  rip)  } 
i=  1  i=l 

=  y  EjhjsjsJ  +  cr2I. 

i= 1 

Then,  inserting  (218)  in  (217)  we  can  express  the  optimization  constraints  as  explicit  functions  of  the  code 
sequence  of  the  secondary  user  ss,  i.e. 


£'s/r2s,R  }  sss^R  }  Si 

sfR-ts.a  P  +n,  i  = 


3,;  it  /  ^  - — - 

*  p/s  1  +  E.h  2s^R 


-1 . 


^p/s°a 

For  notational  simplicity,  define  the  L  x  L  matrix 


2t>  —  1 


B,  =  h^R^^iSi  Rn/jt  -  PihB R 


' p/s 


p/s 


where 


Pi  =  s f  Rp/jSi  -  7i,  *  =  1, 2, . . . ,  K. 


Then,  the  optimization  problem  in  (217)  can  be  rewritten  -for  one  more  time-  as 

xopt  =  arg  max  x1  RJ^  x 


subject  to  x1  B,x  —  Pi  <  0,  i  =  1, 2, . . . ,  K, 


x  x  <  Emax 


(219) 


(220) 

(221) 


(222) 


where  x  is  the  amplitude-including  transmitted  signature  vector  of  the  secondary  user,  x  =  \/Ws ss.  We  notice 
that  for  (222)  to  be  solved  at  the  secondary  transmitter  ST,  the  primary  receiver  PR  must  communicate  the 
matrix  parameters  B,  and  scalar's  p%,  i  =  1,2,...,  K.  Therefore,  no  explicit  communication  of  the  primary 
channel  codes  and  gains  is  required  that  may  directly  compromise  the  privacy/security  of  the  primary  system. 
In  terms  of  the  computational  effort,  however,  (i)  B;  ,  i  =  1,2,...,  K,  arc  not  necessarily  positive  semidefrnite, 
hence  the  problem  in  (222)  is  in  general  a  non-convex  quadratically  constrained  quadratic  program  (non-convex 
QCQP),  and  (ii)  the  complexity  of  a  solver  of  (222)  is  exponential  in  the  dimension  L  (NP-hard  problem). 

To  circumvent  these  two  difficulties,  we  first  observe  that  if  we  use  the  trace  property  of  matrices  U,  V, 
7>{UV}  =  7r(VU),  we  arc  able  to  represent  the  objective  function  in  (222)  as 

(223) 

where  X  =  xx7  .  Thus,  the  optimization  problem  in  (222)  takes  the  new  equivalent  matrix  form 

Xopt  =  arg  max  TriRT1  X} 

xe»ixi  s  s 

subject  to  7r{B,X}  <  Pi,  i  =  1,2, ,  K,  (224) 

7r{X}  <  Emax,  X  y  0,  rank(X)  =  1 

where  X  >z  0  denotes  that  the  matrix  X  is  positive  semidefinite. 

So  far,  we  have  shown  that  the  original  secondary  cognitive  link  design  problem  in  (214)  is  equivalent  to  the 
one  in  (217),  (222),  and  finally  (224),  and  is  non-convex  NP-hard.  To  effectively  attack  the  problem  anyway, 
we  now  propose  to  relax  the  rank  constraint  in  (224)  and  proceed  by  solving  the  following  problem  instead, 

X7  =  arg  max  7r{RjAX} 

xeKixi  s  s 

subject  to  7>{BjX}  <  Pi,  i  =  1,2, ... ,  K,  (225) 

7k{X}  <  Emax,  X^O. 
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Then,  (225)  is  a  convex  polynomial-complexity  problem  that  can  be  solved  using  semidefinite  programming. 
Strictly  speaking,  we  can  solve  (225)  in  polynomial  time  within  an  error  e  >  0  from  its  value  at  the  optimum 

point  X'.  More  specifically,  let  fQ  =  7>{R(/,1sX}jx=x'>  he.  fo  is  the  optimum  value  of  the  constrained 
(affine)  objective  function  in  (225).  Then,  for  any  given  e  >  0,  semidefinite  programming  guarantees  that 
we  can  converge  in  polynomial  time  (polynomial  in  the  input  size  L  and  in  the  error  requirement  function 
log  1/e)  to  a  solution  that  lies  in  ( fQ  —  e,  fQ)  [94].  In  this  paper,  for  the  semidefinite  programming  problem  in 
(225),  we  propose  to  use  a  primal-dual  interior-point  method  [95].  In  particular,  we  consider  the  problem  in 
(225)  as  the  primal  optimization  problem,  we  create  a  differently  parameterized  equivalent  dual  problem,  and 
then  solve  both  problems  iteratively  in  a  coupled  fashion.  Then,  each  iteration  can  be  implemented  in  0(L3) 
and  the  algorithm  converges  after  ()( L  log  1/e)  iterations  to  the  matrix  X"  that  makes  the  objective  function 
7>{R,(//X}  attain  a  value  within  (f0  —  e,  f0).  The  proposed  method  is  outlined  in  Fig.  35.  We  note  that  relaxing 
the  rank  constraint  of  the  non-convex  NP-hard  problem  in  (224)  we  created  the  convex  optimization  problem  in 
(225)  that  can  be  solved  in  0(L 4  log  1/e)  time  (by  semidefinite  programming  methods  as  described  in  Fig.  35). 
Of  course,  because  of  the  constraint  relaxation  itself  the  objective  function  evaluated  at  the  optimum  point  X' 
in  (225)  is  just  an  upper  bound  on  the  value  of  the  objective  function  evaluated  at  the  optimum  point  of  interest 
Xopt  of  (224),  7>{R//lsXopt}  <  Ir/R^X'}.  Moreover,  X'  is  not  available  exactly  either  and  instead  we  have 
X"  with  Tr{R-l  X"}  €  (7r{ R^X'}  -  e,  Tr{ RJ^X'}). 

To  summarize  our  developments  so  far  for  the  cognitive  design  of  a  code-division  secondary  link,  first,  for 
the  given  primary  SINR-QoS  threshold  a  >  0,  we  test  whether  Pi,  i  =  1,2,...  K,  in  (221)  are  all  greater 
than  zero.  If  this  is  not  true,  then  the  SINR-QoS  constraints  for  the  primary  users  cannot  be  met  and  outright 
no  secondary  transmission  is  allowed  (see  (low-chart  in  Fig.  36).  Otherwise,  we  run  the  procedure  of  Fig.  35 
which  returns  matrix  X".  If  the  rank  of  X"  is  1  with  eigenvalue,  eigenvector  pair  Ai,  ai,  then  we  already  have 
our  secondary  link  design  with  signature  ss  =  ai  and  transmission  amplitude  Es  =  Ai.  If  the  rank  of  X"  is  not 
1 ,  further  work  is  needed  as  described  below. 

When  X'  of  (225)  (or  in  practice  X"  returned  by  Fig.  35)  happens  to  be  of  rank  1  with  eigenvalue,  eigenvector 
pair  Ai,  ai,  then  X'  =  ~X°pt  in  (224)  and  x.opt  =  v/Aiai  in  (222).  Otherwise,  there  is  no  direct  path  from  X'  of 
(225)  to  x'W  in  (222).  In  this  case,  we  may  simply  consider  changing  the  search  for  an  optimal  vector  in  (222) 
to  a  search  for  an  optimal  probability  density  function  (pdf)  of  vectors  that  maximizes  the  average  objective 
function  subject  to  average  constraints,  i.e. 

/opi(x)  =  argmax  E{xtR;/1sx} 

subject  to  £{x7  Btx}  <  Pi,  i  =  1, 2, . . , ,  K,  (226) 

E{yi  x}  <  Emax 

where  /(x)  denotes  the  probability  density  function  of  x.  This  switch  to  a  statistical  optimization  problem  has 
been  known  as  the  “randomized  method”  in  semidefinite  programming  literature  [94].  Using  the  commutative 
property  between  trace  and  expectation  operators,  the  pdf  optimization  problem  in  (226)  takes  the  equivalent 
form 


/opt(x)  =  argmax  7>{R  A  U{xxT}} 

/(*)  ' 

subject  to  Tr{RiE{xx.7 }}  <  Pi,  i  =  1,2, ... ,  K,  (227) 

7>-{£{xxt}}  <  Emax. 

We  can  show  that  /op/(x)  is  in  fact  Gaussian  with  0  mean  and  covariance  matrix  X',  /op*(x)  =  2/(0.  X'). 
With  X"  from  Fig.  35  as  a  close  approximation  of  X',  we  can  draw  now  a  sequence  of  samples  xi,  X2, . . . ,  xp 
from  Ah(0,X").  We  test  all  of  them  for  “feasibility”  on  the  constraints  of  (222)  whether  x^B*xp  <  /%, 
Vi  =  1, 2, ... ,  K,  and  x^  xp  <  Emax,  p  =  1,2, ...  ,P,  and  among  the  feasible  vectors  (if  any)  we  choose  the 
one,  say  x^°\  with  maximum  x;R^'x  objective  function  value  (see  flow-chart  in  Fig.  36).  We  could  have 
suggested  at  this  time  a  cognitive  secondary  link  design  with  \/EsSs  =  x((/  Instead,  we  will  use  x':(l)  as  an 
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Formulate  the  pair  of  primal  and  dual  semidefinite  programming  problems: 


Primal 

Dual 

max  7r{R_/X} 

min  bTv 

veRx+l,ZeR^xi. 

subject  to  7r{B,;X}  <  Pi,  i  =  1,  2, . . 

.,K, 

subject  to  Y%=i  +  vK+1I  =  RJ^  +  Z 

r?-{X}  <  Emax,  X  A  0 

ZAO,  Vi>0,i=l,2,...,K+l 

where  b  =  [/3i,  /%,  •  •  •  >  Pk,  Emax\  and  Uj  is  the  ith  element  of  v. 

(i)  Initialize  (X,  v,  Z)  <—  (I,  0, 1). 

(ii)  Calculate  AX,  Av,  AZ  by  method  in  Appendix. 

(iii)  Replace  AX  <-  ±(AX  +  AXT). 

(iv)  Choose  0  <  r  <  1  and  select  steplengths  ni,  «2 

n\  =  min(l,r  •  ni),  h\  =  sup{fi  6  R  :  X  +  fiAX  A  0}, 
n2  =  min(l,  r  •  fi2),  fi2  =  sup{f2  G  R  :  Z  +  ij>AZ  A  0}. 

(v)  Update 

X  <—  X  +  ni  AX, 
v  <—  v  +  ngAv, 

Z  <—  Z  +  r^AZ. 

(vi)  Repeat  (ii)-(v)  until  X  feasible  solution  for  primal  problem,  (v,  Z)  feasible  solution  for  dual  problem, 

and  the  difference  of  primal  and  dual  objective  values  lies  within  (— e,  e). 

(vii)  Return  X"  <-  X. 


Figure  35:  Proposed  interior-point  algorithm  for  solving  (225). 
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s 


Input  a  >  0 


Figure  36:  Flow-chart  of  proposed  power  and  code  allocation  algorithm  for  secondary  link. 


initialization  point  to  an  iterative  procedure  that  will  lead  to  a  much  improved  link  design  vector.  The  iterative 
procedure  is  developed  below  and  its  performance  is  evaluated  by  simulation  studies  in  the  next  section. 

First  we  express  RS//.S  as 

Rs/s  =  SE  ST  +  a2I  (228) 

where  S  =  [si,  S2,  • . . ,  s^]  denotes  the  matrix  with  columns  the  signatures  of  the  primary  users,  and  S  = 
diag(E\q2,  E-tq2-  •  •  • ,  Ek :Qk)-  Using  the  matrix  inversion  lemma, 

r;A  =  41  -  4S(S_1  +  -4sts  )-1St.  (229) 

s/s  crz  cr4  a z 

Substitution  of  (229)  in  the  objective  function  of  (222)  leads  to 

xrR7/] Sx  =  ALxrt  x  _  J_XT  q  x  (230) 

S/S  aZ  <J4 

where  Q  =  S(E_1  +  \ S7  S)”1Si  .  In  (230),  the  first  term  4x7  x  is  a  convex  function  while  the  second  term 
-ArxT  Q  x  is  a  concave  function  (which  implies  that  ^x1  Q  x  is  convex).  Based  on  the  first-order  conditions 
of  convex  functions  [96] ,  we  have 

xTx  >  2x(0)tx  -  x(0)I  x(0)  (231) 

where  x'  ni  denotes  an  initial  feasible  vector.  Then,  we  combine  (230)  and  (231)  and  form  an  optimization 
problem  that  maximizes  the  following  concave  function 

2  X(0)TX  _  1  XTQ  x  _  J_x(0)Tx(0)  (232) 

aA  cr4 

that  leads  to  a  suboptimum  solution  for  our  original  problem  in  (222).  To  maximize  (232)  in  view  of  our 
constraints  in  (222),  we  restrict  all  non-convex  constraints  into  convex  sets  (linearization).  In  particular,  we 
consider  the  non-convex  constraints 


xT  B;x  -  A  <  0,  i  G  Xnc,  (233) 

where  Znc  denotes  the  set  of  all  indices  i  e  {1,2, . . . ,  A'}  for  which  xv  B,x  is  a  non-convex  function.  Then, 
we  decompose  the  matrix  B,  into  its  positive  and  negative  parts 

B i  =  B+  -  Br  (234) 

where  B^  =  h'2 R U s,sj R^,1,  and  Bj"  =  fak2 R”^  arc  all  positive  semidefinite.  Therefore,  the  original 
constraints  (233)  can  be  written  as 

xT B+  x  -  Pi  <  xT B^  x,  i  €  Inc,  (235) 

where  both  sides  of  the  inequality  arc  convex  quadratic  functions.  Linearization  of  the  right-hand  side  of  (235) 
around  the  vector  x((i  leads  to 


xtB+  x  -  ^  <  x(°)tB:  x(°)  +  2x(°)rB,7  (x  -  x^),  i  G  In 


(236) 


In  (236),  the  right-hand  side  is  an  affine  lower  bound  on  the  original  function  x7  B;  x.  It  is  thus  implied  that 
the  resulting  constraints  arc  convex  and  more  conservative  than  the  original  ones,  hence  the  feasible  set  of  the 
linearized  problem  is  a  convex  subset  of  the  original  feasible  set.  Thus,  by  linearizing  the  concave  parts  of  all 
constraints,  we  obtain  a  set  of  convex  constraints  that  are  tighter  than  the  original  non-convex  ones.  Now,  the 
original  optimization  problem  takes  the  form 


r(i)  - 


Ay(0 )Tv  _  _  J_x(0)T  (0) 

n2 


9-xvw,i  X - tX  Qx 

x  a4 


xv  7  =  argmax 

X 

subject  to  x7  B^  x  —  x(0^7  Bg  (2x  —  x^0^)  —  Pi  <  0,  i  G  In 


(237) 


X  B^x  (3i  <0,  i  G 


x  x  <  Emax 
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where  Inc  =  {1,2,...,  K)  —  Znc.  The  problem  in  (237)  is  a  convex  QCQP  problem  and  can  be  solved 
efficiently  by  standard  convex  system  solvers  [97]  to  produce  a  new  feasible  vector  x(1\  The  objective  function 
xtR{^x  in  (222)  evaluated  at  x1 1  j  takes  a  value  that  is  larger  than  or  equal  to  its  value  at  x((l\  Repeating 

iteratively  the  linearization  procedure,  we  can  obtain  a  sequence  of  feasible  vectors  x^°\  xn  x,2{  . . . ,  x:  7  : 
with  non-decreasing  values  of  the  objective  function  in  (222).  This  procedure  converges  after  very  few  (eight 
or  nine)  iterations  as  demonstrated  experimentally  in  the  following  section. 


3.5.3  Simulation  Studies  and  Discussions 

We  consider  a  primary  CDMA  system  with  signature  length  (system  processing  gain)  L  =  16  and  K  syn¬ 
chronous  users.  We  are  interested  in  establishing  a  secondary  code-division  transmitter/receiver  pair  when 
the  primary  system  is  fully  loaded  to  overloaded,  say  K  varies  from  16  to  20.  All  signatures  for  primary 
users  arc  generated  from  a  minimum  total -squared-correlation  optimal  binary  signature  set  which  achieves  the 
Karystinos-Pados  (KP)  bound  for  each  ( K,L )  pair  of  values4  [98]-[100].  The  transmission  SNRs  of  the  K 
primary  users  arc  all  set  equal  to  ^  =  15  dB,  i  =  1,  2, ...  K\  the  maximum  allowable  transmission  SNR  for 
the  secondary  link  is  set  to  Emgx  =  12 dB.  The  channel  coefficients  h%  and  </,,  i  =  1,2, ...  ,K,  (see  Fig.  34) 
are  taken  to  be  the  magnitude  of  independant  complex  Gaussian  random  variables  with  mean  0  and  variance 
4;  the  same  holds  true  for  hs  and  qs.  The  receiver  SINR  threshold  for  primary  users  is  set  to  a  =  2 dB  which 
corresponds  to  an  average  raw  bit-error-rate  (BER)  at  the  output  of  the  maximum  SINR  linear  filter  receiver 
of  about  1CT1.  Ten  thousand  (10,000)  system/secondary-line  optimization  experiments  arc  run  under  the  de¬ 
scribed  (quasi-static)  flat  fading  conditions.  When  random  vector  drawing  is  necessitated  by  the  flow-chart  in 
Fig.  36,  P  =  50  test  vector  points  arc  generated. 

In  Fig.  37,  we  plot  as  a  function  of  the  number  of  primary  system  users  K  the  percentage  of  time  that 
secondary  transmission  is  enabled  directly  under  the  case  rcinki X" )  =  1  or  by  the  iterative  linearized  opti¬ 
mizer  as  well  as  the  “Interference-Minimizing-Code-Assignment”  (IMCA)  scheme  in  [86].  We  observe  that 
significant  opportunity  exists  for  cognitive  secondary  transmission  when  the  primary  system  is  fully  loaded 
(K  =  L  =  16).  As  we  expect.  Fig.  37  shows  that  the  frequency  of  secondary  transmissions  reduces  as  the  pri¬ 
mary  system  load  increases.  We  observe  also  that  our  proposed  scheme  offers  more  opportunities  for  cognitive 
secondary  transmission  than  [86]. 

In  Fig.  38,  we  test  the  quality  of  the  secondary  transmission  line  (the  pre-set  SINR-QoS  of  the  primary 
system  is  -of  course-  guaranteed  by  the  algorithmic  procedure)  and  the  significance  of  the  iterative  linearized 
optimizer  in  (237)  (see  flow-chart  in  Fig.  36).  We  fix  the  primary  system  load  K  =  16  (fully  loaded)  and 
plot  the  secondary  receiver  average  SINR  for  the  experimental  instants  of  rank(X. ")  >  1  as  a  function  of  the 
iteration  of  the  optimizer  initialized  at  the  point/design  x(0)  that  is  the  best  out  of  P  =  50  samples  drawn  from 
the  7V(0,  X")  pdf.  It  is  pleasing  to  observe  that  eight  or  nine  iterations  are  enough  for  effective  convergence. 

Finally,  in  Fig.  39,  to  gain  visual  insight  into  the  operation  of  the  primary/secondary  system  we  plot  the 
instantaneous  receiver  SINR  of  a  primary  signal  and  the  secondary  signal  for  the  case  I\  =  17  over  an  exper¬ 
imental  data  record  sequence  of  1000  Rayleigh  fading  channel  realizations.  Missing  secondary  signal  SINR 
values  indicate  the  instances  when  no  secondary  transmission  was  allowed.  The  proposed  scheme  almost  dou¬ 
bled  the  occurrences  of  secondary  transmission  compared  to  [86].  When  secondary  transmissions  do  occur  for 
both  schemes,  the  joint  power  and  sequence  optimization  executed  by  the  proposed  scheme  results  in  superior 
SINR  performance  for  the  secondary  receiver  over  [86]. 

Appendix:  Calculation  of  AX,  Av,  AZ  of  Fig.  35 

Let  vec{-}  denote  column-by-column  matrix  vectorization  and  mat{-}  the  exact  inverse  operation.  Choose 
0  <  6  <  1  and  define  fx  =  dT'^Z\ 


4For  L  =  16,  when  K  <  L  the  KP-optimal  sequences  coincide  with  the  familiar  Walsh-Hadamard  signature  codes. 
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Figure  37:  Secondary  transmission  percentage  as  a  function  of  the  number  of  primary  users  under 
Cases  rank(X.")  =  1  and  >  1  (the  study  includes  also  the  code  assignment  scheme  in  [86]). 


Figure  38:  Secondary  receiver  SINR  as  a  function  of  the  iteration  step  of  the  linearized  optimizer 
in  (237)  initialized  at  the  best  feasible  sample  out  of  P  =  50  drawings  from  the  jV(0,  X")  pdf. 
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Figure  39:  Instantaneous  output  SINR  of  a  primary  signal  against  SINR-QoS  threshold  a  (thick 
line)  and  instantaneous  output  SINR  of  secondary  signal. 


In  Fig.  35,  AX,  Av,  AZ  are  obtained  by  solving  the  following  linear  system 


where 


'  0  Bt  I 

'  nec{AX}  ' 

nec{Ti} 

BOO 

Av 

= 

f  2 

E  0  F 

nec{AZ} 

_  nec{T3}  _ 

B(A'+1)xL2  — 


(nec{Bi})T 

(vec{BK})T 

(nec{I})T 


(238) 


TiLxl  =  Rs/1s  +  Z-mat{BTv}’t2(Jf+i)xi  =  b  -  Buec{X},  T3lxL  =  y.  I  -  XZ,  EL2xL2  =  Z(g)I,  and 
^L2xL2  =  I  X  (0  denotes  the  standard  Kronecker  product).  Applying  Gauss  elimination,  the  solution  is 


Av  =  (BE_1FB)_1(t2  +  BE_1(Fnec{Ti}  -  nec{T3})),  (239) 

nec{AX}  =  -E~1(F(nec{Ti}  -  BTAv)  -  nec{T3}),  (240) 

nec{AZ}  =  nec{Ti}  —  BrAv.  (241) 


4  Conclusions 

In  this  report,  we  summarized  our  findings  resulting  from  the  project  and  described  in  details  the  system  models 
and  the  proposed  methods  and  procedures.  First,  we  developed,  for  the  first  time,  a  closed-form  asymptotically 
tight  (as  SNR  — >  oo)  approximation  of  the  outage  probability  for  the  DF  cooperative  ARQ  relay  scheme  under 
fast  fading  conditions.  The  analytical  approach  relies  on  two  lemmas  that  we  established  in  this  project.  The 
closed-from  expression  provides  significant  insight  on  the  benefits  of  the  DF  cooperative  ARQ  relaying  relative 
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to  direct  ARQ  scheme  in  fast  fading  scenarios  and  shows  that  the  cooperative  scheme  achieves  diversity  order 
equal  to  2L  —  1  while  the  diversity  order  of  the  direct  scheme  is  only  L.  Simulation  and  numerical  results  show 
that  the  closed-form  approximation  of  the  outage  probability  is  tight  at  high  SNR.  Based  on  the  asymptotically 
tight  approximation  of  the  outage  analysis,  we  were  able  to  determine  an  optimum  power  allocation  for  the 
DF  cooperative  ARQ  relay  scheme.  It  turns  out  that  the  equal  power  allocation  is  not  optimum  in  general 
and  the  optimum  power  allocation  depends  on  the  link  quality  of  the  channels  related  to  the  relay.  It  shows 
that  we  should  allocation  more  power  at  the  source  and  less  at  the  relay.  Simulation  results  show  that  the  DF 
cooperative  ARQ  relay  scheme  with  the  optimum  power  allocation  has  a  performance  improvement  of  more 
than  ldB  compared  to  the  scheme  with  the  equal  power  allocation. 

Second,  we  determined  an  optimal  transmission  power  assignment  strategy  for  the  H-ARQ  protocol  in  quasi¬ 
static  Rayleigh  fading  channels  to  minimize  the  average  total  transmission  power.  A  locally  optimal  power 
sequence  was  found  first  by  the  Lagrangian  multiplier  and  then  it  was  shown  optimal  globally.  The  optimal 
transmission  power  sequence  is  described  by  a  set  of  equations  which  allow  an  exact  recursive  calculation  of 
the  optimal  power  sequence.  To  reduce  calculation  complexity,  we  also  developed  an  approximation  to  the 
optimal  power  sequence  that  is  close  to  the  numerically  calculated  exact  result.  It  is  interesting  to  observe  that 
the  optimal  transmission  power  sequence  is  neither  increasing  nor  decreasing;  its  form  depends  on  given  total 
power  budget  and  targeted  outage  performance  levels.  The  optimal  power  assignment  sequence  reveals  that 
conventional  equal-power  assignment  is  far  from  optimal.  For  example,  for  a  targeted  outage  performance  of 
Hr  5  and  maximum  number  of  transmissions  L  =  5,  the  average  total  transmission  power  by  the  optimum 
assignment  is  about  27  dB  less  than  that  of  using  the  equal-power  assignment.  We  also  observe  that  with 
the  same  targeted  outage  performance,  the  larger  the  number  of  retransmission  rounds  allowed  in  the  H-ARQ 
protocol,  the  larger  the  performance  gain  between  the  optimal  power  assignment  scheme  and  the  equal-power 
assignment  scheme.  Moreover,  the  lower  the  required  outage  probability,  the  more  performance  gain  of  the 
optimal  power  assignment  strategy  compared  to  the  equal-power  strategy. 

Third,  we  analyzed  the  outage  probability  performance  of  the  differential  ampl i  fy-and- forward  (DAF)  relay¬ 
ing  and  obtained  an  exact  formula  that  involves  only  a  single  integral  and  is  much  more  practical  than  the  only 
previously  known  result  that  involves  a  triple  integral.  To  gain  further  insight  in  the  process  of  the  DAF  relaying, 
we  also  developed  a  simple  closed-form  approximation  for  the  outage  probability  which  is  tight  at  high  SNR. 
Based  on  the  tight  outage  probability  approximation,  we  finally  determined  an  asymptotically  optimum  power 
allocation  scheme  for  the  DAF  relaying.  We  have  carried  out  extensive  numerical  calculations  and  simulations 
to  validate  and  illustrate  our  analysis. 

Fourth,  we  designed  and  optimized  a  multi-source  multi-destination  relay  network  where  a  relay  amplifies 
and  forwards  simultaneously  the  signals  received  from  all  sources.  We  developed  two  optimum  power  assign¬ 
ment  schemes.  The  first  scheme  minimizes  the  total  power  consumption  of  all  sources  and  the  relay  under  the 
constraint  that  the  SINR  requirement  of  each  source-destination  pair  is  satisfied,  while  the  second  scheme  max¬ 
imizes  the  minimum  SINR  of  all  source-destination  pairs.  Clearly,  both  optimization  problems  as  stated  above 
involve  K  power  variables,  where  K  is  the  number  of  source-destination  pairs  in  the  network,  which  implies 
that  an  exhaustive  search  approach  is  prohibitive  for  large  K.  In  this  paper,  we  derived  an  asymptotically  tight 
approximation  of  the  SINR  that  allows  us  to  reformulate  the  original  optimization  problems,  and  eventually 
reduce  them  to  single-parameter  optimization  problems,  which  can  be  easily  solved  by  numerical  search  of  the 
single  parameter.  Then,  the  corresponding  optimal  transmission  power  at  each  source  and  at  the  relay  can  be 
calculated  directly.  The  proposed  optimization  scheme  is  scalable  and  the  power  assignment  algorithm  has  the 
same  optimization  complexity  for  any  number  of  source-destination  pairs  in  the  network.  Moreover,  we  also 
applied  the  developed  methodology  to  solve  a  related  max-min  SINR  based  optimization  problem,  where  for 
any  given  total  power  budget,  we  arc  able  to  determine  the  optimum  power  assignment  for  the  sources  and  the 
relay  to  maximize  the  minimum  SINR  of  all  source-destination  pairs. 

Finally,  we  considered  the  problem  of  cognitive  code-division  channelization  for  cognitive  airborne  net¬ 
works.  The  goal  is  to  improve  the  efficiency  of  spectrum  utilization  and  coexisting  of  primary  and  secondary 
users  in  heterogeneous  cognitive  airborne  networks,  by  jointly  assigning  power  and  code-channel  allocation 
for  a  secondary  transmitter/receiver  pair  coexisting  with  a  primary  CDMA  system.  Unfortunately,  in  a  com¬ 
mon  Rayleigh  fading  wireless  environment,  the  formulated  constrained  optimization  problem  is  non-convex 
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and  NP-hard  in  the  code  vector  dimension.  Nevertheless,  in  pursuit  of  a  computationally  manageable  and 
performance-wise  appealing  suboptimal  solution,  we  first  converted  the  amplitude/code-vector  optimization 
problem  to  an  equivalent  matrix  optimization  problem  under  a  rank-1  constraint.  Disregarding  (relaxing  in 
formal  language)  the  rank- 1  constraint  makes  the  problem  amenable  to  an  easy  polynomial-cost  semidefinite 
programming  solution.  When  luckily,  a  rank- 1  matrix  happens  to  be  returned,  optimal  secondary-line  design  is 
achieved.  For  the  common  case  of  a  higher  rank,  we  developed  an  iterative  linearized  polynomial-cost  convex 
optimizer  with  much  appealing  (yet  suboptimal)  amplitude/code-vector  design  solutions  after  a  few  iterations. 
Extensive  numerical  studies  have  validated  our  theoretical  developments  and  the  proposed  iterative  algorithm. 
The  proposed  scheme  almost  doubled  the  occurrences  of  secondary  transmission  compared  to  some  early  work. 
The  joint  power  and  sequence  optimization  executed  by  the  proposed  scheme  results  in  superior  SINR  perfor¬ 
mance  for  the  secondary  receiver  compared  to  prior  work.  The  proposed  scheme  has  great  potential  to  improve 
the  efficiency  of  coexisting  of  primary  and  secondary  users  in  heterogeneous  cognitive  airborne  networks. 
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ACRONYMS 


ACK 

ACKnowIedgement 

AF 

Ampl  i  fy-and-Forward 

ARQ 

Automatic-Repeat-reQuest 

CDMA 

Code-Division  Multiple-Access 

CR 

Cognitive  Radio 

CSI 

Channel  State  Information 

DAF 

Differential  Amplify-and-Forward 

DDF 

Differential  Decode-and-Forward 

DF 

Decode-and- Forward 

H-ARQ 

Hybrid  Automatic-Repeat-reQuest 

IRC 

Interference  Relay  Channel 

MIMO 

Multiple-Input-Multiple-Output 

MRC 

Maximal-Ratio-Combining 

NACK 

Negative-ACKnowledgement 

QoS 

Quality-of-Service 

SINR 

Signal-to-Interference-plus-Noise  Ratio 

SNR 

Signal-to-Noise  Ratio 
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